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EXECUTIVE  SUMMARY 


This  report  describes  model  development  and  application 
activities  undertaken  in  the  investigation  of  the  sensiti¬ 
vity  of  remotely  sensed  radiances  in  the  DMSP  VHR/LF  spectral 
region  (0.4  to  1.1 pm)  to  near  surface  meteorological  para¬ 
meters  such  as  relative  humidity,  wind  speed,  and  aerosol 
characteristics  including  mixing  height  and  number  density. 

A  research  grade  computer  code  developed  in  an  earlier 
study  (Isaacs,  1980)  to  simulate  DMSP  visible  and  near  in¬ 
frared  radiances  has  been  upgraded  and  augmented  with  addi¬ 
tional  capabilities  desirable  in  imagery  analysis  for  tacti¬ 
cal  purposes  including:  (a)  treatment  of  continental  (i.e., 
urban  and  rural)  aerosol  size  distribution  and  composition 
models,  (b)  data  simulation  calculations  based  on  realistic 
sun/sensor  geometries,  and  (c)  approximate  ocean  surface  re¬ 
flection  modeling  for  potential  sunglint  regions.  Addition¬ 
ally,  the  approximate,  analytical  radiative  transfer  submodel 
incorporated  into  the  code  is  fully  described  and  the  results 
of  comparisons  with  accurate  numerical  treatments  are  provi¬ 
ded  for  a  variety  of  scattering  situations.  Comprehensive 
program  documentation  is  included  consisting  of  listings  of 
the  algorithm,  sample  output,  and  operating  instructions  for 
implementation  on  a  desk  top  minicomputer. 

A  variety  of  simulation  studies  employing  the  algorithm 
are  discussed  illustrating  the  effect  of  relative  humidity 
and  mixing  height  on  radiances  in  non  sunglint  situations 
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for  both  rural  and  urban  aerosol  models  as  compared  to  the 
maritime  models  previously  examined  (i.e.  in  Isaacs,  1980). 
Radiance  calculations  in  sunglint  regions  using  the  maritime 
aerosol  model  illustrate  the  effect  of  relative  humidity, 
wind  speed,  and  aerosol  scale  height  on  the  intensity  gra¬ 
dient  observed  along  simulated  sensor  scan  lines  in  the 
vicinity  of  the  primary  specular  point.  Additionally,  the 
model  is  exercised  to  simulate  a  swath  of  DMSP  imagery  for 
illustrative  purposes.  Results  of  these  simulations  con¬ 
firm  the  behavior  generally  observed  in  relevant  imagery 
for  corresponding  meteorological  analysis  case  studies  as 
previously  reported  (cf .  Fett  and  Mitchell,  1973) . 

Comparison  of  model  predicted  values  to  measured  data 
are  presented  in  two  contexts:  (a)  multiwavelength  aerosol 
extinction  coefficient  profiles  calculated  from  field  mea¬ 
sured  size  distributions  are  compared  to  model  predicted 
values  based  on  measured  relative  humidity  profiles  as  the 
required  meteorological  input,  and  (b)  model  calculated  ra¬ 
diances  along  a  simulated  scan  line  in  sunglint  are  compared 
to  data  points  obtained  from  sensometry  for  the  correspond¬ 
ing  hardcopy  imagery. 

In  addition  to  these  imagery  simulation  exercises  po¬ 
tentially  related  to  meteorological  analyses,  application 
of  the  algorithm  to  the  inference  of  surface  propagation 
parameters  is  explored.  Utilizing  the  model's  capability 
to  provide  simulated  radiances  for  specified  meteorology 
and  choice  of  aerosol,  DMSP  radiance/aerosol  optical  depth 
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relationships  are  computed.  Results  provide  a  theoretical 
fit  of  radiance  to  aerosol  optical  depth  for  maritime,  rural, 
and  urban  aerosols  and  relative  humidities  of  50  through 
95%.  These  may  potentially  be  applied  to  the  retrieval  of 
aerosol  optical  depth  (at  visible  wavelengths)  from  appro¬ 
priately  calibrated  DMSP  radiances . 

Since  surface  propagation  parameters  are  desired  at 
other  than  visible  wavelengths,  a  prototype  algorithm  is 
devised  to  infer  optical  range  multispectrally  using  the 
DMSP  radiance/optical  depth  relationships  described  above 
to  provide  the  aerosol  contribution  to  extinction. 

The  following  conclusions  and  observations  are  based 
on  the  work  summarized  above: 

•  Although  limited  in  scope,  the  validation  effort 
associated  with  the  MAGAT  data  set  indicates 
that  the  AFGL  maritime  aerosol  model  performs 
relatively  well  where  expected. 

•  A  variety  of  comparisons  performed  between  the 
model  adopted  approximate  radiative  transfer 
scheme  and  accurate  numerical  treatments  for 

a  range  of  optical  depths,  angular  scattering 
functions,  and  observer/sun  orientations  leads 
to  the  conclusion  that  simulated  radiances  are 
accurate  to  within  10%  overall.  This  level  of 
accuracy  is  consistent  with  the  accuracy  of  input 
data  and  the  intended  purpose  of  the  model,  i.e. 
the  simulation  of  meteorological  influences  in 


the  imagery.  The  results  of  specific  sensitivity 
tests  conclude  that  adoption  of  simple  analytical 
phase  function  representations  and  the  hybrid 
modified  6-Eddington  flux  parameterization  are 
appropriate  choices. 

•  The  wind -driven,  rough  ocean  surface  reflection 
model  implemented  within  the  radiance  simulation 
code  provides  a  good  representation  of  the  sun- 
glint  phenomenon  as  compared  to  accurate  numeri¬ 
cal  treatments,  simulated  radiance  comparisons 

to  hardcopy  imagery,  and  qualitative  observations. 

•  Model  simulations  in  sunglint  regions  suggest  a 
decrease  in  radiance  in  the  vicinity  of  the 
primary  specular  point  due  to  increases  in 
aerosol  optical  depth  (related  to  increased 
relative  humidity  or  scale  height).  Away  from 
the  primary  sunglint  region,  radiances  increase 
with  increasing  optical  depth  as  predicted  in 
previous  studies. 

•  A  comparison  of  the  optical  properties  of  the 
maritime  and  recently  implemented  rural  and 
urban  models  suggests  that  some  care  must  be 
taken  in  the  interpretation  of  meteorological 
influences  on  individual  air  masses.  For  example, 
the  maritime  aerosol  coarse  mode  is  much  more 
hygroscopic  than  the  others,  while  the  carbon 
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component  of  the  urban  aerosol  is  critical  in 
determining  its  radiative  properties. 

•  As  a  consequence  of  the  above,  each  aerosol 
composition  model  and  relative  humidity  range 
exhibit  unique  radiance/optical  depth  relation¬ 
ships.  While  those  for  maritime  and  rural  aero¬ 
sols  are  close  enough  to  permit  a  single  linear 
predictor  equation  to  be  used  operationally,  those 
for  the  urban  aerosol  introduce  a  potential  ambi¬ 
guity.  It  appears  that  it  may  be  important  to 
know  whether  the  air  mass  aerosol  under  investi¬ 
gation  has  any  carbon  component  present  (anthro¬ 
pogenic  or  natural) .  This  may  be  important  in 
the  vicinity  of  inland  seas,  coastal  regions,  or 
open  ocean  in  situations  of  long  range  transport. 

A  parameter  which  may  potentially  be  retrievable 
in  this  regard  is  the  single  scattering  albedo 

(Kaufman  and  Joseph,  1982) . 

* 

•  A  sample  multispectral  range  retrieval  exercise 
illustrates  both  the  dominance  of  aerosol  extinc¬ 
tion  in  constraining  range  at  visible  wavelengths 
and  the  complementary  significance  of  non-aerosol 
extinction  contributions  in  the  infrared.  With 
respect  to  the  latter,  calculated  aerosol  extinc¬ 
tion  accounts  for  less  than  7%  of  the  total  at 
10.6um,  while  that  for  water  vapor  continuum  ab¬ 
sorption  is  75%.  While  visible  radiance  data  may 
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1 .  INTRODUCTION 


1 . 1  Overview 

This  report  describes  the  results  of  applied  research 
efforts  conducted  at  Atmospheric  and  Environmental  Research, 
Inc.  (AER)  in  support  of  the  Naval  Environmental  Prediction 
Research  Facility  (NEPRF)  during  the  period  1  October  1980 
to  30  September  1982.  In  a  previous  study  sponsored  by  NEPRF 
AER  investigators  undertook  a  study  of  the  influence  of  met¬ 
eorological  variables  and  optical  propagation  parameters  on 
remotely  sensed  data  in  the  visible  ( 0.4-0. 7ym)  and  near  in¬ 
frared  (0.7-l.lym)  spectral  regions  such  as  that  potentially 
available  from  the  DMSP  sun  synchronous  polar  orbiting 
satellite  sensor  platforms.  In  the  course  of  that  effort 
(Isaacs,  1980)  a  computationally  efficient  model  was  devel¬ 
oped  to  simulate  the  sensitivity  of  wavelength  dependent, 
satellite  incident  radiances  to  variations  in  characteristic 
meteorological  variables  such  as  relative  humidity,  wind 
speed,  mixing  height,  and  aerosol  number  density  and  size 
range.  The  model  is  based  in  part  on  extant  propagation 
codes  commonly  applied  within  the  DOD  community  to  evaluate 
point-to-point  optical  transmission  and,  in  that  context, 
provides  an  extension  of  these  codes  to  simulate  satellite- 
based  observations.  The  model  developed  may  therefore  pro¬ 
vide  an  informative  link  between  concurrent  in  situ  and 
remotely  sensed  data  bases. 


Although  simulation  results  presented  provided  insights 
into  a  variety  of  imagery  analysis  questions,  a  number  of 
limitations  to  the  model  were  identified.  The  objective  of 
the  work  reported  here  was  to  supplement  and  extend  the  theo 
retical  studies  previously  performed  and  provide  the  founda¬ 
tion  for  an  operationally  useful  analysis  tool  by  validating 
elements  of  the  model. 

1.2  Study  Objectives 

The  objectives  of  this  study  included:  (a)  validation 
of  model  and  submodel  elements  through  comparison  of  model 
predicted  results  based  on  in  situ  meteorological  input  data 
and  field  measured  quantities  where  possible,  (b)  incorpora¬ 
tion  of  non-maritime  aerosol  optical  properties  models  to 
facilitate  treatment  of  coastal  regions  and  areas  over  the 
open  ocean  subject  to  continental  aerosol  sources,  (c)  im¬ 
plementation  of  the  simulation  algorithm  in  a  format  more 
consistent  with  the  operational  analysis  of  imagery,  (d)  ex¬ 
tension  of  the  domain  of  applicable  cases  to  include  those 
in  potential  sunglint  areas,  (e)  application  of  the  model 
to  the  retrieval  of  surface  propagation  parameters  particu¬ 
larly  aerosol  contributions  to  range  degradation,  and 
(f)  provision  for  program  documentation  and  description. 

1.3  Report  Summary 

This  document  is  organized  into  ten  sections.  Section 
2  briefly  reviews  work  previously  completed.  Section  3 


discusses  model  validation  studies  based  on  comparison  of 
measured  and  model  predicted  aerosol  extinction.  Section 
4  describes  the  incorporation  of  continental  (i.e.,  urban 
and  rural)  aerosol  optical  properties  submodels  within  the 
simulation  code.  Section  5  catalogs  the  equations  used  in 
the  model  to  provide  the  necessary  sun/sensor  geometry  for 
the  simulation  of  real  data.  Section  6  discusses  the  radi¬ 
ative  transfer  theory  implemented  within  the  model.  Section 
7  documents  the  treatment  of  ocean  surface  reflection  in 
potential  sunglint  regions.  Section  8  provides  theoretical 
DMSP  radiance/aerosol  optical  depth  relationships  calculated 
by  the  model  which  may  be  used  to  infer  the  aerosol  contri¬ 
bution  to  range  degradation.  Section  9  proposes  an  algorithm 
for  multispectral  range  retrieval  based  on  the  availability 
of  satellite  visible  radiance  data  (using  the  results  of 
Section  8) .  Finally,  Section  10  summarizes  relevant  conclu¬ 
sions  and  recommendations. 

Two  supporting  appendices  are  included.  Appendix  A 
includes  a  derivation  of  the  radiative  transfer  solutions 


used  to  calculate  radiance  and  Appendix  B  incorporates  a 
program  description  and  basic  users'  guide  into  the  report. 
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2. 


REVIEW  OF  PREVIOUS  WORK 


2.1  Simulation  Model  Development 

Mathematical  models  to  simulate  satellite  based  data 
sets  have  been  applied  extensively  in  the  recent  past  both 
a  priori  to  investigate  the  optimum  configuration  of  pro¬ 
posed  sensor  systems  for  observing  specific  phenomona  or 
parameters  and  after  system  deployment  as  tools  to  aid  in 
the  analysis  and  reduction  of  resultant  real  data.  A  sig¬ 
nificant  amount  of  input  data  describing  the  composition 
and  optical  properties  of  the  atmosphere-surface  system 
is  required  to  successfully  implement  these  models.  Mete¬ 
orological  imaging  sensors  (as  contrasted  with  sounders) 
are  generally  designed  to  sense  near  surface  phenomena 
and,  therefore,  simulation  of  such  "window  region"  data 
requires  input  of  meteorologically  dependent  optical  pro¬ 
perties  of  the  lower  atmosphere.  A  number  of  previous 
studies  concerning,  for  example,  the  near  infrared  proper¬ 
ties  of  clouds,  precipitation  monitoring,  and  remote  sensing 
of  surface  properties  (Isaacs  et  al. ,  1974;  Gaut  et  al. ,  1975 
Isaacs  and  Chang,  1975)  employed  standard  optical  propagation 
models  linked  with  radiative  transfer  theory  to  provide  such 
simulations . 

Studies  of  so-called  anomalous  gray  shade  patterns  (see 
Fett  and  Mitchell,  1977),  identifiable  in  wide  band  (0. 4-1.1 
um)  DMSP  VHR  and  LF  imagery  suggested  the  role  of  low  level 
haze  and  moisture  as  one  potential  cause,  highlighting  the 
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role  of  meteorological  variables  such  as  relative  humidity 
in  determining  remotely  sensed  imagery  data.  In  a  subsequent 
study  undertaken  to  investigate  the  relative  wavelength  depen¬ 
dent  effects  of  atmospheric  particulates  and  water  vapor  con¬ 
tent  in  these  gray  shade  phenomena  (Barnes  et  al.,  1979),  a 
simulation  model  was  employed  using  optical  properties  based 
on  LOWTRAN  2  (Selby  and  McClatchey,  1972)  and  a  radiative 
transfer  code.  Results  suggested  the  significant  role  of 
atmospheric  aerosol  content  in  backscattering  solar  radiation. 
A  related  first  order  analysis  (Fett  and  Isaacs,  1979)  empha¬ 
sized  the  role  of  particulates  and  their  characteristic  scat¬ 
tering  properties  especially  in  the  near  infrared  spectral 
region.  Because  of  the  limiting  assumptions  adopted  in  these 
studies,  a  more  comprehensive  modeling  program  was  initiated. 

The  technical  approach  adopted  consisted  of  simulating 
the  sensitivity  of  satellite  incident  visible  and  near  IR 
radiances  to  variations  in  near  surface  ambient  meteorological 
variables  by  modeling  the  physical  mechanisms  linking  the 
meteorologically-dependent  optical  properties  of  the  lower 
atmosphere  to  the  transfer  of  incident  solar  radiation  due 
to  reflection  and  scattering  from  the  atmosphere-ocean  system. 
Mechanisms  simulated  included:  (1)  the  change  in  aerosol 
size  distribution  and  complex  index  of  refraction  due  to  var¬ 
iations  in  relative  humidity,  (2)  enhanced  aerosol  coarse 
mode  loading  due  to  wind  speed  dependent  sea  spray,  (3)  vari¬ 
ations  in  aerosol  optical  properties  due  to  the  preceeding 
effects,  (4)  relative  humidity  dependent  variations  in  non- 
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aerosol  atmospheric  transmission  due  to  molecular  absorption 
by  water  vapor,  and  (5)  radiative  transfer  including  multiple 
scattering,  scattering  anisotropy,  and  surface  reflection. 

With  the  objective  of  providing  an  economical  tool  to  aid  in 
the  analysis  of  DMSP  imagery,  the  modeling  level-of-ef fort 
devoted  to  each  mechanism's  submodel  was  kept  as  simple  as 
possible. 

As  a  result  of  this  work  (Isaacs,  1980) ,  a  model  to 
relate  simulated  satellite  radiances  (in  the  wavelength 
region  0. 4-1.1  um)  to  variations  in  meteorological  variables 
such  as  relative  humidity,  wind  speed,  aerosol  number  density/ 
visual  range,  and  aerosol  scale  height  was  developed  and  im¬ 
plemented.  This  model  was  exercised  to  provide  a  theoretical 
sensitivity  analysis  of  wavelength-dependent  and  DMSP  bandpass 
weighted  radiances  to  the  above  cited  meteorological  variables. 

2.2  Model  Results 

By  consolidating  individual  model  simulations  of  DMSP 
bandpass  weighted  radiances  for  each  of  the  two  maritime 
aerosol  models  adopted,  nomograms  of  radiance,  I (DMSP),  vs. 
relevant  meteorological  variables  were  prepared.  These  are 
illustrated  in  Figures  2-1  and  2-2  for  the  maritime  and  hy¬ 
brid  aerosol  models,  respectively.  As  illustrated  in  the 
sensitivity  analysis  results,  DMSP  radiances  are  dependent 
on  a  number  of  meteorologically  based  variables.  Theoret¬ 
ically,  a  single  DMSP  radiance  data  point  can  be  expected 
to  provide  one  piece  of  information  in  the  inverse  applica- 
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tion.  Conceptually,  if  DMSP  radiances  were  strongly  asso¬ 
ciated  with  the  variation  of  one  meteorological  variable  and 
only  weakly  dependent  on  the  others,  it  would  be  possible  to 
obtain  information  about  the  controlling  independent  variable. 
The  results  indicate,  however,  that  both  relative  humidity  and 
aerosol  scale  height  are  of  equal  importance  within  the  con¬ 
text  of  the  maritime  aerosol  model  (Figure  2-1) .  Therefore, 
either  variable  cannot  be  uniquely  determined  from  DMSP  data 
alone.  The  situation  is  similar  for  the  wind  speed  dependent 
hybrid  model  (Figure  2-2)  although  in  this  case  wind  speed 
effects  are  apparently  of  second  order  with  respect  to  rel¬ 
ative  humidity  and  scale  height  dependence. 

2.3  Model  Applications 

As  written,  the  simulation  model  was  applicable  to  the 
following  cases: 

.  regions  characterized  by  maritime  aerosol  models 
(chemical  composition,  number  size  distribution, 
and  aerosol  growth  laws) , 

.  regions  remote  from  sunglint  (i.e.,  the  surface  is 
treated  as  a  Lambert  reflector  and  no  treatment  of 
surface  reflectance  on  wind  speed  dependent  surface 
slope  geometry  is  included) , 

.  the  spectral  region  extending  from  0.4um  in  the 
visible  to  l.lym  in  the  near  infrared,  i.e.  that 
of  the  DMSP  primary  sensor. 
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MODEL  VALIDATION  STUDIES 


3.1  MAGAT  Data 

Data  from  the  Monterey  Aerosol  Generation  and  Atmospheric 
Turbulence  (MAGAT)  experiment  conducted  in  the  vicinity  of 
Monterey  Bay  in  May  1980  were  employed  in  support  of  the 
evaluation  and  validation  of  the  aerosol  optical  properties 
submodel  of  the  simulation  model.  Seven  (7)  sets  of  MAGAT 
data  were  supplied  by  the  Naval  Environmental  Prediction 
Research  Facility  consisting  of:  (1)  wavelength-dependent 
(.49,  .53,  .63,  .83,  1.03  urn)  aerosol  extinction  (km  1)  cal¬ 
culated  from  Knollenberg  spectrometer  size  spectra  measure¬ 
ments  made  during  aircraft  ladder  profile  (LP)  flights  on  1, 

4  and  7  May  1980,  (2)  coincident  temperature  and  dew  point 
temperature,  (3)  water  vapor  mixing  ratio,  and  (4)  location 
of  each  ladder  profile.  Data  received  for  analysis  are  sum¬ 
marized  in  Table  3-1.  Subsequent  to  receipt  of  the  aerosol 
optical  data,  a  variety  of  corresponding  DMSP  LF  (visible 
and  near  IR)  and  thermal  imagery  was  obtained  in  hardcopy 
(i.e.,  photographic  print)  format.  Due  to  our  inability  to 
extract  correlative  quantitative  radiance  data  from  the  photo¬ 
graphic  print  imagery,  validation  exercises  focused  on  the 
aerosol  extinction  data. 

In  order  to  manipulate  the  aerosol  optical  data  (which 

were  received  in  tabular  form) ,  data  files  were  established 
for  each  data  subset.  The  object  of  the  validation  exercise 
was:  (1)  to  use  the  aerosol  optical  properties  submodel  of 
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Table  3-1  Summary  of  MAGAT  Data 


the  simulation  algorithm  to  predict  wavelength  dependent  ex¬ 
tinction  properties  of  the  atmosphere  given  the  measured  in 
situ  meteorological  variables  (e.g.,  wind  speed,  relative 
humidity)  as  model  input,  and  (2)  to  compare  the  predicted 
aerosol  extinction  to  that  actually  measured^  To  accomplish 
these  steps,  software  was  written  to:  (1)  read  in  each  MAGAT 
data  set,  (2)  evaluate  the  required  relative  humidity  profile 
from  temperature  and  dew  point  data,  (3)  evaluate  correspond¬ 
ing  model  aerosol  extinction  profiles  at  the  required  wave¬ 
lengths  and  relative  humidities  for  tropospheric,  maritime, 
rural,  and  urban  aerosol  models,  and  (4)  provide  plots  com¬ 
paring  measured  and  modeled  aerosol  optical  extinction. 

To  accomplish  the  comparison  step,  model  number  densities 
were  normalized  to  yield  measured  extinction  at  0.53  ym. 
Normalization  of  overall  results  was  necessary  since  no  data 
were  available  to  select  appropriate  total  aerosol  number 
density  values  for  each  data  set.  Thus,  comparison  results 
provide  only  relative  wavelength  consistency  of  measured  to 
modeled  results. 

3.2  Comparison  Results 

Of  the  seven  data  sets  received,  four  were  selected  for 
in-depth  analysis.  These  are  indicated  by  an  asterisk  in 
Table  3-1  and  generally  consist  of  those  ladder  profiles  not 
contaminated  by  the  presence  of  excessive  clcuds.  For  each 

4* 

The  measured  extinction  was  actually  calculated  from  size 
distribution  measurements  using  Mie  theory  (A.  Goroch,  per¬ 
sonal  communication) . 
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data  set  results  are  presented  in  this  section  consisting  of: 
(a)  profiles  of  the  measured  extinction  data  at  each  of  the 
five  available  wavelengths  (0.49,  0.53,  0.63,  0.83  and  1.03 
um) ,  (b)  a  corresponding  profile  of  relative  humidity  calcu¬ 
lated  from  the  temperature  and  dew  point  data  supplied,  (c) 
comparison  of  model  predicted  extinction  profiles  at  selected 
wavelengths  for  each  of  the  four  available  aerosl  models 
(tropospheric,  rural,  urban,  and  maritime)  based  on  the  mea¬ 
sured  relative  humidity  profile  as  model  input+  and  normal¬ 
izing  results  to  the  data  at  0.53  pm,  and  (d)  comparison  of 
the  wavelength  dependence  of  modeled  and  measured  extinction 
at  selected  pressure  levels. 

3.2.1  Ladder  Profile  5  (LP  5) 

Extinction  profile  data  for  LP  5  on  5  May  1980  are 
illustrated  in  Figure  3-1  (note:  pressure  profile  coordin¬ 
ates  (mb)  with  surface  toward  top  of  plot) .  The  correspon¬ 
ding  relative  humidity  profile  is  shown  in  Figure  3-2.  The 
comparison  of  measured  (crosses)  to  simulated  extinction  is 
given  in  Figure  3-3  for  wavelengths  of:  (a)  0.49,  (b)  0.63, 
(c)  0.83,  and  (d)  1.03  um.  The  general  slope  of  the  vertical 
extinction  profile  is  reproduced  by  all  four  available  aero¬ 
sol  models.  At  the  longer  wavelengths  (0.83  uni  -  Figure  3-3c; 
1.03  um  -  Figure  3-3d) ,  extinction  is  overestimated  by  the 

t  -1 

Light  winds  (i.e.,  <5  ms  )  were  reported  for  the  data  sets 
used  in  the  analysis  (A.  Goroch,  personal  communication) 
and  therefore,  0  ms“l  was  used  as  coarse  model  input  when 
required. 
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Figure  3-1.  Extinction  coefficient  (knr1)  Figure  3-2.  Relative  humidity  (%>  profile 

data  for  ladder  profile  15  on  for  ladder  profile  #5  on  050180 

050180.  Wavelengths:  (a)  0.49pm, 

(b)  0.53pm,  (c)  0.63pm,  (d)  0.83pm, 
and  (e)  1.03pm. 
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maritime  model  in  the  layer  near  the  surface  (p  >  990  mb) 
while  it  is  underestimated  by  the  other  models.  Note  that 
this  layer  is  a  region  of  higher  relative  humidity.  Select¬ 
ing  a  pressure  level  within  this  region  and  examining  the 
wavelength  dependence  of  the  modeled  extinction  results  in 
Figure  3-4a.  Here  extinction  data  are  selected  at  p  =  1023 
mb  near  the  surface  where  the  relative  humidity  is  about  80%. 
At  shorter  wavelengths  (<  0.7  um) ,  the  maritime  model  does 
relatively  well  in  predicting  extinction  wavelength  depen¬ 
dence.  In  the  near  IR,  however,  extinction  apparently  de¬ 
creases  with  wavelength  faster  than  the  maritime  model  per¬ 
mits.  The  wavelength  dependence  at  these  longer  wavelengths 
is  closer  to  that  at  the  urban  and  rural  models  which  have 
much  smaller  coarse  mode  number  densities.  As  Figure  3-4b 
indicates,  it  is  difficult  to  connect  some  of  the  predicted 
extinction  curves  with  the  measured  data.  These  data  are 
at  p  =  918  mb  which  is  probably  within  a  haze  layer  perceived 
by  the  surface  observer. 

3.2.2  Ladder  Profile  14  (LP  14) 

A  less  complicated  situation  is  illustrated  by  the  data 
from  LP  14  collected  on  4  May  1980.  Figures  3-5  and  3-6  il¬ 
lustrate  the  extinction  profile  data  and  calculated  relative 
humidity  profiles,  respectively,  for  a  haze  layer  extending 
from  about  600  to  1400  ft.  in  altitude  970-1010  mb) .  The 
relative  humidity  in  this  layer  is  between  80  and  85%.  If 
model  predicted  extinction  profiles  are  compared  to  the  data 
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(a)  Pessure  =  1023  mb.  Relative  Humidity  =  78.7% 


Figure  3-4.  Comparison  of  model  predicted  extinction  coefficients 
and  data  (+)  at  (a)  1023  mb  and  (b)  918.3  mb  pressure 
levels  for:  maritime  (M),  urban  (U),  rural  (R),  and 
tropospheric  (T)  aerosol  models  (normalized  at  0.53um) 
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at  wavelengths  of  0.83  (Figure  3-7a)  and  1.03  um  (Figure 
3-7b) ,  the  urban,  tropospheric,  and  rural  models  all  under¬ 
estimate  extinction,  while  the  maritime  model  performs  rather 
well  particularly  at  the  lower  relative  humidities  (from  970- 
990  mb) .  Examining  the  comparison  from  near  the  bottom  of 
the  perceived  haze  layer  (at  p  =  1013  mb)  in  Figure  3-8,  it 
may  be  noted  that  the  wavelength  dependence  of  the  data  and 
that  of  the  maritime  model  are  more  consistent  than  for  the 
previous  data  set. 

3.2.3  Ladder  Profile  23  (LP  23) 

The  data  from  LP  23  gathered  on  7  May  1980  are  inter¬ 
esting  due  to  the  presence  of  a  cloud  between  1100  and  1400 
ft.  (approximated  970-980  mb) .  The  transition  from  clear  to 
cloudy  can  be  seen  both  in  the  extinction  profile  data  in 
Figure  3-9  and  the  relative  humidity  profile  in  Figure  3-10. 
Note  the  scale  change  in  Figure  3-9  at  990  mb  where  extinc¬ 
tion  changes  by  an  order  of  magnitude.  These  data  provide 
an  excellent  test  of  the  model’s  ability  to  treat  relative 
humidity  dependent  aerosol  growth  by  comparison  in  the  high 
relative  humidity  region  below  the  cloud  (p  >  990  mb)  and  the 
low  relative  humidity  region  above  the  cloud  (p  <  990  mb) . 
Such  a  comparison  is  illustrated  in  Figure  3-11  for  extinc¬ 
tion  data  at  1.03  um  (again  note  the  scale  change  at  990  mb). 
The  maritime  model  appears  to  reproduce  the  data  well  above 
990  mb,  while  it  overestimates  the  data  near  the  surface. 

The  other  models  all  underestimate  the  measured  extinction. 
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Figure  3-7.  Comparison  of  model  predicted  extinction  coefficient  profiles  to  data(+)  at 
wavelengths  of:  (a)  0.83pm  and  (b)  1.03pm  (normalized  at  0,53pm). 


Extinction  coefficient  (km 


Figure  3-8. 


Comparison  of  model  predicted  extinction 
coefficient  and  data(+)  at  the  1013  mb 
pressure  level  for:  maritime  (M) ,  urban  (U), 
rural  (R),  and  tropospheric  (T)  aerosol 
models  (normalized  at  0.53ym). 
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Figure  3-11.  Comparison  of  model  predicted  extinction 
coefficient  profiles  with  data(+)  at 
wavelength  of  1.03ym  for:  maritime  (M), 
urban  (U),  rural  (R),  and  tropospheric  (T) 
aerosol  models  (note  scale  change  at  990mb) . 


In  Figure  3-12  data  is  compared  to  the  modeled  extinction 
coefficients  at  one  altitude  within  the  low  relative  humidity 
region  (at  p  =  916  mb,  RH  =  33%,  Figure  3-12a)  and  one  within 
the  high  relative  humidity  region  (at  p  =  1019  mb,  RH  =  84%, 
Figure  3-12b) .  The  comparison  at  high  relative  humidity  is 
like  that  observed  in  previous  cases.  Within  the  dry  layer 
aloft,  the  maritime  model  and  data  are  quite  consistent. 

3.2.4  Ladder  Profile  24  (LP  24) 

A  final  example  which  also  illustrates  the  effect  of 
cloud  is  the  LP  24  from  7  May  1980.  As  the  extinction  co¬ 
efficient  data  indicate  (Figure  3-13) ,  there  is  a  cloud 
between  about  1200  and  1500  ft.  (960-990  mb).  In  this 
region,  cloud  extinction  coefficients  are  uniformly  greater 
than  1.0  km  1  and,  hence,  are  off  scale  in  Figure  3-13. 
Corresponding  to  this  stratus  cloud  is  a  region  of  enhanced 
relative  humidity  (Figure  3-14).  Near  the  base  of  the  cloud, 
the  relative  humidity  approaches  Saturation.  This  can  be 
seen  in  Figure  3-15  where  the  relative  humidity  profile  is 
reproduced  for  a  pressure  scale  expanded  in  the  region  below 
the  cloud  deck  (i.e.,  p  >  980  mb).  The  modeled  extinction 
at  1.03pm  wavelength  significantly  underestimates  the  data 
in  this  region,  except  for  the  maritime  model  which  again 
characteristically  overestimates  (Figure  3-16).  In  the 
drier  region  above  the  cloud  (p  <  960  mb) ,  the  maritime 
model  fits  the  data  well  (Figure  3-17).  This  behavior  is 
further  supported  by  a  comparison  of  wavelength  fitting 
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Figure  3-12.  Comparison  of  model  predicted  extinction  coefficient 
and  data(+)  at:  (a)  916mb  and  (b)  1019. 2mb  pressure 
levels  for:  maritime  (M),  urban  (U),  rural  (R)  and 
tropospheric  (T)  aerosol  models  (normalized  at  0.53um). 
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Figure  3-13.  Extinction  coefficient  (km  )  data  Figure  3-14.  Relative  humidity  (%)  profile 
for  ladder  profile  #24  on  050780.  for  ladder  profile  #24  on 

Wavelengths:  (a)  0.49pm,  (b)  0.53pm,  050780. 

(c)  0.63pm,  (d)  0.83pm,  and  (e)  1.03pm 
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Figure  3-17.  Comparison  of  model  predicted  extinction 
coefficient  profiles  to  data (+)  at  wave¬ 
length  of  1.03ym.  (Low  relative  humidity 
region,  RH<45%,  normalized  at  0.53  m.) 
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(b)  Pressure  =  1002mb,  Relative  Humidity  =  91.3% 

Figure  3-18.  Comparison  of  model  predicted  extinction  coefficient 
and  data(+)  at:  (a)  932. lmb  and  (b)  1002mb  pressure 
levels  for:  maritime  (M) ,  urban  (U),  rural  (R),  and 
tropospheric  (T)  aerosol  models  (normalized  at  0.53um). 
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in  the  low  RH  region  above  the  cloud  at  p  =  932  mb  and  RH  = 

32%  (Figure  3-18a)  to  that  in  the  high  RH  region  below  the 
cloud  at  p  =  1002  mb  and  RH  =  91%  (Figure  3-18b) . 

3.3  Discussion 

Prediction  of  wavelength  dependent  aerosol  optical  prop¬ 
erties  (including  extinction)  based  on  meteorological  input 
data  is  the  kernel  of  existing  optical  propagation  and  trans¬ 
mission  algorithms.  Furthermore,  as  illustrated  in  Figure  2-1 
(and  of  particular  relevance  to  the  current  study) ,  it  is  a 
key  step  in  the  approach  to  simulating  remotely  sensed  radi¬ 
ances  from  meteorological  input  parameters.  For  this  reason, 
the  comparisons  of  measured  and  model  predicted  aerosol  ex¬ 
tinction  discussed  above  provide  valuable  insight  into  the 
simulation  problem. 

Based  on  the  relatively  small  sample  of  comparisons  pre¬ 
sented,  it  is  gratifying  to  see  that  the  maritime  aerosol 
model  performs  relatively  well  where  expected.  It  is  clear 
that  for  the  MAGAT  data  sets  examined  predictions  based  on 
the  rural,  tropospheric,  and  urban  models  are  not  consistent 
with  the  measured  extinction  in  wavelength  dependence.  This 
is  almost  certainly  due  to  their  lack  of  appropriate  oceanic 
origin  coarse  mode  particles. 

It  is  noteworthy,  however,  that  the  maritime  model  also 
exhibits  deficiencies  which  are  manifested  in  overestimation 
of  extinction  at  longer  wavelengths  (as  contrasted  with  under¬ 
estimation  by  the  other  models)  particularly  at  high  relative 


humidity  (e.g.,  greater  than  80%).  This  may  be  attributable 
to  either  too  many  particles  in  the  oceanic  origin  coarse 
mode  or  to  too  large  a  relative  humidity  dependent  aerosol 
growth  factor  (see  Figure  4-2b)  assumed  in  the  model. 


4.  CONTINENTAL  AEROSOL  MODELS 


4.1  Non-Maritime  Origin  Aerosols 

While  investigating  the  effect  of  low  level  maritime 
haze  in  DMSP  visible  imagery,  two  oceanic  origin  aerosol 
size  models  were  treated:  (a)  the  relative  humidity  depen¬ 
dent  maritime  model  of  Shettle  and  Fann  (1979) ,  and  (b)  a 
hybrid  model  consisting  of  a  fine  mode  identical  to  that  of 
Shettle  and  Fenn' s  (1979)  tropospheric  model  and  a  wind  speed 
dependent,  sea  spray  origin  coarse  mode  adopted  from  the  Munn- 
Katz  model  (cf.  Wells  et  al.,  1977).  The  physical,  chemical, 
and  optical  properties  of  these  aerosol  models  are  described 
in  Sections  4  and  5  of  Isaacs  (1980)  and  are  implemented  with¬ 
in  the  computer  code  described  in  Section  7  of  this  report. 

It  is  notable  that  regions  over  the  open  ocean  are  often 
under  the  influence  of  air  masses  of  continental  origin.  This 
is  particularly  true  of  coastal  regions.  A  variety  of  examples 
is  presented  in  Fett  and  Mitchell  (1979)  illustrating  the  ef¬ 
fect  of  continental  origin  dust,  smoke,  and  other  atmospheric 
pollutants  on  DMSP  VHR  imagery.  These  include,  for  example, 

scirocco  conditions  in  North  Africa  and  Santa  Ana  outflows 
over  the  Los  Angeles  bight.  Other  fugitive  dust  conditions 
over  the  open  oceans  have  been  documented  (cf . ,  Fraser,  1976; 
Duce  et  al. ,  1980) .  Recognizing  that  analysis  of  imagery  in 
coastal  regions  is  often  required  in  support  of  tactical 
naval  operations,  it  seems  appropriate  to  include  the  capa¬ 
bility  to  treat  such  continental  origin  aerosols  within  the 


<# 


simulation  model.  For  this  purpose,  the  domain  of  applicable 
situations  was  extended  by  introducing  relative  humidity  de¬ 
pendent  aerosol  optical  data  into  the  simulation  algorithm 
for  both  rural  and  urban  aerosol  models.  The  rural  and  urban 
aerosol  models  differ  fundamentally  from  the  maritime  model 
both  in  chemical  composition  and  size  distribution:  for  ex¬ 
ample,  the  composition  dependent  complex  index  of  refraction 
of  each  model  is  quite  distinct  (see  Figure  4-1) .  In  the 
following  sections,  the  specific  optical  data  implemented  for 
each  aerosol  model  and  relevant  simulation  results  based  on 
them  are  presented.  Details  of  the  physical  and  chemical 
basis  of  these  models  are  described  in  Shettle  and  Fenn  (1979) . 
A  review  of  the  theory  relating  physical  aerosol  models  and 
their  respective  optical  properties  is  given  in  Isaacs  (1980) . 

4.2  Rural  Aerosol  Model 

The  rural  aerosol  model  adopted  from  Shettle  and  Fenn 
(1979)  is  described  es  applicable  to  situations  not  directly 
influenced  by  urban  and/or  industrial  aerosol  sources.  Its 
chemical  composition  includes  both  water  soluble  salts  and 
dust-like  components  in  a  bimodal,  log  normal  size  distribu¬ 
tion.  The  size  distribution  is  illustrated  in  Figure  4-2 
with  relative  humidity  dependent  mode  radii  in  fine  and 
coarse  modes  given  by  the  growth  factors  illustrated  in 
Figures  4-3a  and  4-3b,  respectively. 

The  relevant  wavelength-dependent  aerosol  optical  prop¬ 
erties  based  on  the  rural  model  and  normalized  to  a  total 
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number  density  of  1.5  x  10  cm  which  are  incorporated  with¬ 
in  the  simulation  model  are  the:  (a)  extinction  coefficient 


d 
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(a)  Real  part  of  the  complex  index  of  refrac¬ 
tion  vs.  relative  humidity. 


(b)  Imaginary  part  of  the  complex  index  of 
refraction  vs.  relative  humidity. 


Figure  4-1.  Dependence  of  the  complex  index  of  refraction 
of  maritime  (-) ,  urban  (o) ,  and  rural  (+) 
aerosol  models  on  relative  humidity. 
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RURAL  AEROSOL  MOOELS 


Figure  4-2.  Aerosol  number  distribution 
(cm”  3  ym“l)  for  rural  model  as  a  functic 
of  relative  humidity  with  N  =  1. 5  x 104  c 
(from  Shettle  and  Fenn,  1979) . 


(a)  Fine  mode  growth  factors  for 

rural  (+)  and  urban  (o)  models. 
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(b)  Coarse  mode  growth  factors 
for  rural  (-) ,  urban  (o) , 
and  maritime  (+)  models. 

Figure  4-3.  Aerosol  model  relative  humidity 
growth  factors  (Shettle  and  Fenn,  1979) . 
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and  (o'  phase  function  asymmetry  parameter  (Figure  4-4c) . 

These  are  given  as  functions  of  relative  humidity  for  discrete 
values  of  0,  50,  70,  80,  90,  95,  98  and  99%.  Employing  the 
simulation  model,  sensitivity  analyses  were  conducted  in 
analogy  to  those  performed  for  the  oceanic  aerosol  models 
in  Isaacs  (1980) . 

Figures  4- 5a  through  4-5d  illustrate  the  wavelength- 
dependent  radiance  spectra  for  relative  humidities  between 
50  and  95%  assuming  aerosol  scale  heights  of  (a)  0.5,  (b)  1.0, 

(c)  1.5,  and  (d)  2.0  km,  respectively.  The  sun/ sensor  geom¬ 
etry  corresponds  to  a  nadir  viewing  sensor  (y  =  1.0)  and  a 
solar  zenith  angle  cosine  of  y0  =  0*5  as  in  the  previous 
sensitivity  studies.  The  results  are  qualitatively  similar 
to  those  obtained  for  the  maritime  aerosol  model  although  on 
closer  inspection  it  is  noted  (see  Figures  4-6a  and  b)  for  a 
given  scale  height,  sensor  incident  radiances  for  the  rural 
model  (particularly  in  the  near  IR  spectral  region)  are  much 
less  sensitive  to  relative  humidity.  This  is  primarily  due 
to  the  larger  coarse  mode  relative  humidity  growth  factor 
of  the  maritime  model  compared  to  that  of  the  rural  model 
(see  Figure  4-3b)  . 

4.3  Urban  Aerosol  Model 

The  urban  aerosol  model  considers  the  modification  of 
aerosol  optical  properties  when  combustion  products  from 
various  anthropogenic  and  perhaps  natural  sources  introduce 
soot-like  aerosol  components. 
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Figure  4- 4a. 

Extinction  coefficient 
(krrrl)  -  Rural  model 
(from  Shettle  and 
Fenn,  1979) 
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Figure  4-4b. 
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Figure  4-4c. 

Asymmetry  Parameter 
Rural  model  (from 
Shettle  and  Fenn, 
1979) 
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Figure  4- 5a.  Dependence  on  Relative  Humidity 
(Rural  Model,  H  =  0.5) 
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Figure  4- 5b.  Dependence  on  Relative  Humidity 
(Rural  Model ,  H  =  1.0) 
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Figure  4-6 .  Comparison  of  radiance  spectra  for 
(a)  maritime,  (b)  rural,  and  (c)  urban  aerosol 
models  as  functions  of  relative  humidity  (50  to 
95%)  for  H  =  1.5  km,  y  =  1.0,  and  y0  =  0.5. 
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Its  size  distribution  is  assumed  identical  to  that  of  the 

4 

rural  model,  however,  the  total  number  density  is  2.0x10 
cm  ^  (one  third  higher  than  that  of  the  rural  model  and  five 
times  higher  than  the  maritime  model) .  But  for  the  number 
density  difference,  the  distinguishing  characteristic  of  the 
urban  model  is  chemical,  i.e.,  the  presence  of  carbon.  Not¬ 
ably,  recent  studies  suggest  that  soot-like  aerosol  components 
are  quite  pervasive  in  the  atmosphere  (cf.  Rosen  et  al.,  1981) 
and  thus,  the  optical  properties  of  the  urban  aerosol  model 
may  be  appropriate  over  regions  hitherto  assumed  remote  from 
urban  influence. 

As  a  consequence  of  its  soot  component,  two  important 
effects  are  noted.  As  indicated  in  Figure  4-3,  the  relative 
humidity  growth  factor  for  both  the  fine  (Figure  4-3  a)  and 
coarse  (Figure  4-3 b)  modes  are  distinct  from  the  other  models 
considered.  This  is  due  to  the  dependence  of  aerosol  growth 
on  chemical  composition.  Secondly,  carbon  is  quite  absorbing 
(Bergstrom,  1973) .  This  can  be  seen  by  examining  the  extinc¬ 
tion  coefficient  and  single  scattering  albedo  at  low  relative 
humidities  in  Figures  4-7a  and  4-7b,  respectively.  (The  im¬ 
plemented  optical  properties  of  the  urban  model  are  given  in 
Figures  4-7a,b,c,  respectively.)  Such  high  absorption  sug¬ 
gests  a  potential  adverse  impact  on  optical  transmission. 

The  sensitivity  analyses  performed  on  scale  height  and 
relative  humidity  for  a  nadir  viewing  sensor  (y  =  1.0)  and 
solar  zenith  angle  cosine  of  (u  =  0.5)  are  illustrated  in 
Figure  4-3  for  scale  heights  of:  (a)  0.5,  (b)  1.0,  (c)  1.5, 
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Figure  4-7 a. 

•  Extinction  coefficient 
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Figure  4 -7  b. 

&  3 


Single  scattering 

albedo  -  Urban  model 

.  3  5 

(from  Shettle  and 

Fenn,  1979) 

o 

t3  ~  (T 

aj  i 

-Q 

5  .  7 

.  6  5 


95 


4  .5  .6  .7  .8  9  1  1.1 


Wavelength  (um) 


Figure  4-7c. 

Asymmetry  parameter 
Urban  model  (from 
Shettle  and  Fenn, 

,,  1979) 


t  < 


S- 

0) 


OJ 


E 


<o 

s~ 

fO 

CL 


2? 

+-> 

<U 


CO 

< 


74 
7  2 

i*' 

.  6  8 

.  *  6 

64 
.  62 

.  6 


« 


4-12 


ance 


Wavelength  (urn) 

Figure  4- 8a.  Dependence  on  Relative  Humidity 

(Urban  Model,  H  =  0.5) 
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Figure  4-8b.  Dependence  on  Relative  Humidity 

(Urban  Model ,  H  =  1.0) 
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Figure  4-8d 


Dependence  on  Relative  Humidity 
(Urban  Model ,  H  =  2.0) 


and  (d)  2.0  km,  respectively.  Results  show  wavelength  depen- 
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dent  radiance  spectra  for  relative  humidities  of  50,  70,  80, 

90  and  95%.  Qualitatively,  they  are  not  remarkably  different 
from  those  for  the  maritime  and  rural  models.  Inspection  of 
Figure  4-6,  however,  suggests  a  few  notable  distinctions. 

The  relative  humidity  dependence  of  the  urban  aerosol  fine 
and  coarse  model  growth  factors  is  somewhat  higher  than  that 
of  the  rural  model  although  it  is  not  as  large  as  that  of 
the  maritime  model  (see  Figures  4-3a,b) .  Nevertheless,  in 
the  resultant  radiance  calculation,  it  is  apparent  that  radi¬ 
ance  spectra  sensitivity  to  increases  in  relative  humidity 
given  the  urban  aerosol  (Figure  4-6c)  are  much  closer  to  those 
of  the  maritime  aerosol  (Figure  4-6a)  than  are  those  for  the 
rural  aerosol  model  (Figure  4-6b) .  This  is  true  both  in  the 
visible  and  near  IR  regions.  As  a  consequence  of  the  presence 
of  absorbing  carbon,  however,  radiance  spectrum  levels  in  the 
visible  at  low  relative  humidities  (i.e.,  50%)  are  appreciably 
lower  when  backscattered  from  the  urban  atmosphere.  At  higher 
relative  humidities,  single  scattering  albedos  increase  due  to 
the  additional  water  volume  of  the  aerosol  and  the  effect  is 
not  as  drastic.  It  can  be  seen,  however,  in  Figure  4-9  which 
directly  compares  radiance  spectra  for  the  four  adopted  aero¬ 
sol  models. 

4.4  Aerosol  Model  Intercomparisons. 

The  distinct  wavelength  and  relative  humidity  dependent 
optical  properties  of  the  hybrid,  maritime,  rural  and  urban 
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Figure  4-9.  Simulated  Radiances  for:  (1)  hybrid, 

(2)  maritime,  (3)  rural,  and  (4)  urban, 
aerosol  models.  [ RH  =  50%,  WS  =  5  ms  . 
HI  =  1.2,  H2  =  0.8,  y  =  1.0,  pQ  =  0.5, 

A<t>  =  0] 
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aerosol  models  are  manifested  in  their  respective  radiance 
spectra  signatures  when  viewed  over  the  ocean  surface.  In 
the  previous  comparisons,  the  degree  of  sensitivity  to  in¬ 
creases  in  relative  humidity  has  been  noted.  This  section 
compares  the  wavelength  dependence  of  the  various  aerosol 
models  discussed  above.  Radiance  spectra  are  presented  for 
a  relative  humidity  of  50%  and  fixed  scale  height  of  1.2  km 
for  the  fine  mode  of  the  hybrid  model  and  both  modes  of  other 
models.  For  the  hybrid  model  a  scale  height  of  0.8  km  is 
assumed  for  the  coarse  mode  and  a  wind  speed  of  5  ms  ^  is 
used  to  evaluate  the  loading  in  this  size  range.  Table  4-1 
presents  the  number  densities  and  corresponding  visual  ranges 
represented  by  each  curve  in  Figure  4-9. 

A  number  of  observations  can  be  made  based  on  this 
aerosol  model  intercomparison.  Given  the  variey  of  aerosol 
chemical  and  size  distribution  models  and  the  range  of  number 
densities  and  visual  ranges  indicated,  the  radiance  spectra 
in  Figure  4-9  are  remarkably  similar.  Although  it  might  be 
expected  that  at  a  given  wavelength  the  radiance  would  be  pro¬ 
portional  to  loading  (i.e.,  total  number  density),  this  is 
not  so  for  a  variety  of  reasons.  In  the  visible,  for  example, 
although  the  urban  model  has  the  largest  particle  density,  its 
radiance  is  lowest  of  the  four  models  due  to  absorption.  The 
highest  radiance  is  from  backscattering  from  the  model  with 
the  next  largest  loading,  the  rural  model.  In  the  near  IR, 
the  number  of  larger  particles  (i.e.,  in  the  coarse  mode) 
determines  to  a  great  extent  the  backscattered  radiance.  As 


Model 


Visual  Range 
(tan) 


Total  Number  Density 
(cm  “3) 


hybrid 

5.0  x  10 3 

59.4 

maritime 

4.0  x  103 

39.6 

rural 

1.5  x  104 

23.9 

urban 

2.0  x  104 

19.9 

Table  4-1.  Number  densities  (cm  3)  and  visual 
ranges  corresponding  to  aerosol  models  (1-4)  in 
Figure  4-8. 


Figure  4-9  indicates,  the  hybrid  model  is  predicted  to  have 
the  smallest  number  of  coarse  mode  particles  which  is,  in 
fact,  the  case.  As  the  discussion  suggests,  much  information 
concerning  the  properties  of  the  ambient  aerosol  is  poten¬ 
tially  available  in  the  wavelength  dependent  radiance  spectra 
This  additional  information  content  is  essentially  destroyed 
when  integrated  by  a  broad  band  sensor  such  as  DMSP. 
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5.  SUN/SENSOR  GEOMETRY  FOR  DATA 
SIMULATION  APPLICATIONS 
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Previous  sensitivity  analyses  results  based  on  the  simu¬ 
lation  model  were  evaluated  for  a  fixed  generic  sun/sensor 
geometry  (cf .  Isaacs,  1980) .  For  the  purpose  of  providing 
an  applications  oriented  analysis  tool,  however,  the  present 
study  focuses  on  program  modifications  which  directly  model 
the  sun/sensor  geometry  characteristic  of  specific  polar- 
orbiting  satellite  platforms  (such  as  DMSP) .  Thus  given 
basic  information  regarding  the  orbital  parameters  of  the 
satellite  and  the  desired  time  (since  the  imagery  is  annotated 
in  GMT  time) ,  a  method  is  sought  to  both  locate  (i.e.  lati¬ 
tude/longitude)  specific  surface  elements  (pixels)  along  the 
sensor  scan  line  perpendicular  to  the  satellite  subpoint  track 
and  to  provide  information  regarding  the  solar  and  sensor 
positions  relative  to  each  of  the  individual  pixels.  These 
relative  solar  and  sensor  positions  are  then  used  as  the 
basis  of  the  radiative  transfer  calculations  to  follow. 

The  problem  is  simplified  considerably  if  the  satellite 
subtrack  is  assumed  to  be  a  great  circle.  In  this  case  con¬ 
cepts  from  spherical  trigonometry  may  be  applied  to  both  sub¬ 
track  and  scan  line  location  calculations.  As  discussed  in 
Tsui  and  Fett  (1980) ,  who  used  the  same  approach  in  locating 
the  primary  specular  point  (PSPX,  the  magnitude  of  the  errors 
introduced  by  adopting  this  assumption  should  be  acceptable 
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for  practical  purposes,  especially  in  light  of  other  sources 
of  error.  The  calculation  may  be  divided  into  three  funda¬ 
mental  steps:  (1)  location  of  the  satellite  subtrack  posi¬ 
tion  for  a  given  time,  (2)  specification  of  the  sensor  scan 
line  pixels  (with  arbitrary  density  along  the  scan  line) 
corresponding  to  that  time,  and  (3)  calculation  (for  each 
pixel)  of  the  spherical  coordinates  (relative  to  a  fixed 
coordinate  system)  determining  the  local  sun  and  sensor 
positions.  These  calculations  are  discussed  in  the  following 
sections  with  reference  to  Figure  5-1  which  illustrates  the 
satellite  subtrack  (AB)  and  sensor  scan  line  (BE)  geometry 
for  a  generic  polar  orbiting  satellite  with  ascending  node 
A.  Indicated  are  the  positions  of  the  pole,  C(P),  equator 
(E) ,  and  prime  meridian  (M)  . 

5.1  Satellite  Subtrack  Position 

Assuming  the  satellite  subtrack  (SST)  in  Figure  5-1 
(arc  AB)  to  be  a  great  circle,  the  position  in  latitude 
(9  )  and  longitude  ($  )  of  point  B  at  time  t.  can  be  calcu- 
lated  given  the  time,  t  ,  and  longitude,  <j>Q,  of  the  ascending 
node  (crossing  of  the  equator,  point  A)  provided  an  orbital 
period,  p,  and  inclination  of  the  orbit  to  the  equatorial 
plane,  e,  are  assumed.  For  DMSP ,  a  nominal  circular  orbit 
(although  rarely  achieved  in  practice)  would  be  characterized 
by  a  period  of  101.35m  (or  1.6982h),  an  orbital  inclination 
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of  98.7°  (making  it  sun  synchronous),  and  an  altitude  of 
833km. 

Referring  to  the  figure,  the  length  of  side  C  (in 
degrees)  will  be: 

C  =  360 (tj-t  )/p  (5-1) 

The  angle  at  A  is  determined  by  the  orbital  inclination  or 

A  =  e  -  90  =  8.7°  (5-2) 

The  arc  length  (in  degrees)  from  point  A  (on  the  equator) 
to  the  pole  C(p)  is  side  b,  clearly: 

b  =  90°  (5-3) 

The  1<  titude  of  the  subsatellite  point  at  t-^  will  be  given 
by  side  a: 


9  =  90-a  (5-4) 

s 

while  the  longitude  will  depend  on  the  angle  at  point  C(p) 
or : 


+  c  +  ‘W15 


(5-5) 


(The  final  term  above  accounts  for  the  earth's  rotation  at  a 
rate  of  15°h  ^  during  the  time  interval  from  tQ  to  t^) .  Solu¬ 
tion  for  the  unknown  quantities  (a  and  C,  above)  is  based  on 
Napier's  analogies  for  a  spherical  triangle  (Selby,  196  7; 
p.  223)  : 
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a  =  2  tan"1  [  h(b,c,B,C)  ] 

C  =  tan  1  [f(b,c,A)J  -  tan"1  (g(b,c,A)] 

where : 

B  =  tan  1  lf(b,c,A)]  +  tan  1  [g(b,c,A)J 
and  the  functions  f,  g,  and  h  are  given  by: 


(5-6) 

(5-7) 


(5-8) 


f (b,c,A)  =  tan  j  (B  +  C)  =  cos  ~  (b-c)  sec  j  (b+c)  cot  jA 

(5-9) 

g(b,c,A)  =  tan  j  (B-C)  =  sin  ~  (b-c)  esc  j  (b+c)  cot  j  A 

(5-10) 

h(b,c/B,C)  =  tan  j  a  =  tan  j  (b-c)  sin  j  (B+C)  esc  j  (B-C) 

(5-11) 

Upon  substitution  of  the  known  quantities  c  (5-1) ,  A  (5-2) , 
and  b  (5-3)  into  (5-8)  through  (5-11) ,  the  unknowns  a  (5-6) 
and  C  (5-7)  can  be  evaluated.  Equations  (5-4)  and  (5-5) 
then  yield  the  desired  satellite  subtrack  point  latitude  and 
longitude,  respectively.  Sample  satellite  subpoint  track 
latitude/longitude  calculated  as  a  function  of  time  using 
the  procedure  described  above  are  presented  in  Table  5-1. 


5.2  Sensor  Scan  Line  Position 

Applying  similar  considerations,  the  location  of  speci¬ 
fic  surface  elements  along  a  sensor  scan  line  can  be  evaluated. 
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Table  5-1 


SATELLITE  SUBPOINT  TRACK  LATITUDE  AND  LONGITUDE 
FOR  DMSP  SATELLITE  IMAGE  FOR  22  AUGUST  1978 


Input  Data: 


<fr0  =  113. 5W 

t  =  18.92  GMT 
o 


Time 

Latitude/Longitude 

Time 

Latitude/Longitude 

185827 

10.3/115.8 

190427 

31.3/121 .1 

185927 

13.8/116.6 

190527 

34.8/122.1 

190027 

17.3/117.5 

190627 

38.3/123.2 

190127 

20.8/118.3 

190727 

41.8/124.3 

190227 

24.3/119.2 

190827 

45.3/125.6 

190327 

27.8/120.1 

190927 

48.7/127.0 

Table  5-2 

SCANLINE  PIXEL  LOCATIONS  CORRESPONDING  TO 
SUBTRACK  POINT  AT  185827  GMT 


Pixel  Latitude/Longitude  Pixel _ Latitude/Longitude 


10.3/115.8 
10.4/115.2 
10.5/114.5 
10.6/113.8 
10.7/113.2 
10.8/112.5 
10.9/111.8 
1 1 .0/1 11.1 
11.1/110.5 
11.1/109.8 


11.2/109.1 
11 .3/108.5 
1 1.4/107.8 
11 .5/107.1 
1 1 .6/106.4 
11 .6/105.8 
11.7/105.1 
1 1 .8/104.4 
11 .9/103.7 
1 1 .9/103.1 
12.0/102.4 


* 


It  is  assumed  that  the  location  and  time  (8  ,  <j>  ,  t  )  of  the 

s  s  s 

subtrack  point  (B  in  Figure  5-1)  associated  with  the  desired 
sensor  scan  line  are  known  based  on  the  procedure  described 
in  the  previous  section  and  that  the  scan  line  including 
point  B  (i.e.  arc  BE)  is  also  a  great  circle  perpendicular 
to  the  subpoint  track.  Considering  spherical  triangle 
BEC(P) ,  the  length  of  side  A  (deg)  from  (5-4)  is: 

a  «  90-0  (5-12) 

s 

The  angle  at  B  is  determined  by  the  orbital  inclination  and 
thus  from  the  figure: 

B  =  81.3°  (5-13) 

The  length  of  the  scan  line  (in  degrees) ,  side  f ,  is 
determined  from  the  length  of  the  scan  from  subpoint  to 
horizon  (1479km)  and  the  earth's  radius,  R^,  (6370km): 

f  =  360  (1479) /2tt  (6370)  (5-14) 

=  13.303° 

The  ith  surface  element  pixel  is  chosen  such  that  the  length 
of  the  arc  from  the  subpoint  B  to  the  pixel  is  given  by: 

fL  (deg)  =  13.303  (5-15) 

for  i  =  1  to  n  +  1.  Thus,  the  first  pixel  (i  =  1)  corres¬ 
ponds  to  the  subtrack  point  (B)  and  the  last  pixel  (n  +  1) 
corresponds  to  the  point  on  the  easterly  horizon  (E) .  The 
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latitude  of  the  ith  pixel  along  the  scan  line  will  be: 


0i(deg)  =  90-di  (5-16) 

and  the  corresponding  longitude  will  be: 

<j> .  (deg)  =  <j>  -C.  (5-17) 

1  S  X 

Again  using  Napier's  analogies  to  solve: 

di  »  2tan“1[h(a,fi,E,Ci)]  (5-18) 

C±  =  tan"1 [f (a,fi/B) ]  -  tan"1 [g (a, f± ,B)]  (5-19) 

where : 

El  =  tan"1^ (a^^B)  ]  +  tan"1  [g  (a,  f±  ,B)]  (5-20) 


and  the  functions  f,  g,  and  h  are  given  by: 
f  (a,fi,B)  =  tan  ~  (Ei  + 

=  cos  j  (a-f^)  sec  j  (a+f^)  cot  j  B  (5-21) 
gUff^B)  =  tan  ~ 

=  sin  j  (a-f^)  esc  j  (a+f^)  cot  j  B  (5-22) 

h (a, f . ,E . ,C. )  =  tan  \  d 

=  tan  j  (a-fi)  sinj(Ei+Ci)  esc  j  (E^C.p 

(5-23) 

Substitution  of  (5-18)  through  (5-23)  into  (5-16)  and 
(5-17)  yield  pixel  latitude  and  longitude.  Calculations  for 
the  scan  lire  beginning  at  the  subtrack  point  corresponding 
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to  185827  GMT  in  Table  5-1  and  considering  21  equally  spaced 
pixels  are  given  in  Table  5-2. 


5 . 3  Pixel-related  Sun/Sensor  Geometry 

Evaluation  of  pixel  location  (0^,4^)  using  the  procedures 
outlined  in  the  previous  sections  provides  the  necessary  data 
to  calculate  the  solar  zenith  angle,  0^,  sensor  zenith  angle, 
0^,  and  azimuth  angle  difference,  A4- ,  relative  to  each  pixel 
location.  These  sun/sensor  geometric  parameters  are  required 
as  input  to  the  radiative  transfer  calculation  performed  in 
the  simulation  and  additionally  specify  the  sun  and  sensor 
position  vectors  (s,  o,  respectively;  see  Figure  7-1,  Section 
7.3)  used  in  the  surface  reflection  model. 

The  solar  zenith  angle  at  the  ith  pixel,  0*,  and  its 
associated  zenith  angle  cosine,  y^  =  cos0^  are  given  by: 


01  =  cos*"^[y1]  =  90-sin  1  [sin0  .  sin<5  +  cos0  •  cosScosn-  ]  (5-24) 

O  *  O  j,  1  x 


.•  -lr_; 


where  6  is  the  solar  declination  angle  and  is  the  hour 
angle  (deg)  given  approximately  by  the  difference  between 
the  pixel  longitude  (4>^)  and  the  longitude  (4>SSp)  of  the 
subsolar  point  (SSP) ,  i.e.:"*" 


ni  *  *i  “  ^SSP 


(5-25) 


+For  morning  satellites,  nj.will  generally  be  positive?  that 
is  the  scan  line  pixels  will  be  at  earlier  solar  times  than 
the  SSP. 
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The  solar  declination  angle,  6,  is  a  function  of  the  date 
(zero  at  the  equinoxes  and  ±23°27'  at  the  solstices)  and 
may  be  obtained  from  a  solar  ephemeris.  It  is  given  approx¬ 
imately  by  the  equation: 


5(d)  =  sin-1{ .3978  sin  [. 9863  (d-80) j }  (5-26) 

where  d  is  the  Julian  day  of  the  year  (i.e.  January  1  is  001 
and  December  31  is  365) . 

The  solar  azimuth  angle,  cu  (the  azimuth  of  the  sun 
measured  clockwise  from  north,  see  Figure  7-1) ,  at  pixel 
location  (9^,^)  is  given  by: 


ai  =  cos-1 (1— 2q5)  for  ni > 0  (morning) 

(5-27) 

360-cos_1(l-2q5)  < 0  (afternoon) 

where : 

q5  =  sin(q4-q2)  sin^-q^ 
sinq2  sinq^ 

q4  =  (ql  +  q2  +  q3)/2 
and:  q^  =  90-9^ 

<*2  *  6o 
q3  -  90-6 


Since  the  sensor  azimuth  looking  eastward  is  fixed  at  261.3°, 
the  azimuth  angle  difference  is: 


The  satellite  zenith  angle  at  the  ith  pixel,  0^,  and 
its  associated  zenith  angle  cosine,  y1  =  cos0^  are  given  by: 

0j  =  cos_1y1  =  tan"1!  1479/833]  (5-29) 

k  1  n 

for  i  =  1 ,2 , 3 . . . (n+1)  pixels  and  assuming  a  sensor  altitude 
of  833km  for  a  nominal  circular  orbit. 


t 
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6.  RADIATIVE  TRANSFER 


6.1  Technique  Selection  Criteria 

Radiative  transfer  theory  provides  a  mathematical  de¬ 
scription  of  the  interaction  between  incident  solar  radiation 
and  the  relevant  optically  active  constituents  of  the  atmo¬ 
sphere  including  both  ambient  gases  and  aerosols.  The  inter¬ 
action  mechanisms  treated  by  radiative  transfer  theory  have 
been  alluded  to  in  the  preceeding  discussions  of  atmospheric 
optical  properties  and  include  molecular  (Rayleigh)  scatter¬ 
ing,  gaseous  absorption,  and  aerosol  absorption  and  scatter¬ 
ing.  In  the  context  of  the  present  radiance  simulation  model, 
the  radiative  transfer  sub-package  facilitates  a  quantitative 
link  between  the  meteorologically  dependent  optical  properties 
of  the  atmosphere  (especially  aerosol  content)  and  sensor 
incident  backscattered  solar  radiation.  For  this  purpose, 
it  is  minimally  required  that  the  radiative  transfer  treat¬ 
ment  which  is  implemented  is  capable  of  handling:  (a)  mul¬ 
tiple  scattering,  (b)  the  inherent  anisotropic  (i.e.  highly 
directional)  scattering  characteristic  of  aerosols,  (c)  re¬ 
flection  of  radiation  at  the  atmosphere-ocean  interface,  and 
(d)  the  azimuthal  dependence  given  by  the  sun/sensor/geometry. 

The  general  attributes  cited  above  require  a  solution 
to  the  radiative  transfer  equation  ( RTE)  for  a  multiple 
scattering  atmosphere  (cf.  Chandrasekhar,  1950;  Goody,  1964; 
Liou,  1980) .  An  extensive  hierarchy  of  potential  computa¬ 
tional  methods  exists  to  treat  the  RTE  ranging  from  simple. 


6-1 


single  scattering  analyses  (such  as  that  in  Fett  and  Isaacs, 
1979)  to  a  variety  of  highly  numerical  computational  algo¬ 
rithms  (Hansen  and  Travis,  1974;  Lenoble,  1977).  Numerical 
solutions  generally  provide  high  accuracies  (less  than  1% 
difference  between  competing  state-of-the-art  methods  for 
identical  input  data)  at  the  expense  of  considerable  compu¬ 
ter  resources.  Highly  accurate  modeled  radiances  are  useful, 
of  course,  only  when  the  accuracies  of  input  atmospheric 
optical  data  are  consistently  high.  In  general,  the  tradeoff 
analysis  is  between  desired  accuracy  and  computational  effi¬ 
ciency  (Filmore  and  Karp,  1980) . 

A  variety  of  analytical  solutions  to  the  RTE  are  avail¬ 
able  for  simulation  model  implementation  which  require  very 
little  computational  effort.  However,  their  use  does  incur 
a  penalty  in  the  form  of  associated  errors  and  uncertainties 
(Isaacs  and  Ozkaynak,  1980)  .  Often,  application  of  these 
methods  is  precluded  by  their  low  levels  of  achievable 
accuracy. 

Implementation  of  radiative  transfer  theory  within  lar¬ 
ger  operational  models  introduces  severe  constraints  on  the 
available  computational  level  of  effort.  For  this  reason 
investigators  have  increasingly  relied  on  approximate,  ana¬ 
lytical  solutions  to  the  RTE  to  provide  tractable  simulation 
algorithms  (Kaufman,  1979,  1981;  Hering,  1981;  Isaacs,  1981). 

In  selecting  an  appropriate  approach  in  this  investiga¬ 
tion,  criteria  have  included:  (a)  ability  to  satisfactorily 
treat  mechanisms  cited  above,  (b)  maximization  of  the  number 


of  sensitivity  analyses  which  could  be  undertaken,  (c)  opti¬ 
mization  of  available  computer  resources,  and  (d)  consistency 
of  approach  accuracies  with  those  in  other  components  of  the 
simulation  model  (such  as  in  specifying  the  physical  aerosol 
models) . 

An  additional  consideration  in  the  context  of  the  current 
study  was  the  possible  application  of  the  radiance  simulation 
model  discussed  here  to  the  operational  analysis  of  DMSP 
imagery  for  the  purpose  of  inferring  surface  propagation 
properties.  For  this  reason  an  algorithm  implementable  on 
a  desk  top  type  computer  was  desired.  These  considerations 
led  to  selection  of  an  approximate,  analytical  solution  to 
the  RTE . 

In  the  following  sections  the  adopted  radiative  trans¬ 
fer  formulation  is  discussed  and  accuracies  are  assessed 
by  comparison  with  highly  accurate  numerical  methods.  These 
comparisons  are  performed  for  a  representative  set  of  cases 
including:  (a)  variations  in  aerosol  optical  depth  represen¬ 

tative  of  optically  thin  (x  =  0.1)  and  optically  thick 
(x  =  1.0)  situations,  (b)  isotropic,  Rayleigh,  and  aerosol 
(i.e.,  anisotropic)  scattering  (c)  conservative  (i.e.  non¬ 
absorbing,  u)Q  =  1.0)  and  absorbing  aerosols  (oqo  ^  1.0)  and 
(d)  differing  geometries  (variations  in  observer  and  solar 
orientations) . 


t 
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6.2  Analytical  Approaches 


From  a  practical  perspective,  the  essential  difference 
between  numerical  and  analytical  approaches  is  that  the  lat¬ 
ter  require  very  little  computational  effort  and  hence  a  sig¬ 
nificant  amount  of  computer  time  may  be  saved  if  they  are 
implemented.  Analytical  approaches  are  defined  as  those 
which  are  implementable  using  only  simple  algebraic  equa¬ 
tions  with  specifiable  parameters.  There  are  no  recurrent 
or  iterative  calculations  and  no  convergence  criteria.  Un¬ 
fortunately,  exact  analytical  treatments  are  available  only 
for  a  few  cases  which  are  not  immediately  applicable  to 
geophysical  remote  sensing  problems.  They  are  of  interest 
mathematically,  however,  to  compare  with  corresponding  cases 
of  more  computationally  complex  numerical  techniques.  These 
cases  include,  for  the  most  part,  approaches  based  on  Chan¬ 
drasekhar's  H  functions  for  semi-infinite  atmospheres 
(t  =  »)  and  X  and  Y  functions  for  finite  atmospheres  (Chan¬ 
drasekhar,  1960) .  For  the  former  set  of  problems  solutions 
are  available  for  isotropic  (Stibbs  and  Weir,  1959;  Abhyankar 
and  Fymat,  1971) ,  Rayleigh  (Chandrasekhar,  1960;  Lenoble, 
1970)  and  various  anistropic  phase  functions  (Chandrasekhar, 
1960;  Sobolev,  1956;  Kolesov,  1972).  For  finite  atmospheres, 
solutions  are  available  only  for  isotropic  (Carlstedt  and 
Mullikan,  1966)  and  Rayleigh  scattering  (Sekera  and  Kahle, 
1966)  . 

For  remote  sensing  problems  relevant  to  radiative  trans¬ 
fer  in  the  atmosphere-ocean  system,  methods  are  required 
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which  treat  a  finite  atmosphere  (i.e.,  the  total  optical 
depth  is  not  infinite)  with  scattering  properties  which  are 
anisotropic  in  order  to  simulate  aerosol  scattering.  Although 
exact  analytical  treatments  are  not  available  in  such  cases, 
a  variety  of  approximate  analytical  methods  with  quantifiable 
accuracy  may  be  employed.  The  utility  of  approximate  analy¬ 
tical  treatments  lies  in  their  extreme  computational  effi¬ 
ciency  while  retaining  treatment  of  the  salient  physical 
mechanisms  of  radiative  transfer  involved. 

Approximate  analytical  methods  include  those  based  on 
taking  the  first  few  tractable  orders  of  more  extensive 
numerical  treatments  and  those  formulated  specifically  as 
approximate  treatments.  Examples  falling  into  the  first 
category  include  first  (Deirmendjian,  1969)  and  second  (Ho- 
venier,  1971)  order  of  scattering  treatments  explicitly  for¬ 
mulated  from  successive  order  of  scattering  approaches 
(Irvine,  1965;  Nagel  et  al.,  1978),  analytic  two  and  four 
stream  (Liou,  1974)  approximations  based  on  discrete  ordinate 
methods  (Liou,  1973) ,  and  the  two-step  function  approach 
(Burke  and  Sze,  1977)  derived  from  more  general  variational 
methods  (Sze,  1976) .  Treatments  explicitly  formulated  as 
approximate  approaches  include  hybrid  methods  evaluating 
single  scattering  exactly  and  approximating  multiply  scat¬ 
tered  terms  (Bergstrom  et  al.,  1980;  Hering,  1981)  and  var¬ 
ious  general  finite  stream  methods  (Meador  and  Weaver,  1981) 
when  applied  as  suggested  by  the  work  of  Kaufman  (1979)  . 
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The  inherent  degree  of  accuracy  associated  with  a  spe¬ 
cific  approximate  analytical  treatment  varies  considerably 
with  the  relevant  optical  propagation  parameters  involved. 
Thus  a  domain  of  applicability  must  be  ascertained.  Expec¬ 
ted  errors  may  be  quantified  by  examining  standard  cases  and 
comparing  accuracies  either  with  available  exact  solutions 
or  with  numerical  solutions  of  specified  precision.  For 
example.  Table  6-1  lists  percentage  error  for  a  comparison 
between  upward  and  downward  radiances  derived  from  single 
scattering  vs.  multiple  scattering  (Coulson  et  al.,  1960) 
treatments  in  a  Rayleigh  atmosphere  with  the  sun  at  zenith 
(y^  =  1.0)  and  zero  surface  albedo  (AT  =  0.0). +  Results 
are  presented  for  various  observer  zenith  angles  (y)  and 
total  atmospheric  optical  depths  (x*) .  Note  that  even  for 
this  fundamental  approximate  analytical  approach,  errors 
are  highly  dependent  on  atmospheric  optical  depth,  emphasiz¬ 
ing  the  need  to  quantify  the  behavior  of  such  treatments 
a  priori. 

6 . 3  Adopted  Approach 

The  approximate  solution  to  the  radiative  transfer  equa¬ 
tion  adopted  in  the  current  radiance  simulation  model  is 
based  on  that  suggested  by  Kaufman  (1979)  and  employed  in 
a  variety  of  remote  sensing  contexts  (Kaufman,  1981,  1982; 

+These  parameters  are  defined  in  §6.4. 


Upward  Radiance 


Downward  Radiance 


y 

0.10 

* 

T 

0.25 

0.50 

0.10 

* 

T 

0.25 

0.50 

1.00 

13.2 

24.6 

37.1 

12.3 

24.7 

37.9 

0.72 

13.8 

25.6 

39.3 

13.4 

26.1 

40.3 

0.52 

14.3 

26.9 

40.8 

13.9 

27.2 

42.3 

0.28 

14.7 

27.8 

41.9 

14.6 

28.5 

44.6 

0.10 

14.8 

27.6 

40.2 

14.9 

29.3 

46.1 

Table  6-1.  Percent  error:  single  scattering  vs.  multiple 

scattering  with  Rayleigh  phase  function  (y  =  1.0, 
Al  =  0.0)  .  0 


) 


Kaufman  and  Joseph,  1982;  Mekler  and  Kaufman,  1980,  1982)  . 
The  approximate  solution  for  the  radiance  is  an  extension 
of  existing  methods  for  flux  calculation  (Meador  and  Weaver 
1980)  based  on  the  assumption  that  the  diffuse  incident 
flux  may  be  used  instead  of  the  actual  directionally  depend 
ent  radiance  in  the  calculation  of  the  scattered  radiance 
contribution.  Originally  formulated  using  the  two  stream 
(Coakley  and  Chylek,  1975;  Acquista  et  al.,  1981)  flux 
parameterization,  the  method  is  generally  applicable  to 
any  two  parameter  flux  parameterization  including  the  Ed¬ 
dington,  <5 -Eddington,  and  others  discussed  in  Meador  and 
Weaver  (1981).  The  accuracy  of  derived  radiances  is,  of 
course,  dependent  in  turn  on  the  accuracy  of  the  flux  para¬ 
meterization  used.  By  employing  calculated  backscatter 
fractions  (Wiscombe  and  Grams,  1976)  for  general  aerosol- 
phase  functions  in  the  analysis,  the  general  scattering 
anisotropy  of  the  real  atmosphere-ocean  system  may  be  cal¬ 
culated  . 

6 . 4  Thee 

In  order  to  simulate  the  upward  radiance  at  the  top  of 
the  atmosphere  by  a  meteorological  satellite  viewing  the 
atmosphere-ocean  system,  it  is  necessary  to  solve  the  RTE 
for  a  plane-parallel,  multiple  scattering  atmosphere  of  the 
form  (Chandrasekhar,  1960) : 
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(6-1) 


(x  ,  y  ,  <J>)  =  I  (x  , y ,  <J>)  -  coQJ  (x  ,  y ,  4>) 


where : 


2tt+1 


(t  ,y ,  <j>)  =  j  P(0o)e"T/y°  +  ±  J  j  P  ( 0 )  I  ( x ,  y  '  ,  <J> ' )  dy ' 


0  -1 


(6-2) 


Expression  (6-2)  is  the  scattering  source  function  with  terms 
on  the  right-hand  side  corresponding  to  single  scattered  and 
multiply  scattered  (i.e.,  the  integral  over  angles)  radiances, 
respectively.  The  following  definitions  apply: 


I  (x  ,  y  ,  <J>) 
J(x,y  ,<J>) 


P(9) 

pR(e) 

Pa(9) 


=  radiance  (mW  cm  2  ym  ^  sr  1) 

=  source  function  due  to  scattering  by 
atmosphere 

=  optical  depth,  total  optical  depth  f°°$(z)dz 

=  total  extinction  coefficient  (km  )  at 
height  z  (km) 

=  +  Bg  +  PX  +  BJ 

=  Rayleigh  scattering  coefficient  (km  ) 

=  molecular  absorption  coefficient  (km  ) 

=  aerosol  scattering,  absorption  coefficients 

(km"l) 

=  single  scattering  albedo 

■  (e5  +  ba)/s 

=  angular  scattering  (phase)  function 

=  [ S®  PR(0)  +  PA(0) 1/S 

=  Rayleigh  phase  function 
=  aerosol  phase  function 

=  incident  solar  irradiance  (see  Figure  6-1) 


All  of  the  above  are  wavelength  dependent  optical  properties 
of  the  atmosphere.  The  sun-sensor  geometry  is  defined  by 


the  angles : 


<• 


6-9 


4 


t  /  \ 


solar  flux 
V  S  .  .’t 


e  130 

o 


.3  .9  1  1.1 


Wavelength  (urn) 


Figure  6-1.  Incident  Solar  Irradiance 
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Figure  6-2.  Ocean  Surface  Albedo  (Ramsey,  1968) 
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=  cos0  (8  is  the  sensor  zenith  angle) 

=  cos0Q  (8q  is  the  sun  zenith  angle) 

=  azimuth  angle  difference  between  sensor  and  sun 

=  scattering  angle  for  incident  direction  y  and 
scattered  direction  y'  with  azimuth  angle 
difference  <j> 

=  yy'  +  (l-y^)1*  (l-y'^)1*  cos<J> 


Analytical  radiance  solutions  to  (6-1,  6-2)  are  obtained 

as  described  above  by  assuming  that  the  angular  dependent 

intensity  factor  in  the  integral  sum  of  the  source  function 

may  be  taken  as  isotropic  (i.e.,  independent  of  angle)  to 

first  order  (see  Appendix  A,  A.l).  This  approximates  the 

multiply  scattered  field.  A  quantity  independent  of  angle 

related  to  radiance  is  flux  (or  irradiance)  is  defined  as: 
2tt+1 


F(x)  = 


I  (t  ,y  '  ,  <j) ' )  y  'dy  ’d<J>  ’ 


(6-3) 


0  -1 


=  TTl  (t) 


The  last  step  follows  by  assuming  radiance  is  independent  of 

angle.  The  result  above  (6-3)  is  substituted  back  into  (6-1, 

6-2)  [cf.  A. 8]  to  relate  radiance  to  flux.  Thus  by  solving 

the  resulting  coupled  differential  equations  (A. 9,  A. 10), 

radiance  solutions  are  obtained  dependent  on  the  fluxes. 

The  boundary  conditions  required  for  solution  at  the  surface 
* 

(t  =  t  )  and  top  of  the  atmosphere  (x  =  0)  are  imposed  on 
the  fluxes: 


F  (x  =  0)  =  0 

F+(x  =  x*)  =  Ax  [irFy  e  T  //u°  +  F  (x  =  x*)  ]  .  (6-4) 

Li  0 
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and  state  that  there  is  no  diffuse  irradiance  from  space 
at  the  top  boundary  and  that  the  upward  irradiance  at  the 
surface  is  equal  to  the  product  of  a  constant,  AT  (wave- 

length-dependent) ,  and  the  downward  irradiance  due  to  the 

"f* 

direct  solar  beam  and  the  multiply  scattered  radiance. 

Here  A  is  the  wavelength-dependent  Lambert  surface  albedo 
Li 

*t* 

from  Ramsey  (1968)  [as  given  in  Curran  (1972)].  Ocean 
surface  albedo  as  a  function  of  wavelength  is  illustrated 
in  Figure  6-2. 

Solutions  to  (6-1,  6-2)  for  upward  radiance  are  obtained 
of  the  form: 


.  0.ekT  U2e-kT  U,e~T/lJ° 
kt.u.*)  =  Dle  +  -T+u/u- 


where : 


D  =  r  ii 

(T*) 

i  * 

TT  kT 

Uj^e 

,,  -kT* 

V 

1 

TT 

1-yk 

1+yk 

Aok 

U,  " 

1  TT 

1-6 (y ) 

,  6(u) 
Y2 

i - 1 

M 

1 

i — 1 

Bco 

1-6 (y) 

+  8  (y ) 
y2 

■  (Yj+k)  ] 

Cod 

U-,  =  - — 

3  IT 
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e’T*/y  (6-6) 


(6-7) 


(6-8) 


(6-9) 


Discussion  of  the  specular  (non-Lambert)  surface  reflec¬ 
tance  due  to  sunglint  is  deferred  until  Section  7. 


6-12 


(6-10) 


Y  =  c  (Yi  +  ~)  -  irF  a)  6  (U  ) 

y  o  oo 


,  ,  2  2. 
k  =  (yL  -  Y2) 


1  2  IT 


8  (u)  = 


P(y  ,  y  <j> ,  <£  1 )  dy  '  dc{> ' 


(6-11) 


(6-12) 


Values  of  A,  B,  C  and  the  upward  flux  at  the  surface, 

+  * 

F  (t  )  are  given  by  expressions  in  Appendix  A  (A. 15,  16,  17, 

and  11a,  respectively).  The  backscatter  fractions,  B(y), 

defined  in  (6-12)  above  (see  also  A. 4,  5,  6)  are  obtained 

from  results  of  Wiscombe  and  Grains  (1976)  and  are  dependent 

on  the  phase  function  asymmetry  parameter  defined  as  Joseph 

et  al.  (1976) : 

+1 


g  =  j  j~  yP(y)dy 


(6-13) 


Figure  6-3  illustrates  the  dependence  of  B(y)  and  B '  on  g 
as  used  in  the  radiance  simulation  model.  The  correct 
choice  of  asymmetry  parameter  used  in  the  evaluation  of 
8(y)  is  that  obtained  using  the  effective  phase  function 
(see  6-14)  in  the  definition  of  g  (6-13)  above.  Generally 
B'  varies  from  0.5  for  isotropic  or  Rayleigh  scattering 
(g  =  0)  to  0.0  for  total  forward  scattering  (g  =  1.0) 

[i.e.,  there  is  no  backscatter].  The  range  of  8(y)  values 
is  similar,  however,  the  rate  of  charge  in  the  domain 
0.0  <  g  <  1.0  is  dependent  on  the  specific  value  of 
y(=  cos8)  chosen. 
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Figure  6-3.  Backscatter  Fractions  as  *  Function 

of  Asymmetry  Parameter 
(Wiscombe  and  Grams,  1976) 
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6.4.1  Phase  function  parameterization 

The  phase  function  P(9)  defined  above  appears  explicitly 
in  the  solution  (see  6-9)  where  it  determines  primary  scat¬ 
tering.  As  noted  immediately  above,  the  backscatter  frac¬ 
tions  8(u),  8'  are  also  dependent  on  the  phase  function. 

The  appropriate  phase  function  above  is  the  effective  phase 
function  at  a  given  wavelength  weighted  for  Rayleigh  scat¬ 
tering  and  aerosol  scattering  contributions.  For  a  one 
layer  model  this  is  given  by: 

Peffl0)  "  V0)  +  ta  V0)1/t'  <6-14) 

where : 

XR  =  6R  Hr  (6-15) 

tA  =  eA  HA  (6-16) 

PR ( 9 )  =  |  (1  +  cos2e)  (6-17) 

and  H  ,  H  are  molecular  and  aerosol  scale  heights,  respec- 
tively.  Figure  6-4  illustrates  the  variation  of  the  effec¬ 
tive,  weighted  phase  function  at  0.55ym  wavelength  (where 
s 

tr  =  0.1)  with  increasing  aerosol  optical  depth.  The  aero¬ 
sol  phase  function  P^O)  is  the  Shettle  and  Fenn  (1979) 
maritime  aerosol  model  at  70%  relative  humidity.  A  question 
naturally  arises  concerning  the  degree  of  accuracy  required 
in  specifying  the  aerosol  phase  function  in  order  to  provide 
useful  radiance  simulations.  Two  alternatives  were  investi¬ 
gated:  (a)  exact  phase  function  digitized  at  9  increments 


<0 
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Angle  (deg) 


Change  of  weighted  atmospheric  phase 

function  at  0.55ym  wavelength  with 

increasing  aerosol  optical  depth. 

Illustrated  are  Raleigh  phase  function 

(....)  t  :  0.1,  x  =  0.0  and  curves  for 
r  a 

t  =  0.1,  0.5,  1.0,  and  5.0  (in  descend 

cl 

ing  magnitude  order  at  120°)  . 


sufficiently  dense  to  characterize  the  shape  and  (b)  the 
analytical  Henyey-Greenstein  (1941)  phase  function  (H-G) 
given  by 

P(0,g)  =  (l-g2)/(l  +  g2  -  2g  cos0)  3/"2  (6-18) 

Figure  6-5  illustrates*  a  comparison  between  exact  and  H-G 
phase  functions  for  the  Shettle  and  Fenn  (1979)  maritime 
aerosol  at  relative  humidities  of  70%  and  80%  within  the 
DMSP  spectral  interval;  (a)  A  =  0.55ym,  (b)  A  =  0.69ym,  and 
(c)  A  =  1.06ym.  In  the  context  of  model  parameterization 
it  is  noted  that  the  apparent  dependence  on  relative  humi¬ 
dity  is  greater  than  that  on  wavelength.  This  is  illustrated 
by  Figures  6-6  and  6-7.  Figures  6-6a  and  b  illustrate  the 
exact  maritime  phase  functions  at  relative  humidities  of 
70%  and  80%,  respectively.  All  three  wavelengths  are  plotted 
in  each  and  it  is  obviously  difficult  to  distinguish  three 
distinct  curves.  (There  are  certainly  differences  in  the 
forward  (9  =  0°)  and  backscatter  (9  =  180°)  directions 
which  may  be  important  in  some  cases.)  In  contrast,  Figure 
6-7  illustrates  the  exact  phase  functions  at  0.55ym  wave¬ 
length  for  70%  and  80%.  Considerable  difference  is  noted 
especially  at  intermediate  scattering  angles.  The  conclu¬ 
sion  (which  may  also  be  apparent  from  the  respective  g 
factors  provided  in  Figure  6-5)  is  that  for  a  given  rela¬ 
tive  humidity/aerosol  model,  one  phase  function  may  suffice 
for  all  wavelengths  within  the  chosen  spectral  interval. 
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Log  (Phase  Function)  Log  (Phase  Function)  Log  (Phase  Function) 


Figure  6-5.  Comparison  of  exact  (-)  and  approximate  H-G(...) 

phase  functions  for  maritime  aerosol  at  wave¬ 
lengths  of:  (a)  0.55um,  (b)  0.69um,  and  (c) 
1.06ym  and  relative  humidities  of  70%  and  80%. 
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Figure  6-6.  Maritime  aerosol  phase  function  for  relative 
humidities  of:  (a)  70%  and  (b)  80%.  Plotted 
are  values  at  wavelengths  of  0.55,  0.68,  1.06um 
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Figure  6-7.  Maritime  aerosol  phase  function  at  wavelength 
of  0.55ym  for  relative  humidities  of:  70%  (+) 
and  80%  (o)  . 
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It  is  also  apparent  from  Figure  6-5  that  the  analytical 
H-G  parameterization  does  not  reproduce  the  exact  phase 
function  in  the  forward  and  backscatter  directions  although 
it  does  well  in  the  approximate  range  from  40°  <  0  < 140°. 

One  possible  remedy  (Hering,  1981)  is  adoption  of  a  two 
term  H-G  phase  function  (TTHG)  (Irvine,  1968)  given  by  the 
linear  combination  of  single  term  H-G  phase  functions  (STHG) : 

P(9/g1,g2/c)  =  c  Ph-G^'5!*  +  (1-c)  PH-G(9,g2)  (6-19) 

Figure  6-8  compares:  (a)  an  exact  continental  aerosol 
scattering  function  (from  Selby  and  McClatchey,  1972) , 

(b)  a  suggested  TTHG  phase  function  (Hering,  1981)  with 
g^  =  .714,  g2  =  -.613,  and  c  =  .963,  (c)  the  corresponding 
STHG  [i.e.,  overall  g  factor  same  as  (c)  above  g  =  cg^  + 

(1-c)  g2  =  .664]  and  (d)  a  STHG  used  in  Liou  (1973)  to 
characterize  aerosol  scattering  with  g  =  0.62.  The  TTHG 
phase  function  (b)  does  a  much  better  job  at  providing  the 
overall  shape  of  a  typical  exact  phase  function  than  does 
either  of  the  STHG  parameterizations .  In  order  to  investi¬ 
gate  the  impact  of  the  phase  function  parameterization  on 
simulation  of  upward  radiances  using  the  approximate  stream 
method  (SA)  described  above,  calculations  were  performed 
for  a  range  of  optical  depths  (0.125  to  1.0)  and  sensor 
zenith  angles  from  0  to  75°  for  an  overhead  sun  (uq  =  1.0) 
using  both  the  STHG  (curve  d,  Figure  6-8)  and  TTHG  (curve 
b,  Figure  6-8).  For  optical  depth  0.125  a  single  scattering 
(SS)  calculation  was  done  using  the  TTHG  phase  function 


(single  scattering  would  not  be  appropriate  for  larger  opti 
cal  depths).  The  single  scattering  albedo,  goq  ,  for  these 
calculations  was  0.987.  Results  are  illustrated  in  Figure 
6-8 .  The  single  scattering  result  uniformly  underestimates 
the  exact  radiance  due  to  the  neglect  of  multiply  scattered 

contributions.  The  best  fit  overall  (i.e.,  for  all  t  and 
0)  is  obtained  using  the  TTHG.  The  STHG  particularly  under 
estimates  backscattering  (in  this  case,  the  vicinity  of 
y  =  1.0)  due  to  its  failure  to  provide  the  backscatter  lobe 
of  the  exact  phase  function  (see  Figure  6-8) .  For  scatter¬ 
ing  angles  in  the  vicinity  of  120°  (zenith  angles  of  60°) 
both  STHG  and  TTHG  provide  similar  results.  This  is  as 
expected  since  the  two  parameterizations  cross  over  near 
this  angle.  Coincidentally  a  nadir  viewing,  sun  synchro¬ 
nous  (morning  or  afternoon)  satellite  sensor  will  have  a 
scattering  geometry  with  backscatter  angle  in  this  vicinity 
Figure  6-10  compares  simulated  wavelength  dependent  DMSP 
radiances  using  exact  and  STHG  (denoted  P  and  HG,  respec¬ 
tively)  phase  functions  assuming  the  70%  maritime  model. 

The  scattering  geometry  corresponds  to  about  110°  backscat¬ 
ter.  Note  that  the  radiance  spectra  and  bandpass  weighted 
DMSP  radiances  are  virtually  identical.  As  illustrated 
here  (and  also  in  Hering,  1981)  the  TTHG  gives  a  better 
overall  fit  to  the  angular  scattering  function  and  result¬ 
ant  radiances.  However,  for  the  purpose  of  DMSP  simulation 
within  the  present  context  a  STHG  approach  may  be  adequate. 
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Figure  6-9 


Comparison  of  Stream  Approx.  Intensities  Using 
Single  ( STHG)  and  Two  Term  (TTHG)  Henyey-Greenstein 
Phase  Functions 


9 

0, 


Zenith  Angle  0  (y  =  cos  0) 


Log  (Radiance) 


Figure  6-10.  Comparison  of  wavelength  dependent 
radiances  calculated  using  the  com¬ 
plete  phase  function  P(0)  and  the 
corresponding  approximate  H-G  phase 
function  for  the  maritime  aerosol 
model  at  70%  relative  humdity. 
Corresponding  simulated  DMSP  radi¬ 
ances  are  0.6  985  and  0,7054  mW  cnT2 
sr”i'  respectively. 
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6.4.2  Flux  Parameterization 

As  noted  previously,  any  two  parameter  flux  parameteri¬ 
zation  may  be  utilized  in  the  evaluation  of  the  radiances 
using  (6-5,  6-12).  The  specific  parameterization  used  en¬ 
ters  the  analysis  through  the  values  Y^  and  y2  in  the  expres¬ 
sions  cited  above.  Table  6-2  from  Meador  and  Weaver  (1980) 
presents  two  of  the  possible  choices.  As  noted  the  Appen¬ 
dix  is  based  on  the  hemisphere  constant  (HC)  parameterization 
which  was  used  by  Kaufman  in  his  original  paper  (1979)  . 

Since  the  validity  of  the  radiance  approximation  depends  on 
the  degree  to  which  fluxes  are  accurately  calculated,  various 
applicable  flux  parameterizations  described  in  Meador  and 
Weaver  (1980)  were  tested  in  comparison  to  accurate  numeri¬ 
cal  calculations.  Table  6-3  compares  reflection  (R)  and 
transmission  (T)  functions: 

r  =  f+(t  =  0)/yoTrF  (6-20) 

x*,y0  * 

T  =  e  /0  +  F(t  =  t  )/uqttF  (6-21) 

from  accurate  discrete  ordinate  method  (E)  calculations 
(Liou,  1973)  to  analytic  flux  parameterizations  including: 
Liou's  (1974)  two  stream  (TS) ,  the  hemisphere  constant  (HC) 
approach  of  Coakley  and  Chylek  (1975)  and  the  hybrid  modi¬ 
fied  Eddington-delta  function  (HM6E)  (Meador  and  Weaver, 

1980)  (see  Table  6-2).  Calculations  assume  g  =  0.75, 
y  =  1.0  and  total  optical  depths  of  0.25  and  1.0,  single 
scatter  albedos  of  0.8  and  1.0,  and  solar  zenith  angle 
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Table  6-2.  Coefficients  and  y2  determining  two  stream  flux  parameterization  (from 
Meador  and  Weaver,  1980)  . 


cosines  of  0.1,  0.5,  and  0.9.  The  R  and  T  functions  determine 
the  upward  and  total  downward  fluxes  given  by  a  particular 
flux  calculation.  Examining  the  values  in  Table  6-3  the 
following  conclusions  are  drawn:  (1)  the  TS  parameterization 
is  the  least  reliable  sometimes  predicting  negative  values 
for  the  smaller  optical  depth,  (2)  the  HC  approach  is  more 
accurate  for  the  smaller  optical  depth  (it  is  formulated  to 
have  the  correct  limit  for  small  optical  depths)  than  either 
of  the  others,  however,  it  is  less  accurate  by  far  than  the 
HM6E  for  the  larger  optical  depth,  and  (3)  overall  the  HM6E 
gives  the  best  representation  for  the  domain  illustrated. 

Based  on  these  conclusions  the  hybrid  modified  Eddington 
delta  function  flux  parameterization  of  Meador  and  Weaver 
(1980)  [see  Table  6-2]  was  adopted  for  the  overall  radiance 
simulation  model.  In  particular  some  inacurracies  may  be 
expected  for  small  optical  depth  situations,  however,  these 
cases  correspond  to  high  visual  ranges.  Since  from  a  prac¬ 
tical  perspective  it  is  probably  more  important  to  have  higher 
operational  accuracy  for  low  visual  range  (i.e.,  higher  opti¬ 
cal  depth)  cases  this  seems  an  appropriate  choice.  Future 
refinements  may  include  an  option  to  switch  to  the  HC  para¬ 
meterization  for  increased  accuracy  at  small  optical  depths. 

6.5  Comparison  to  Exact  Results 

The  approximate  solution  to  the  RTE  described  above 
based  on  Kaufman  (1979)  was  compared  to  a  variety  of  exact 
treatments. 
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1.0 

0.41610 

0.07179 

0.02250 

0.5839 

0.9282 

0.9775 

E 

0.41433 

0.07635 

0.01294 

0.5887 

0.9237 

0.0129 

TS 

0.37221 

0.0681 

0.0255 

0.6102 

0.9279 

0.9721 

HC 

0.38574 

0.0634 

0.0198 

0.6141 

0.9365 

0.9802 

HM6E 

0.8 

0.28961 

0.04855 

0.01547 

.4302 

.8476 

.9266 

E 

0.31802 

0.05739 

-.01125 

.4657 

.8498 

.9540 

TS 

.28416 

0.04976 

0.017861 

.4800 

.8530 

.9210 

HC 

.28567 

0.046186 

0.01421 

.4739 

.8588 

.9298 

1 

HM6E 

1.0 

0.58148 

0.24048 

0.09672 

0.5839 

0.7595 

0.9033 

E 

0.51962 

0.22559 

0.02389 

0.4822 

0.7744 

0.9761 

TS 

0.51753 

0.2392 

0.1381 

0.4648 

0.7575 

0.8567 

HC 

0.5746 

0.2026 

0.0848 

0.4250 

0.7971 

0.9150 

HM6E 

0.8 

0.35487 

0.12342 

0.04929 

.2055 

.5161 

.7177 

E 

0.37519 

0.14279 

-.00064 

.2902 

.5527 

.7633 

TS 

0.36987 

.14511 

.07754 

.2695 

.5280 

.6739 

HC 

0.36111 

.11606 

.04726 

.21540 

.56020 

.7321 

HM6E 

Table  6-3.  Comparison  of  two  stream  (TS)  ,  hemisphere  con¬ 
stant  (HC) ,  and  hybrid  modified  Eddington  delta 
function  (HM6E)  flux  parameterizations  to  exact 
(E)  discrete  ordinate  method  calculations  (from 
Liou,  1973) . 


6.5.1  Isotropic  scattering 

Isotropic  scattering  is  defined  by  the  angular  scatter¬ 
ing  function: 

P  ( 6)  =  o>0  (6-22 

i.e.,  there  is  no  angular  dependence.  While  this  is  not 
applicable  to  the  realistic  atmosphere  exact  results  exist 
for  comparison.  Table  6-4  summarizes  a  comparison  between 
the  two-steam  model  and  exact  results  for  isotropic  scatter¬ 
ing  obtainable  using  Chandrasekhar ' s  X  and  Y  function  tabu¬ 
lated  in  Carlstedt  and  Mullikin  (1966) .  These  calculations 
were  done  using  the  HC  flux  parameterization  and  assume 
<*>0  =  0.5,  yQ  =  1.0.  It  is  particularly  notable  that  errors 
are  less  than  5%  over  much  of  this  domain,  and  particularly 
for  emergent  intensities  with  zenith  angles  approaching 
unity  (y-*1.0)  .  This  geometry  simulates  a  nadir-pointing 
satellite. 

6.5.2  Rayleigh  Scattering 

Analogous  results  for  Rayleigh  scattering  are  provided 
in  Table  6-5.  Rayleigh  scattering  characterizes  the  pure 
molecular  atmosphere  without  aerosols.  In  this  case,  the 
solar  zenith  angle  is  fixed  at  57°  and  the  observer  zenith 
at  13.5°  to  simulate  a  satellite  field  of  view  for  a  polar- 
orbiting  sensor.  Percent  errors  are  given,  comparing  two- 
stream  results  to  the  exact  calculation  of  Rayleigh  scatter¬ 
ing  by  Coulson  et  al.  (1960)  for  a  variety  of  optical  depths 


Upward  Radiance 


Downward  Radiance 


y 

4.8  1.8  2.5  1 

0.2 

-14.6  2.7  2.7 

0.6 

-24.8  -4.6  1.8 

5.0  2.6  .2 

•  0.9 

1 

-  3.0  .1  1.4 

4.9  2.7  -  .3 

Table  6-4.  Percent  error:  Two-stream  approximation  vs.  exact 
result  for  isotropic  scattering  evaluated  from 
Chandrasekhar's  X  and  Y  functions.  (u>  =  0.5, 

=  1.0) 


Table  6-5. 


Percent  error:  Two-stream  approximation  vs. 
exact  result  for  Rayleigh  scattering. 

(u  =  .98,  u  =  .60) 
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surface  albedos  of  0.0  and  0.25,  and  azimuth  distances  of 
30°,  90°,  and  150°.  A  comparison  of  these  results  to  the 
single  scattering  approximation  results  in  Table  6-1,  indi¬ 
cates  that  a  much  higher  degree  of  accuracy  is  achievable 
using  the  two-stream  approach.  Further  examples  are  provi¬ 
ded  in  Isaacs  (1981). 

6.5.3  Aerosol  Scattering 

The  most  stringent  test  of  an  approximate  multiple¬ 
scattering  radiative  transfer  model  is  highly  anisotropic 
aerosol  scattering.  Using  the  hemispherical  constant  flux 
parameterization  Kaufman  (1979)  has  estimated  the  error  for 
a  nadir  viewing  satellite  to  be  a  few  percent.  Calculations 
based  on  the  implementation  of  the  approach  described  above 
have  achieved  errors  between  5%  and  16%  in  radiance  for 
typical  cases. 

6. 5. 3.1  Optically  thin  case 

Table  6-6  summarizes  a  comparison  between  exact  Gauss- 
Seidel  (Dave  and  Gazdag,  1970)  calculations  and  the  adopted 
approach  for  an  optically  thin  (t*  =  .11026)  situation 
taken  from  Dave  (1972) .  The  effective  phase  function  is 
illustrated  in  Figure  6-11  along  with  the  H-G  parameteriza¬ 
tion  used  in  the  approximate  calculation  (g  =  .6862).  Two 
values  of  the  surface  albedo  (0.0,  0.4)  were  investigated. 
Comparisons  of  exact  and  approximate  radiances  are  illustra¬ 
ted  in  Figures  6-12,  13,  and  14.  Upward  radiance  errors 
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Table  6-6.  Comparison  of  upward  and  downward  diffuse  fluxes 
F+ (o) ,  F  (x*)  and 
upward  total  flux  F+  (x*)  for 
various  parameterizations . 


(Figure  6-12)  for  zero  surface  albedo  (r  =  0.0)  and  nadir 
viewing  (zenith  angle  of  0°)  are  about  5%.  This  level  of 
error  is  consistent  with  Rayleigh  scattering  results  (Table 
6-5)  since  in  backscatter  the  Rayleigh  phase  function  domi¬ 
nates  for  low  aerosol  loading  (see  Figure  6-4) .  Error  is 
minimal  for  the  larger  surface  albedo  since  reflected  radi¬ 
ance  is  significant.  Maximum  errors  (20%)  are  evident  in 
the  forward  scattering  direction  (Figure  6-13,  9  =  60°). 

6. 5. 3. 2  Optically  thick  cases 

Comparisons  for  three  optical  depth  unity  cases  from 
Lenoble  (1977)  are  summarized  in  Table  6-7.  Flux  parame¬ 
terization  results  are  compared  with  exact  spherical  har¬ 
monic  results.  The  superiority  of  the  HMSE  approach  for 
larger  optical  depths  is  confirmed  for  cases  1  and  2  repre¬ 
senting  overhead  sun.  For  case  3  (pq  =  0.5),  results  are 
somewhat  ambiguous.  The  exact  phase  function  and  its  HG 
counterpart  are  illustrated  in  Figure  6-15.  Figures  6-16 
and  6-17  illustrate  radiance  comparisons  with  the  exact 
calculations  for  three  approximate  approaches:  (a)  HM6E 
flux  parameterization,  exact  phase  function;  (b)  HC  flux 
parameterization,  HG  phase  function;  and  (c)  HC  flux  para¬ 
meterization,  exact  phase  function.  For  cases  1  and  2 
(Figure  6-16) ,  the  first  of  these  approaches  does  best 
overall.  When  the  sun  is  not  at  zenith  (case  3,  Figure 
6-17)  this  approach  is  superior  only  for  nadir  viewing. 
Illustrated  are  cases  for  azimuth  angles  of  0,  90,  and 


(a) 

Case 

1 

T*  = 

1.0 

.1736 

1.8124 

Spherical  Harmonics'*' 

y .  = 

1.0 

.  3279 

1.6580 

2 

Hemi-const 

ii 

o  o 

3 

1.0 

.1615 

1.8240 

H-Mod- 6-Edd3 

(b) 

Case 

2 

t*  = 

1.0 

.1236 

1.5155 

Spherical  Harmonics 

= 

1.0 

.2362 

1.3470 

Hemi-const 

0)  = 

0 

0.9 

.1202 

1.5324 

H-Mod- 6-Edd 

(c) 

Case 

3 

T*  = 

1.0 

.2255 

.8032 

Spherical  Harmonics 

Uo  = 

0.5 

.2403 

.8327 

Hemi-const 

U)  = 

0 

0.9 

.1962 

i _ 

.8803 

H-Mod- 6-Edd 

1Lenoble  (1977) 

2Kaufman  (1979) 

3 Meador  and  Weaver  (1980) 


Table  6-7.  Comparison  of  upward  and  downward  diffuse 
fluxes  for  Haze  L  cases  using  various 
parameterizations . 
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Figure  6-16.  Comparison  of  Methods  for  Case  1 
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Figure  6-16 .Comparison  of  Methods  for  Case  2 
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180°  (3a,  3b,  3c,  respectively).  The  accuracies  for  nadir 
viewing  for  cases  1,  2,  and  3  are  12,  16,  and  27%,  respec¬ 
tively. 


6. 5. 3. 3  Comparison  to  Look-up  Table  Kesults 

Griggs  (1975,  1978a, b,  1979a, b)  has  performed  extensive 
multiwavelength  calculations  based  on  the  highly  accurate 
Dave  (1972)  radiative  transfer  code  to  relate  aerosol  load¬ 
ing  to  remotely  sensed  radiances  for  a  variety  of  sensors. 

A  look-up  table  was  formulated  (Griggs,  1978b)  to  provide 
radiance  values  given  the  following  input:  (1)  sensor 
zenith  angle,  (2)  azimuth  angle  difference,  (3)  sun  zenith 
angle,  and  (4)  aerosol  loading  factor.  This  latter  value 
is  related  to  the  aerosol  optical  depth  through  the  aerosol 
model  properties  (Griggs,  1978a).  Assuming  the  aerosol  model 
illustrated  .in  Figure  6-11,  surface  albedo,  AL  =  0.015,  and 
single  scattering  albedo,  o)Q  =  1.0,  approximate  radiances 
were  evaluated  (using  the  HM5E  flux  parameterization) .  A 
nadir  viewing  sensor  (y  =  1.0)  was  chosen  and  calculations 
were  performed  for  aerosol  loading  factors  of  N  =  1,2, 4, 6, 
8,10  corresponding  to  0.65ym  aerosol  optical  depths  of 
(Elterman,  1970) : 

T,  (0.65)  =  0.2318N  (6-23) 

The  required  Rayleigh  optical  depth  of  xR  =  .05  was  chosen. 

-2  -1  -1 

Radiance  (mW  cm  sr  ym  )  results  for  exact  (E)  and 

approximate  (A)  radiances  (assuming  an  incident  irradiance 
-2  -1 

of  150  mW  cm  ym  at  0.65ym)  and  percent  error  (A)  figures 
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are  provided  in  Table  6-8  for  a  variety  of  solar  zenith  angles. 
The  exact  results  are  apparently  reproduced  to  within  10%  given 
that  the  assumed  input  variables  are  appropriate. 


6.6  DMSP  Bandpass  Weighted  Radiances 

The  solution  to  equation  (6-1,  6-2)  for  each  set  of  me¬ 
teorological  variables  investigated  provides  a  wavelength- 
dependent  intensity  spectrum,  I  (X) .  The  intensity  measured 
by  the  DMSP  sensor  at  the  satellite,  I  (DMSP) ,  however,  is 
given  by  weighting  these  monochromatic  intensities  by  the 
DMSP  spectral  bandpass  function  given  in  Figure  6-18  accord¬ 
ing  to : 


I  (DMSP) 


4>  (X)  I  (X)  dX 


AX 


<HX)dX 


(6-24) 


where  cp  is  the  sensor  response  function  over  wavelength  in¬ 
terval  AX  (0.4-l.lym  for  DMSP  VHR,  LF) . 
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1.06  0.97  1.45  1.33  2.01  1.93  2.41  2.49  2.74  3.00  3.04  3.45 

8.5  8.3  4.0  3.3  9.5  13.5 

0.69  0.61  0.84  0.72  1.01  0.97  1.15  1.22  1.28  1.45  1.39  1.65  R 

13.1  14.3  4.0 


Table  6-8.  Comparison  of  approximate  (A)  and  exact  (E) 

-2  -1  -1 

radiances  (R)  in  mW  cm  sr  ym  and 
resulting  percentage  errors  (-A)  in  %. 
results  from  Griggs  (1978b) . 
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7.  REFLECTION  FROM  THE  OCEAN  SURFACE 


4# 


4# 
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7 . 1  Sunglint 

Satellite-borne  visible  imaging  radiometers  receive 
reflected  solar  radiation  from  the  earth's  surface.  Due  to 
the  particular  constraints  imposed  by  their  intrinsic  observ¬ 
ing  geometries,  such  sensors  see  a  relatively  narrow  solid 
angle  field  of  view.  In  order  to  interpret  these  measure¬ 
ments,  therefore,  some  understanding  of  the  directional 
reflectance  properties  of  the  surface  is  required  (Tanre 
et  al. ,  1979).  For  many  natural  surfaces,  it  is  often  pos¬ 
sible  to  assume  that  to  a  first  approximation  radiation  is 
reflected  uniformly  in  all  directions  (or  isotropically) . 

This  Lambertian  reflection  is  particularly  useful  for  a 
variety  of  land  surfaces  although  exceptions  may  be  noted 
(cf.  Eaton  and  Dirmhim,  1979).  The  most  notable  natural, 
non-Lambertian  surfaces  are  bodies  of  water.  Under  appro¬ 
priately  calm  conditions,  water  surfaces  (including  for 
example,  oceans,  lakes,  rivers,  swamps,  etc.)  obey  Fresnel's 
reflection  laws  and  produce  mirror-like  images  of  the  sun 
for  properly  oriented  sensors.  As  discussed  in  Fett  and 
Mitchell  (1977) ,  for  example,  such  sunglint  (or  "sun-glitter") 
is  regularly  observed  in  the  visible  imagery  obtained  from 
both  geostationary  and  polar  orbiting  satellites.  Since 
the  surfaces  of  natural  bodies  of  water,  however,  are  usually 
rippled  and  roughened  due  to  wind  stress  and  other  factors, 


a  perfectly  specular  reflection  is  rarely  observed.  Rather, 
the  sunglint  pattern  is  modified  by  the  existing  meteorologi¬ 
cal  conditions  affecting  the  surface  and  some  directional  re¬ 
flectance  intermediate  between  the  Lambertian  and  specular  is 
often  evidenced.  The  size,  shape,  and  intensity  of  the  result¬ 
ant  sunglint  pattern  thus  depends  on  a  variety  of  factors 
including  (Fett  and  Mitchell,  1979):  (1)  sun-sensor  geometry, 

(2)  characteristics  of  the  image  forming  system,  (3)  sea  state 
(i.e.  wind  speed)  and  (4)  the  low  level  distribution  of  atmo¬ 
spheric  aerosols  and  moisture.  A  number  of  investigators  have 
exploited  this  environmental  modification  of  surface  reflection 
properties  in  order  to  aid  in  the  inference  of  surface  wind 
conditions  (McClain  and  Strong,  1979;  Wylie  et  al.,  1981),  wind 
speed  and  direction  (Kornfield,  1974) ,  and  surface  ridge  lines 
(Anderson,  1974) .  Additional  meteorological  analyses  potentially 
available  from  sunglint  region  data  for  solar  orbiting  satel¬ 
lites  are  discussed  in  Fett  and  Mitchell  (1977) . 

In  this  section,  a  simplified  approach  is  described  to 
incorporate  sunglint  phenomenon  within  the  context  of  the  cur¬ 
rent  DMSP  imagery  simulation  model.  Due  to  the  approximate 
nature  of  the  radiance  calculation  discussed  in  Section  6,  it 
is  not  possible  to  treat  such  effects  in  a  vigorous  manner. 

The  methodology  adopted,  however,  provides  an  excellent  quali¬ 
tative  description  based  on  comparison  with  numerical  results 
and  should  provide  practical  guidance  for  imagery  analysis 


purposes . 


7.2  Surface  Modeling 

As  discussed  above,  static  representations  of  wavelength- 
dependent  ocean  surface  reflectances  based  on  simple  Lamber¬ 
tian  treatments  are  not  necessarily  appropriate  due  to  the 
dynamic  response  of  the  ocean  surface  to  meteorological  forc¬ 
ing  by  wind/wave  interaction.  Three  distinct  considerations 
may  be  cited  as  necessary  for  physically  realistic  simulation 
of  the  reflection  of  radiation  from  the  ocean  surface:  (1) 
geometric  aspects  for  calm  sea  conditions  obtainable  from 
Fresnel's  Laws,  (2)  modifications  due  to  wind  induced,  time- 

i 

dependent  wave  slope  geometries,  and  (3)  incorporation  of  the 
reflected  radiation  within  the  adopted  radiative  transfer 
algorithm. 

In  the  basic  calm  state,  the  sea  surface  may  be  treated 
as  a  fundamental  dielectric  medium  according  to  the  classical 
Fresnel  relations  (Towne,  1967) .  For  given  geometries  between 
the  direction  of  propagation  and  the  normal  to  the  interface 
these  relationships  provide  the  reflection  and  transmission 
coefficients  for  an  incident  plane  wave  of  arbitrary  polari¬ 
zation  and  thus  the  required  explicit  bidirectional  depend¬ 
ence  of  surface  reflectance.  This  reflectance  is  fully 
specular,  however.  The  distinction  between  Lambert  and 
specular  surface  reflection  properties  may  be  clarified  by 
defining  the  bidirectional  reflection  distribution  function, 
f  (u,4>;y'  ,4)' )  .  The  reflected  solar  radiance,  Ir,  in  direc¬ 
tion  (y',4)')  from  a  surface  with  reflectance  characterized 
by  f  is : 


2  TT  1 


Ir (v*  f  4>' )  = 

0  0 

where  the  radiance  directed  toward  the  surface  is : 

F 

I  (y  ,4>)  =  -f-  <S(y-yQ)  6(4)  (7-2) 

and  the  surface  incident  solar  flux  is : 

Fc  =  yotrFe“T*/uo  (7-3) 

For  a  Lambert  surface,  f  is  isotropic  (independent  of  angle), 
thus : 

r  al 

f£  (v,<fr;u'  ,<f>')  «  Al  =  albedo  (7-4) 

Substituting  into  Equation  (7-1)  above  and  integrating  yields 
the  relationship 

I(u,<i>)  =  A  F  /tt  (7-5) 

r  L  o 

or  the  ratio  of  reflected  flux  (ttIj.)  to  incident  flux  (Fq) 

is  the  albedo,  AT .  Except  for  the  diffuse  (or  multiply  scat- 

h 

tered)  term  (i.e.  F-)  neglected  in  the  definition  of  I(y,<}>) 
above,  this  is  the  boundary  condition  used  in  the  analysis  of 
non-sunglint  regions  (see  Isaacs,  1980,  Equation  38): 

F+(t  =  x*)  =  AL[7TFyQ  e~T*/Vo  +  F”(x  *  T*)l  (7-6) 

For  a  specular  surface,  the  bidirectional  reflection 
distribution  function  is: 

f£(y,<f>;y'  ,*’)  =~^-6(u-u')  S(4>'-$+*)  C7 -7) 


ufr(y,4>;u' 


I  (y  ,4>)  dyd<() 


(7-1) 
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where  p(y)  is  the  Fresnel  reflectivity  for  incidence  angle, 
cos  ^y.  Substitution  into  (6-1)  above  yields  a  reflected 
radiance  of: 


(  p<Wyo 


cr(y  ,<P  )  *  < 


for  (y ,  <j>)  =  (yo,4>o+ir) 


\  for  (other  y,<j>) 


(7-8) 


Thus  there  is  a  mirror  (or  specular)  reflection  at  observa¬ 
tion  angle  cos  ^"y  =  cos  ^yQ  and  looking  toward  the  sun. 

This  occurs  in  satellite  imagery  at  the  primary  specular 
point  (PSP)  which  falls  between  the  satellite  subpoint  and 
the  solar  subpoint.  Methods  for  locating  the  PSP  are  de¬ 
scribed  in  Fett  and  Mitchell  (1979)  and  Tsui  and  Fett  (1980) . 

The  nature  of  the  problem  changes  considerably  for  wind- 
ruffled  seas  (cf .  Burt,  1954) .  In  this  case,  wind/wave  inter¬ 
actions  produce  a  sea  surface  slope  distribution  which  may 
be  characterized  explicitly  (Cox  and  Munk,  1954)  or  spectrally 
(Pierson  and  Moskowitz ,  1964;  Pierson  and  Stacy,  1973).  The 
resultant  wave  slope  geometry  modifies  in  a  statistical  sense 
the  relationship  between  the  direction  of  incident  radiation 
and  the  normal  to  the  time-averaged  calm  surface.  Exam¬ 
ination  of  relevant  imagery  for  these  cases  indicates  a  sun- 
glint  region  rather  than  a  specular  point  and  the  surface 
reflectivity  may  best  be  characterized  as  intermediate 
between  Lambert  and  specular.  One  potential  approach  to 


7-5 


treating  the  reflected  radiance  using  Equation  (7-1)  is  to 
define  an  effective  bidirectional  distribution  function, 
f  which  is  a  linear  combination  of  Lambert  and  specular 
functions)  i.e.: 

fE  =  ctfL  +  (1-a)  fE  (7-9) 

r  r  r 

where  ct  is  some  constant. 

An  alternative  approach  consists  of:  (1)  determining 
the  surface  orientation  (i.e.,  crosswind  and  upwind  slopes) 
necessary  to  produce  a  specular  reflection  or  glint  given 
the  sun/sensor  geometry,  (2)  evaluating  the  angle  of  inci¬ 
dence  for  this  surface  orientation  in  the  required  rotated 
(from  a  calm  sea)  system,  (3)  calculating  the  Fresnel  reflec 
tivity  for  this  angle  of  incidence  using  the  optical  proper¬ 
ties  of  sea  water,  and  (4)  finding  the  probability  of  this 
surface  orientation  for  the  given  wind  conditions  based  on 
the  Cox-Munk  sea  slope  probability  distribution  function 
(PDF) .  The  effective  specular  reflectivity  is  then  defined 
as  the  product  of  the  computed  Fresnel  reflectance  for  the 
given  geometry  and  the  probability  of  the  required  sea  sur¬ 
face  orientation  for  the  given  wind  condition.  A  variety 
of  investigators  have  used  variations  of  this  approach  to 
incorporate  the  effects  of  reflection  from  the  rough  ocean 
surface  (Mullamaa,  1964;  Raschke ,  1972;  Fowler  et  al.,  1977; 
Takashima  and  Takayama ,  1981;  Ahmad  and  Fraser,  1982). 

The  formulation  used  in  this  work  is  based  on  identical 
physical  principles  but  is  much  more  heuristic  in  implementa 
tion  providing,  however,  an  analytical  result.  Details  of 


the  approach  and  relevant  results  are  included  in  the  follow¬ 
ing  sections. 

7.3  Wind  Driven,  Rough  Ocean  Surface  Reflection  Geometry 
As  previously  discussed  (Section  5) ,  the  location  of 
particular  satellite  radiometer  resolution  elements  (pixels) 
on  the  earth's  surface  (i.e.,  along  the  satellite  subtrack 
and  associated  scan  lines)  is  fixed  by  spacecraft  orbital 
parameters.  Consequently,  the  geometric  orientation  of 
sensor/ surface  pixel/sun  is  predictable  from  spacecraft  and 
solar  ephemeris  data.  For  a  hypothetical  calm  sea,  only 
specific  points  along  the  orbit  will  allow  for  specular 
reflection.  For  a  wind  driven,  rough  ocean  surface,  however, 
it  is  apparent  that  orientations  are  probable  which  permit 
reflection  of  specular  "glints"  to  the  sensor.  This  requires 
that  the  surface  is  oriented  such  that:  (1)  the  normal  to 
the  surface  is  coplanar  with  both  the  sun  and  sensor  position 
vectors  (i.e.  in  the  principal  plane)  and  (2)  the  normal 
bisects  the  angle  formed  by  the  intersection  of  these  posi¬ 
tion  vectors.  Given  the  sun/sensor  position  vectors  relative 
to  a  calm  ocean  surface,  it  is  possible  to  evaluate  the  sur¬ 
face  orientation  fulfilling  these  requirements  and  character¬ 
izing  this  orientation  by  the  surface  slopes  in  two  chosen 
orthogonal  planes  (which  depend  on  the  local  surface  normal) . 
In  practice  these  are  chosen  in  the  upwind  and  crosswind 
directions.  Evaluation  of  the  probability  of  such  a  surface 
orientation  is  based  on  available  empirically  derived  wave 


slope  statistics  defined  in  terms  of  surface  slope  components 
chosen  (arbitrarily)  in  crosswind  and  upwind  directions.  The 
geometry  of  the  wind  driven,  rough  ocean  surface  is  illustrated 
in  Figure  7-1.  The  parameters  defining  the  sun/pixel/sensor 
geometry  are  listed  in  Table  7-1  with  accompanying  explana¬ 
tion.  Supplementary  discussions  of  the  rough  ocean  reflection 
geometry  may  be  found  in  Takashima  and  Takayama  (1981)  and 
Ahmad  and  Fraser  (1982) . 

A  surface  capable  of  producing  a  specular  return  to  the 
sensor  will  have  a  surface  normal  (in  the  rotated  or  primed 
coordinate  system) : 

n'  =  (cose/2,  sine/2,  0)  (7-10) 


where  the  angle  e  is  that  between  sun  and  sensor  position 
vectors  or: 


-1 ,  ~  ~ i 

e  =  cos  (o-s) 


(7-11) 


The  coordinates  of  n'  in  the  fixed  (unprimed)  (x,y,z)  system 
(required  to  evaluate  the  surface  slope  components)  are 
obtained  by  transforming  from  primed  to  unprimed  systems 
employing  the  transpose  of  transformation  matrix  M,  i.e.: 


n  =  M  n' 


(7-12) 


The  elements  of  the  matrix,  M,  are  given  by  the  components 
of  the  defined  coordinate  systems  or: 


(7-13) 
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s  Solar  unit  position  vector  with  respect  to  pixel 

o  Sensor  unit  position  vector  with  respect  to  pixel 


(x,y,z)  Cartesian  coordinate  system  with  respect  to  calm 
ocean  surface  with  origin  at  pixel,  x  oriented 
along  direction  of  prevailing  wind,  and  z  in 
direction  of  local  zenith 

yQ  Solar  zenith  angle  cosine  with  respect  to  2 

\i  Sensor  zenith  angle  cosine  with  respect  to  2 


WD  Azimuth  of  wind  direction  measured  clockwise  from 

north 


a  Azimuth  of  solar  position  vector  measured  clockwise 

from  north 

4>  Azimuth  angle  difference  between  solar  and  sensor 

position  vectors 

(x' ,y* ,z ') Rotated  Cartesian  coordinate  system  on  sloping 

ocean  surface  with  origin  at  pixel,  x'  in  direction 
of  s,  and  z  in  direction  of  (s  x  6) 

e  Angle  in  (s  x  6)  plane  between  s  and  6  (e/2  is 

the  angle  if  incidence) 

n  Unit  normal  to  surface  in  (x,y,z)  system 

n'  Unit  normal  to  surface  in  (x',y',z')  system 

A6  Azimuth  angle  difference  from  specular  direction 


Table  7-1.  Sun/Pixel/Sensor  Geometric  Factors  Defined 
in  Figure  7-1 
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<• 


(• 


where:  e ^ ( j  =  1,2,3)  are  vectors  x,y,z 


and 


e^(i  =  1,2,3)  are  vectors  x,,y',z'. 


The  rotated  (primed)  coordinate  system  is  simply  defined  in 
terms  of  the  unit  position  vectors  s  and  o,  i.e.: 


~  ~  ^  ^  .  A  A 

x  =  s,  z  =  s  x  o,  y  =  z  x  x 


where  from  spherical  geometry: 

•  a 

1 

(1-y2)2  cos (180  +  AS) 


(7-14) 


s  = 


(1-U^)2  sin(180  +  AS) 


(7-15) 


and: 


o  = 


(1-U2)  2  cos  (180  +  AS  +  A<j>) 
1 

(1-U2)2  sin(180  +  AS  +  A<j>) 


(7-16) 


Considering  a  sea  surface  height  distribution  £(x,y) 
given  by: 

z  =  C (x,y)  (7-17) 

or 

f  (x,y,z)  =  z  -€  (x,y)  =  0, 

the  normal  to  the  surface  will  be  defined  by  the  gradient 
function : 


-  _  -<H  -6C 

N  -  (  -  «7  '  11 


(7-18) 
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k. . 


The  required  surface  slopes  in  upwind  (i.e.  x)  and  cross- 
wind  directions  then  are: 


=  -«£  z  =  -§£ 

u  6x  '  zc  6y 


(7-19) 


Based  on  the  above  definition  of  the  surface  normal: 


2u  =  “ni/n3  '*  zc  =  “n2/n3 


(7-20) 


where  n^  are  the  components  of  the  unit  surface  normal  with 
respect  to  the  fixed  coordinate  system. 

Using  Equations  (7-10)  through  (7-20)  the  ocean  surface 
slope  components  (zu»zc)  and  the  angle  of  incidence  (i.e. 
e/2)  can  be  calculated  given  the  sun  and  sensor  position 
vectors  (s,o). 


7.4  Effective  Fresnel  Reflectivity  Formulation 

As  a  consequence  of  the  relationships  derived  in  the 
previous  section,  sea  surface  slope  components  in  the  upwind 
and  crosswind  directions  (zu  and  zc,  respectively)  required 
to  produce  a  specular  glint  can  be  evaluated  based  on  satel¬ 
lite  and  sun  ephemeris  data.  Assuming  such  a  surface  orien¬ 
tation,  the  magnitude  of  the  ratio  of  reflected  radiance  to 
the  incident  solar  radiance  with  incidence  angle,  c/2,  is 


(Towne,  1967;  Kraus,  1972): 


P (c/2)  =  \ 


sin2 ( e/2-<)  tan2 (e/2-<) 

o  "**  o  '  1  ~ 2.  L) 

sin  (c/2+k)  J  tan  (e/2+<)J 


where 


and  the  indices  of  refraction  for  air  and  sea  water  are: 

n  =  1.000  (7-23) 

cl 

n  =  1.338 
s 

The  first  term  in  (7-21)  is  the  Fresnel  reflectivity  coeffi¬ 
cient  for  the  vertically  polarized  radiance,  Iv,  while  the 
second  is  that  for  the  horizontally  polarized  radiance,  1^. 
The  polarization  (P)  is  defined  as  their  difference  divided 
by  their  sum.  The  mean  reflectivity  is  identical  to  Equa¬ 
tion  (7-21)  above.  These  values  are  illustrated  for  an 
air/sea  water  interface  as  a  function  of  incidence  angle 
(0^  =  e/2)  in  Figures  7-2  (a,b,c,d).  The  observed  increase 
in  mean  reflectivity  for  low  solar  elevation  angles,  for 
example,  is  easily  explained  in  terms  of  simple  arguments 
based  on  Fresnel's  relations  (cf.  Figure  7-2d) . 

Given  that  the  surface  orientation  is  characterized 
by  slopes  (zu'zc)  and  using  Equation  (7-21)  to  evaluate 
the  Fresnel  reflectivity,  a  relationship  such  as  (7-8) 
may  be  used  to  calculate  the  reflected  radiation.  However, 
as  discussed  above,  the  probability  of  such  a  surface  orien¬ 
tation  is  dependent  on  surface  environmental  conditions,  such 
as  wind  speed  and  thus,  the  frequency  of  occurrence  of  such 
an  orientation  given  the  wind  conditions  is  sought.  By  co¬ 
ordinating  aerial  photographs  with  surface  wind  measurements 
Cox  and  Munk  (1954)  developed  statistics  about  surface  re¬ 
flection  derived  from  observed  sunglint  and  hence  inferred 
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Figure  7-2.  Calm  ocean  surface  reflectivities  based  on  Fresnel  formulae  for  air/sea  water  interface: 

(a)  vertical^  polarization,  (b)  horizontal  polarization,  (c)  degree  of  polarization,  and 


the  probability  density  distribution  function  (PDF)  for  sur 

face  slopes  (z  , z  )  with  wind  speed  w(ms-^) : 
u  c 


P<Vzc>  * 


2ira  cr 
u  c 


expf-  1  j(^)2  +  (p) 
Leu 


(7-24) 


where  the  variances  in  the  crosswind  and  upwind  directions, 
respectively,  are  (Guinn  et  al.,  1979): 


a  2  =  .003  +  . 00192w 
c 

a  2  =  .00316w 
u 


(7-25) 


Figure  7-3  (a,b) ,  respectively,  illustrate  the  variation  of 
the  crosswind  and  upwind  probability  density  distribution 
functions  for  various  wind  speeds  assuming  zero  slope  in  the 
orthogonal  direction.  As  might  be  expected,  there  is  a  higher 
probability  of  non  zero  surface  slopes  with  increasing  wind 
speeds  and,  conversely,  a  higher  probability  of  near  zero 
surface  slopes  for  lower  wind  speeds.  When  the  slope  in  the 
orthogonal  direction  is  already  non  zero  (cf.  Figure  7-4 
(a,b)  when  it  is  0.5),  the  PDF  is  predisposed  to  higher 
wind  speeds  and  the  behavior  near  zero  slopes  illustrated 
in  Figure  7-3  is  not  observed. 

It  is  important  to  note  that  the  surface  slope  orienta¬ 
tion  statistics  given  by  Equations  (7-24,  7-25)  represent  a 
time  averaged  picture  of  the  actual  instantaneously  changing 
surface  slope  orientation  and  are  somewhat  simplified.  For 
example,  it  can  be  shown  that  the  true  sea  height  and  its 
derivatives  cannot  be  truly  Gaussian  (for  a  discussion,  see 


[probability]  Log  [probability] 


M '  4  r>  k 


Crosswind  slope,  z 


(a)  Crosswind  probability  density  distribution  (z 


-  1 

-1.3 


Cox- M  u  n k 


iC  =  u 


i  \  *.>;v 
i  \  \\'v 

I  •  \  *  V 


■■  ///  i  i  \  \\\\ 

"////  i 1  M*\\\\ 

////  I  rwx 

I  /  /  /  \  ^S\2S 


I  \  \  \ 

1  \  %S\25 


Upwind  slope,  zu 
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Figure  7-3.  Sea  surface  slope  probability  density  function  (PDF) 
according  to  Cox  and  Munk  (1954)  in:  (a)  crosswind 
and  (b)  upwind  directions. 
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'•  Figure  7-4.  Sea  surface  slope  probability  density  function  (PDF)  according 

to  Cox  and  Munk  (1954)  in:  (a)  crosswind  and  (b)  upwind 
directions. 


7-17 


Barrick,  1972) .  This  may  be  remedied  by  adopting  a  product 
function  Gr(a,b)  to  modify  the  distribution  in  (7-24)  above 
where  Gr(a,b)  is  the  sum  of  products  of  Hermite  polynomials, 
the  Gram-Charlier  sum  (Guinn  et  al.,  1979).  The  difference 
is  very  slight,  however,  and  for  most  analytical  purposes, 
the  Gaussian  is  entirely  adequate  (see  MacKay,  1959) . 

Given  the  Fresnel  reflectivity  (7-21)  and  the  probabil¬ 
ity  distribution  function  (7-24),  an  effective  Fresnel  re¬ 
flectivity,  Rf,  is  defined  proportional  to  the  product  of  the 
Fresnel  reflectivity  for  a  given  surface  orientation  and  the 
probability  of  occurrence  of  that  orientation  [as  character¬ 
ized  by  its  surface  slope  components,  (z^,zc)]  dependent 
upon  the  surface  wind  speed  through  the  Cox-Munk  statistics, 
i ,e.  : 


Rf  p  (e/2)  p  (zc/Zu) 


(7-26) 


Thus  in  the  absence  of  an  atmosphere,  the  reflected  radiance, 
Iref,  above  a  rough  ocean  surface  due  to  incident  irradiance 
ttF  is  (Mullamaa,  1964)  : 


I 


ref 


ttFp  ( e/2)  p  ( z  ,  z  ) 

U 

4y  (ft  •  2)  ^ 


(7-27) 


The  effective  Fresnel  reflectivity  defined  above  is  used 
in  this  work  to  evaluate  surface  reflected  contributions 
to  the  radiance  in  potential  sunglint  regions.  The  behavior 
of  this  function  is  illustrated  in  Figures  7-5  and  7-6  for 
a  hypothetical  sun  at  a  zenith  angle  at  45°.  For  strict 
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specular  reflection  (Equation  (7-7)),  the  predicted  mean 
reflectivity  (Equation  (7-21))  at  the  mirror  angle  (45°) 
is  0.0285.  Figure  7-5  shows  the  wind  speed  dependent  Gaus¬ 
sian  spreading  of  the  reflectivity  due  to  the  Cox  and  Munk 
PDF  about  the  calm  sea  primary  specular  point  (PSP)  at  45° 
zenith  angle.  As  wind  speed  increases  and  the  sea  roughens 
(for  example  at  15ms  ^) ,  the  angular  dependence  of  reflection 
becomes  less  specular  and  more  Lambertian.  Figure  7-6 
contrasts  the  dependence  of  the  effective  reflectivity  on 
wind  speed  at  the  predicted  PSP  and  away  from  the  PSP  but 
in  a  region  of  potential  sunglint  (here  chosen  at  6  =  22.5°). 
Away  from  the  PSP  reflectivity  actually  increases  with  wind 
speed,  while  at  the  PSP  rougher  seas  decrease  predicted 
return  from  the  surface.  The  behavior  of  the  effective 
reflectivity  defined  above  is  qualitatively  consistent  with 
that  based  on  examination  of  actual  imagery  (cf.  Fett  and 
Mitchell,  1977) .  In  order  to  provide  an  example  of  the 
variation  of  surface  reflectivity  along  an  actual  DMSP 
visible  scan  line,  a  case  used  by  Tsui  and  Fett  (1980)  in 
their  discussion  of  the  primary  specular  point  is  adopted 
for  analysis.  The  required  input  parameters  characterizing 
the  example  are  given  in  Table  7-2.  These  data  are  used 
in  calculating  solar  and  sensor  position  vectors  for  each 
of  eleven  surface  pixels  using  the  methodology  described 
in  Section  5.  The  specific  scan  line  chosen  on  this 

revolution  is  that  located  at  185827  GMT  which  appears  to 
contain  the  approximate  position  of  the  primary  specular 


point  at  (10. 7N,  113. 9W).  Results  are  shown  in  Figure 
7-7  (a)  and  (b  )  illustrating  the  dependence  on  wind  speed 
and  wind  direction,  respectively. 


Table  7-2 

Input  Background  Parameters  for 
Surface  Reflectivity  Calculation 


Solar  Declination  Angle  (22  Aug  78) 
Longitude  of  Ascending  Node 
Time  of  Ascending  Node 
Satellite  Orbital  Inclination 
Satellite  Altitude 
Earth  Radius 


12.0° 

113. 5°W 
18:55:31  GMT 
98.7° 

833  km 
6370  km 


Desired  Scan  Line 


18:58:27  GMT 


Figures  7-8  and  7-9  illustrate  the  application  of  (7-27) 
to  the  evaluation  of  reflected  radiance  above  the  rough  ocean 
surface  in  the  absence  of  an  atmosphere  (cf .  Mullamaa,  1964) . 
Assuming  tt  units  of  incident  flux  and  a  solar  zenith  angle 
of  58.7°,  Figure  7-8  shows  the  variation  of  the  shape  of  the 
sunglint  region  in  the  principal  plane  (A<p  =  180°)  with 
zenith  angle  (or  equivalently  sensor  nadir  angle) .  As  wind 
speed  increases  the  glint  region  spreads  over  a  much  broader 
area  as  the  radiances  near  the  PSP  (say  60°)  decrease  and 
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Figure  7-7. 


Variation  of  effective  reflectivity  along  simu¬ 
lated  scan  line  at  185827  GMT  for  assumed  satel 
lite  background  parameters  (LN=113.5W,  TN=18:55 
31,  SD=12.0°)  illustrating  dependence  on:  (a) 
wind  speed  and  (b)  wind  direction. 
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of  tt)  with  zenith  angle  c.nd  wind  of  tt )  with  wind  speed  an 

speed  in  the  principal  plane  (no  angle  about  the  primary 

atmosphere).  point  (6=58.7°,  A<f>=180°) 


those  away  from  the  PSP  (say  30°)  increase.  Similar  varia¬ 
tions  are  notable  in  azimuth  angle  from  Figure  7-9  where  the 
spread  of  surface  reflected  energy  from  the  vicinity  of  the 
PSP  ( A4>  =  180°,  0  =  60°)  is  illustrated  for  the  wind  speed 
transition  from  2  ms  ^  to  10  ms  \ 

7.5  Incorporation  Within  the  Radiance  Simulation 

In  the  presence  of  an  atmosphere,  incorporation  of  the 
rough  ocean  surface  reflectance  model  within  the  radiance 
simulation  requires  some  formal  coupling  with  the  radiative 
transfer  treatment.  Unfortunately,  the  stream  approximation 
approach  to  the  solution  of  the  radiative  transfer  equation 
described  in  Section  6  is  not  ideally  suited  to  the  treatment 
of  specular  reflection.  This  is  a  consequence  of  the  boundary 
condition  on  fluxes  at  the  surface  given  by  Equation  (7-6) . 
Since  reflected  intensities  for  a  Lambert  (i.e.  non-specular) 
surface  are  independent  of  direction,  the  flux  boundary  con¬ 
dition  is  sufficient.  For  the  specular  case,  however,  re¬ 
flected  intensities  are  a  strong  function  of  direction  (7-8, 
7-9)  and  a  rigorous  boundary  condition  on  the  intensities 
should  be  used  in  the  radiative  transfer  solution  (Raschke, 
1972;  Plass  et  al.,  1975;  Guinn  et  al.,  1979;  Ahmad  and 
Fraser,  1981) .  For  example,  for  a  surface  with  general  bi¬ 
directional  reflection  distribution  function,  f^,  the  upward 
radiance  at  the  surface  I+(x*,u,ci>)  is  related  to  the  down¬ 
ward  radiance  l”(x*  ,u ,  <J>)  in  the  direction  (u,<J>)  by: 


Within  the  context  of  the  stream  approximation  this  is  not 
possible.  One  suggested  approach  is  to  treat  the  specular 
reflection  of  the  direct  solar  beam  and  Lambert  reflectance 
of  the  multiply  scattered  field  independently.  Thus,  for 
specular  cases,  an  appropriate  radiative  transfer  boundary 
condition  may  be  formulated  in  analogy  to  (7-6)  above  as 
(Kerschgens  et  al.,  1978): 

F+(t  =  x*)  -  A-JirFy  le“T  /Mo  +  ATF*(x  =  t*)  (7-27) 

where  AL  is  the  wavelength  dependent  Lambert  surface  reflec¬ 
tance  and  AF  is  the  specular  reflectance  of  the  direct  solar 
beam  determined  by  Fresnel's  formulae. 

The  approach  adopted  in  the  current  radiance  simulation 
model  is  based  on  further  simplification  of  an  approximate 
method  for  treating  non-Lambertian  surfaces  suggested  by 
Kaufman  and  Joseph  (1982) .  The  satellite  incident  radiance 
in  sunglint  regions,  Ir,  is  approximated  as  the  sum  of  two 
contributions:  (a)  radiance  due  to  the  specular  reflection 

of  the  direct  solar  beam  by  the  wind-roughened  ocean  surface, 
Ig,  and  (b)  radiance  due  to  back  scattering  by  the  atmosphere, 

I.  ,  i .e .  : 

A 
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[r  =  W  +  W 


(7-28) 


The  former  term  is  evaluated  based  on  the  effective  Fresnel 
reflectivity  defined  earlier,  while  the  latter  term  depends 
on  the  usual  Lambert  surface  albedo.  In  the  immediate  vicin¬ 
ity  of  the  PSP  for  near  calm  conditions  the  contribution 
from  (a)  above  clearly  dominates  that  from  (b) .  In  regions 

remote  from  sunglint,  the  converse  condition  holds.  In  an 

* 

atmosphere  of  total  optical  depth  t  ,  the  surface  reflected 
direct  solar  term  is: 


Is(Rf)  =  ttF  Rf  exp (-t  /uQ)  exp(-t  /u) 


(7-29) 


where  irF  is  the  wavelength  dependent  incident  solar  radiation 
Rf  is  the  effective  Fresnel  reflectivity  (7-26)  and  the  expo¬ 
nentials  are  attenuation  factors  for  air  masses  of  yQ  ^ ,  y  1 
in  the  solar  and  sensor  directions,  respectively.  The  atmo¬ 
spheric  term  is  evaluated  as  described  in  Section  6. 

This  approximation  essentially  decouples  the  surface 
reflected  radiance  from  that  due  to  atmospheric  scattering, 
an  approach  often  used  in  the  formulation  of  remote  sensing 
problems  over  land  (cf .  Fraser  and  Curran,  1976) .  It  also 
ignores  the  relation  of  diffuse  and  specularly  reflected 
components  available  from  (7-1),  i.e.: 

2  7T  1 

A.  =  r  f  f  R*  Udyd<{>  (7-30) 
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However,  as  a  practical  matter,  comparison  of  a  variety  of 
results  based  on  (7-28)  with  accurate  numerical  treatments 
(such  as  Plass  et  al.,  1975;  Ahmad  and  Fraser,  1982)  indicate 
that  this  approach  provides  a  quite  reasonable  simulation  of 
the  qualitative  behavior  of  the  sunglint  phenomenon.  For 
example.  Figure  7-10  is  a  comparison  of  results  using  (7-28) 
to  a  numerical  solution  (Ahmad  and  Fraser,  1982)  for  wind 
speeds  of  5.15  and  10.30  ms  ^  and  two  solar  zenith  angles: 

(a)  20°  and  (b)  57°.  Radiances  are  shown  in  the  principal 
plane  for  a  variety  of  sensor  zenith  angles  extending  from 
directions  away  from  sunglint  (<{>  =  18  0)  to  those  in  the 
vicinity  of  the  primary  specular  points  (fj>  =  0;  9  =  20°, 

57°,  respectively).  Note  that  the  shape  of  the  sunglint 
region  ( 4>  =  0 )  and  its  variation  with  wind  speed,  solar 
position,  and  sensor  zenith  angle  are  simulated  quite  well 
by  the  approximation.  Magnitudes  of  the  radiances  away  from 
the  respective  primary  specular  points  and  particularly  in 
the  anti-solar  direction  ($  =  180°)  underestimate  the  accurate 
calculations  due  to  the  particular  treatment  of  the  atmospheric 
scattered  contribution  used  in  the  approximation  results 
illustrated  (a  zero  Lambert  surface  albedo  was  used) . 

7.6  Radiance  Simulations  in  Sunglint  Regions 

Incorporating  the  surface  reflection  submodel  within 
the  overall  simulation  algorithm  provides  the  capability  to 
model  radiances  in  sunglint  regions.  As  an  example,  the 
sensitivity  of  simulated  radiance  to  variations  in  relative 
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Figure  7-10.  Comparison  of  Radiance  Results  from  Numerical  Simulation 
of  Ahmad  and  Fraser  (1982)  to  Approximate  Results  Summing 
Surface  Specular  Reflected  and  Atmospheric  Scattered  Con¬ 
tributions  Only.  [Exact:  5.15  ms"^-  ( - ),  10.3  ms--'- 

( - );  Approximate:  5.15  ms--'-  (©>  ,  10.3  ms--'-  (•)  ]  . 

Solar  Zenith  Angles  of  (a)  20°  and  (b)  57°. 
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humidity  and  wind  speed  are  illustrated  in  Figures  7-11  and 
7-12,  respectively.  These  calculations  are  based  on  the  same 
input  data  and  hence  scan  line  surface  reflectivity  calcula¬ 
tion  provided  in  Figure  7-7.  Figure  7-11  illustrates  the 
predicted  effect  of  an  increase  in  relative  humidity  from 
50%  to  95%  in  the  vicinity  of  the  primary  specular  point+ 

(pixel  number  2  in  Figure  7-7)  assuming  a  wind  speed  of  5.0 
ms  1  and  a  maritime  aerosol  model  with  a  0.4  km  scale  height. 
Radiances  decrease  by  about  15%  as  relative  humidity  increases 
within  the  indicated  range.  Previous  simulations  for  non- 
sunglint  conditions  (Isaacs,  1980)  predicted  increased  radi¬ 
ances  with  increased  relative  humidity  (see  Figure  2-1)  pri¬ 
marily  due  to  enhanced  aerosol  backscattering .  While  scattered 
contributions  are  also  enhanced  in  the  present  situation,  the 
gross  effect  is  determined  by  a  decrease  in  the  magnitude  of 
the  surface  reflected  solar  radiation  due  to  extinction  of 
the  solar  beam  both  prior  to  and  subsequent  to  surface  re¬ 
flection. 

The  sensitivity  of  simulated  radiances  to  variations 
in  wind  speed  along  the  chosen  scan  line  is  illustrated  in 
Figure  7-12.  The  base  case  is  a  relative  humidity  of  50% 
and  a  maritime  aerosol  model  with  scale  height  of  1.2  km. 

Pixel  2  is  closest  to  the  PSP  while  pixels  1  and  3  are  on 
either  side  (pixel  1  is  on  the  subsatellite  track) .  At  low 

^The  location  calculated  for  this  pixel  assuming  nominal 
satellite  orbital  parameters  is  10.5°N  latitude,  114. 5°W 
longitude . 
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Relative  Humidity 

Dependence  of  simulated  DMSP  radiance  (mW 
cm“2  ym“l  sr~l)  in  vicinity  of  primary  spe¬ 
cular  point  on  relative  humidity.  (Maritime 
aerosoi  model,  H=0.4  km,  ws=5.0  ms“l) 
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Figure  7-12 


Dependence  of  simulated  DMSP  radiance  on 
wind  speed  in  the  vicinity  of  the  primary 
specular  point  (Pixel  #2)  on  scan  line  at 
185830  GMT. 
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wind  speeds,  the  radiance  returned  to  the  sensor  is  quite 
specular.  As  wind  speed  increases,  however,  the  pixels 
along  the  scan  line  become  almost  indistinguishable  due  to 
the  surface  roughness.  Note  that  the  radiances  from  pixels 
neighboring  the  PSP  actually  increase  somewhat  as  their 
probability  of  returning  sunglint  to  the  sensor  and,  hence, 
effective  reflectivities  increase  with  increasing  wind  speed 
In  each  of  the  cases  described  above,  the  aerosol  scale 
height  has  remained  constant  while  either  wind  speed  or  rel¬ 
ative  humidity  have  been  varied.  In  the  example  to  follow 
the  modifying  effect  of  increased  aerosol  loading  in  a  sun- 
glint  region  is  investigated.^  The  case  examined  is  for  a 
maritime  aerosol  with  relative  humidity  of  50%.  The  scan 
line  chosen  is  at  15:06:10  GMT  on  7/27/75  for  an  ascending 
pass  with  ascending  node  of  94.0°  West  longitude  at  15:00: 

08  GMT.  The  solar  declination  angle  for  this  date  is  19.2°. 
Figures  7-13  and  7-14  illustrate  the  effect  of  increasing 
aerosol  scale  height  on  the  scan  line  sunglint  pattern  assum 
ing  surface  wind  speeds  of  1.5  and  5.0  ms  ^ ,  respectively. 

As  the  scale  height  increases  the  dominant  effect  is  a  radi¬ 
ance  decrease  near  the  primary  specular  point  (located  near 
pixel  number  6) .  There  is  an  accompanying  slight  increase 
for  pixels  remote  from  the  glint  region  (numbers  1-2,  10-11) 
Thus,  the  aerosol  layer  acts  to  smooth  out  the  radiance 
gradient  due  to  surface  reflection. 

This  case  is  the  subject  of  Section  7.7. 
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Figure  7-13.  Dependence  of  radiance  on  aerosol  scale  height 
in  sunglint  area.  (Maritime  aerosol,  RH  =  50%, 

WS  =  1.5  ms-1,  AN  =  94°,  TN  =  15008,  GMT  =  150610 
SD  =  19.2° .) 
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Figure  7-14.  Dependence  of  radiance  on  aerosol  scale  height 
in  sunglint  area.  (Maritime  aerosol,  RH  =  50%, 

WS  =  5’ms-1,  AN  =  94°,  TN  =  150008,  GMT  =  150610, 
SD  =  19.2°.) 


7.7  Data  Intercomparison  Case  Study 

Hardcopy  imagery  illustrating  the  typical  sunglint  pat- 
tern  observed  xn  DMSP  visible  data  was  obtained  for  the 
purpose  of  conducting  a  validation  of  the  modeling  metho¬ 
dology  described  in  previous  sections.  The  data  utilized 
was  for  a  morning  ascending  revolution  of  FTV  7529  on  27 
July  1975  located  over  the  Gulf  of  Mexico  with  ascending 
node  at  94.0°  West  longitude  at  15:00:08  GMT.  The  imagery 
is  illustrated  in  Figure  7-15. 

The  simulation  model  was  exercised  for  a  scan  line 
chosen  at  150610  GMT  in  a  relatively  cloud  free  area  of  the 
Gulf  of  Mexico  in  the  vicinity  of  the  primary  specular  point 
(located  at  18.7°N,  92.3°W).  Concurrent  meteorological  data 
used  as  input  for  the  simulation  indicated  variable  surface 
winds  of  between  0  and  5.0  ms  ^ ,  relative  humidities  of 
90-95%,  and  surface  visual  ranges  of  90-100  km.  The  surface 
visual  range  and  relative  humidity  values  were  jointly  used 
in  conjunction  with  the  maritime  aerosol  model  to  derive  an 
appropriate  scale  height  value  (0.08  km).  Results  for  the 
surface  effective  refleOtivity  and  scan  line  radiance  vs. 
pixel  number  are  shown  in  Figures  7-16a  and  b,  respectively. 

In  order  to  accomplish  a  comparison  with  the  hardcopy 
imagery,  film  transmission  density  was  read  directly  from 
the  data  using  a  30ym  square  sampler  along  the  150610  GMT 
scan  line.  Readings  were  taken  every  5mm  for  eventual 

+R.  Fett,  NEPRF . 
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alignment  with  simulated  pixels.  Using  the  gray  shade  scale 
provided  on  the  imagery,  density  readings  were  expressed  in 
terms  of  one  of  the  64  available  gray  shades.  These  in  turn 
were  converted  back  to  digital  word  inputs  by  reversing  the 
contrast  enhancement  process  using  the  appropriate  (in  this 
case  MODl-Low)  mapping  curve  digitized  from  Fett  and  Mitchell 
(1977)  (see  Figure  7-17) .  In  this  manner,  the  film  density 
trace  was  converted  to  values  roughly  proportional  to  the 
original  radiance  values  used  to  construct  the  image  matrix. 
Since  these  values  are  relative,  it  is  not  possible  to  assign 
them  absolute  radiance  values.  Instead  it  was  decided  to 
calibrate  the  data  using  the  model  simulated  result  at  pixel 
number  6  for  scaling,  since  the  model  performs  best  close 
to  the  predicted  PSP.  The  result  of  this  comparison  is  shown 
in  Figure  7-18.  Considering  the  assumptions  adopted  in  the 
analysis  and  the  difficulty  in  treating  surface  features 
and  cloud  elements  readily  apparent  in  the  imagery,  the 
performance  of  the  simulation  is  reasonable. 

7.8  DMSP  Imagery  Simulation 

By  simulating  scan  lines  (as  in  Section  6.5)  in  an 
appropriate  successive  manner,  segments  of  DMSP  imagery  may 
be  simulated.  In  order  to  demonstrate  this  capability, 
six  sample  cases  (see  Table  7-3)  of  input  meteorological 
data  (consisting  of  wind  speed,  relative  humidity,  and 
scale  height)  were  used  to  simulate  portions  of  DMSP  imagery 
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Figure  7-17.  DMSP  enhancement  mode  mapping 

curves  (from  Fett  and  Mitchell, 
1977) . 
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Figure  7-18. 


Comparison  of  model  simulated  radiances 
(-)  and  data  (®)  obtained  from  sensometry 
of  corresponding  hard  copy  imagery  in  the 
vicinity  of  sunglint.  (FTV  7529,  7/27/75 
150610  GMT,  AN  =  94°,  TN  =  150008,  SD  = 
19.2°,  MODI -LOW) . 


corresponding  to  the  satellite  parameters  given  in  Table 
7-2.  The  maritime  aerosol  model  was  employed  in  these  cal¬ 
culations.  Results  for  these  imagery  simulation  cases  con- 

-2  -1  -1 

toured  in  radiance  units  of  mWcm  ym  sr  are  illustrated 
in  Figures  7-19  (a-f) . 


Table  7-3 

DMSP  Imagery  Radiance  Simulation  Cases 
Case  Wind  Speed  (WS)  Relative  Humidity  Scale  Height  Figure 


# 

ms-1 

(%) 

(km) 

# 

1 

0.0 

50 

1.2 

7-19a 

2 

5.0 

50 

1.2 

7-19b 

3 

10.0 

50 

1.2 

7-19c 

4 

5.0 

50 

0.4 

7-19d 

5 

5.0 

80 

0.4 

7-19e 

6 

5.0 

90 

0.4 

7-19f 

The  ordinate  corresponds  to  the  subsatellite  track  and  is 
divided  into  five  segments  hatch-marked  at  30  second  inter¬ 
vals  commencing  at  185730  GMT.  The  corresponding  abscissa 
represents  the  eastward  scan  portion  for  this  morning  satel¬ 
lite  and  is  divided  into  ten  equal  segments  (i.e.  eleven 
pixels)  using  (5-29)  (e.g.,  the  odd  pixels  in  Table  5-2). 

Each  figure  thus  simulates  2h  minutes  of  imagery  data  assum¬ 
ing  uniform  meteorology  within  the  field  of  view.  (This 
restriction  is  easily  removed.)  Cases  are  selected  to 
provide  a  representative  variation  of  near  surface  meteor¬ 
ology  with  a  calm  sea  base  case  (Table  7-3,  case  1),  cases 
1-3  illustrating  the  effect  of  a  wind  speed  increase,  cases 


4-6  that  of  an  increase  in  relative  humidity,  and  a  compari¬ 
son  of  cases  2  and  4  the  effect  of  a  decrease  in  aerosol 
scale  height. 

For  the  calm  sea  (case  1,  Figure  7-19a) ,  a  bright  spe¬ 
cular  point  (at  asterisk)  appears  at  the  approximate  position 
of  the  PSP.  Away  from  the  PSP,  radiance  levels  are  determined 
by  aerosol  backscatter  and  Lambert  surface  reflectance.  The 
gradual  increase  in  radiance  from  the  subsatellite  track  to 
the  eastern  edge  of  the  data  is  attributable  to  atmospheric 
scattering  effects.  As  wind  speed  increases  (Figures  7-19b, 
c)  and  the  sea  surface  roughens,  the  intensity  of  the  PSP 
decreases  drastically  as  radiance  levels  increase  for  the 
pixels  in  the  immediate  vicinity.  The  decrease  of  radiance 
at  the  PSP  and  the  distance  from  the  PSP  of  pixels  with 
increased  radiance  levels  are  proportional  to  wind  speed. 

At  a  wind  speed  of  10  ms  1  (Figure  7-19c) ,  the  radiance 
field  is  remarkably  uniform  differing  from  maximum  to  mini¬ 
mum  by  less  than  a  factor  of  2.  (This  factor  is  unity  for 
a  Lambert  surface . ) 

Cases  4,  5,  and  6  (Figures  7-19d,  e,  f,  respectively) 
illustrate  the  effect  of  an  increase  in  relative  humidity 
for  fixed  surface  wind  speed  and  aerosol  scale  height.  In 
general,  as  relative  humidity  increases,  radiances  in  the 
sunglint  region  decrease  due  to  increases  in  the  attenuation 
factors  of  the  dominant  surface  reflected  direct  term  (see 
7-29).  Away  from  the  sunglint  region,  increases  in  relative 
humidity  increase  simulated  radiance  due  to  enhanced  aerosol 


7-40 


1  2 

3  4  5  6  7  8 

9  10 

Figure  7-19e. 

Simulated  DMSP  radiances. 

case  5  : 

WS  =  5  ms-1,  RH  =  80%,  H 
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Figure  7-I9f.  Simulated  DMSP  radiances,  case  6: 

WS  =  5  ms-1,  RH  =  90%,  H  =  0.4  km. 
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backscatter.  Thus  increasing  wind  speed  and  relative  humidity 
have  similar  impacts  on  the  appearance  of  the  sunglint  region. 

Additional  ambiguity  is  introduced  by  the  effect  of  aero¬ 
sol  scale  height  variations.  Careful  comparison  of  cases  2 
(Figure  7-19b)  and  4  (Figure  7-19d)  illustrates  that  scale 
height  increases  likewise  decrease  the  glint  region  radiances 
while  increasing  those  in  regions  remote  from  the  glint  region. 
These  calculations  are  summarized  in  Figure  7-20  where  a  rep¬ 
resentative  scan  line  (t  =  185830  GMT)  of  radiances  is  plotted 
for  each  case. 
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Figure  7-20.  Radiance  Simulation  Scan  Lines  (see  text) 
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8.  DMSP  RADIANCE/OPTICAL  DEPTH  RELATIONSHIPS 


A  necessary  first  step  in  the  inference  of  surface  range 
from  satellite  data  is  the  retrieval  of  vertical  aerosol  opti¬ 
cal  depth  from  measured  satellite  incident  radiances  (Goroch, 
1982) .  The  meteorologically  dependent  radiance  simulation 
model  described  here  evaluates  aerosol  optical  depth  as  an 
intermediate  step  and,  therefore,  radiance/aerosol  optical 
depth  relationships  are  provided  as  an  adjunct  calculation. 
Additionally,  the  capability  exists  to  explore  both  the 
meteorological  and  aerosol  model  dependence  of  these  rela¬ 
tionships.  Finally,  radiance  results  may  be  bandpass  weighted 
for  the  DMSP  sensor  (or  other  sensors  in  the  0.4  to  l.lym 
spectral  region) . 

In  this  section  DMSP-bandpass  weighted  radiance/optical 
depth  relationships  are  evaluated  using  the  overall  simulation 
algorithm.  The  methodology  employed  is  based  on  first  cal¬ 
culating  wavelength  dependent  upward  radiance  spectra  for 
maritime,  rural,  and  urban  aerosol  models  for  relative  humi¬ 
dity  values  of  50,  70,  80,  90  and  95%.  For  each  input  rela¬ 
tive  humidity,  aerosol  optical  depth  is  parameterized  by 
aerosol  layer  depth  increases  in  increments  of  0.5  km  from 
0.0  to  2.0  km.  For  a  given  aerosol  layer  depth,  H^,  the 
aerosol  optical  depth  at  a  chosen  wavelength,  A  is: 

TA(A,Hi}  =  Hi8A  (A'RH)  (8-1) 


8-1 


where  8A  (X,RH)  is  the  aerosol  extinction  coefficient  corres¬ 
ponding  to  the  appropriate  choice  of  model  and  relative  humi¬ 
dity.  (Extinction  coefficients  are  obtained  from  the  tabula¬ 
tions  in  Shettle  and  Fenn,  1979)  . 

The  radiance  spectra  are  bandpass-weighted  according 
to  (6-24)  to  provide  DMSP  simulated  radiances.  Finally,  the 
bandpass-weighted  radiances  are  plotted  against  the  optical 
depths  given  by  (8-1) .  For  applications  oriented  purposes, 
linear  regressions  are  provided  based  on  the  derived  radiance/ 
optical  depth  relationships.  Calculations  are  performed  for 
a  nadir  viewing  sensor  and  solar  zenith  angle  of  50°.  The 
optical  depth  is  evaluated  at  0.55]jm  for  direct  comparison 
to  visual  range  as  required.  Optical  depths  at  other  wave¬ 
lengths  are  available  by  reference  to  the  aerosol  model  and 
employing  (8-1)  above. 

8 . 1  Radiance  Spectra 

Figures  8-1,  8-2,  and  8-3  illustrate  wavelength  depend¬ 
ent,  simulated  upward  radiance  spectra  for  relative  humidi¬ 
ties  of  (a)  50,  (b)  70,  (c)  80,  (d)  90,  and  (e)  95%,  for  the 
maritime,  rural,  and  urban  aerosol  models,  respectively. 

Each  plot  consists  of  five  curves  corresponding  to  aerosol 
layer  depths  of  0.0,  0.5,  1.0,  1.5,  and  2.0  km  (e.g.,  lowest 
radiance,  H  =  0.0;  highest  radiance,  H  =  2.0).  These  results 
provide  the  capability  to  evaluate  radiance/optical  depth 
relationships  at  discrete  wavelengths  or  for  narrow  bands 
within  the  DMSP  0.4-l.lum  spectral  interval  corresponding 
to  other  sensors. 
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Ficure  8-1. 


Wavelength-dependent  radiances  for  maritime 
aerosol  and  relative  humidities  of:  (a)  50%, 
(b)  70%,  (c)  80%,  (d)  90%,  and  (e)  95%.  Each 
plot  consists  of  five  curves  corresponding  to 
aerosol  layer  depths  of  0.0,  0.5,  1.0,  1.5,  a 
2.0  km,  respectively. 


igure  8-1.  (continued) 
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Figure  8-2  .  Favelcncth-dcpcndcnt  radiances  for  rural  aerosol 
and  relative  buTr>iditi^«*  of:  (a)  r'0S;,  (h)  70%, 
(c)  90%,  (d)  90%,  and  (o)  95%.  Farh  plot  con¬ 
sists  of  five  curves  corresponding  to  aerosol 
layer  depths  of  0.0,  0.5,  1.0,  1.5,  and  2.0  km, 
respectively. 
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Figure  8-3.  Wavelength-dependent  radiances  for  urban  aerosol 
and  relative  humidities  of:  (a)  r>0% ,  (b)  70%, 

(c)  80%,  (d)  90%,  and  (e)  95%.  Each  plot  con¬ 
sists  of  five  curves  corresponding  to  aerosol 
layer  depths  of  0.0,  0.5,  1.0,  1.5,  and  2.0  km, 
respectively. 
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Figure  8-3 


(continued) 


8.2  Bandpass  Weighting  and  Linear  Relationships 

For  each  aerosol  model,  relative  humidity,  and  aerosol 
layer  depth,  the  radiance  spectra  presented  in  Section  8.1 
are  DMSP  bandpass-weighted  using  (6-24)  and  aerosol  optical 
depths  are  calculated  at  0.55ym.  Tables  8-1,  8-2,  and  8-3 
give  these  values  for  the  maritime,  rural,  and  urban  models, 
respectively.  (The  total  optical  depth  t  includes  a  Rayleigh 

optical  depth  contribution  of  .096  at  0.55um.)  Radiances 
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are  in  engineering  units  of  (mW  cm  ym  sr  ) . 

Simulated  DMSP  radiance  vs.  t  values  tabulated  here  are 
plotted  for  each  relative  humidity  in  Figures  8-4,  8-5,  and 
8-6  for  the  three  aerosol  models.  For  applications  purposes, 
linear  regressions  are  provided  for  each  curve  in  Tables 
8-4,  8-5,  8-6,  for  maritime,  rural,  and  urban  models,  respec¬ 
tively  and  for  all  relative  humidities  taken  together. 

8.3  Discussion  of  Radiance/Optical  Depth  Relationships 

The  DMSP  radiance/optical  depth  relationships  illustrated 

in  Figures  8-4,  8-5,  and  8-6  are  remarkably  linear  regardless 

of  oarcsol  model  and  relative  humidity.  This  is  especially 

2 

apparent  from  the  r  values  for  individual  relative  humidity 
curves  given  in  Tables  8-4  through  8-6.  This  linear  behavior 
relating  aerosol  optical  depth  and  sensed  radiance  is  consis¬ 
tent  with  theoretical  calculations  and  experimental  verifi¬ 
cation  performed  for  a  variety  of  visible  sensors  including 
Landsat,  GOES,  and  NOAA-5  (Griggs,  1979;  Norton  et  al . ,  1980). 
h ’ure  8-8,  for  example,  illustrates  the  dependence  of  radiance 


Table  8-1.  Simulated  DMSP  radiances,  I(mWcm  ym  sr  ),  assuming  a  uniform  maritime 
aerosol  layer  of  depth  H{km)  and  indicated  relative  humidity  (RH) .  In 
each  column,  x  is  the  total  optical  depth  at  0.55  ym  wavelength. 


50%  70%  80%  90%  95% 
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Table  8-3.  Simulated  DMSP  radiances,  I (mWcm  ym  sr  ),  assuming  a  uniform  urban 
aerosol  layer  of  depth  H(km)  and  indicated  relative  humidity.  In  eac 
column,  x  is  the  total  optical  depth  at  0.55  ym  wavelength. 
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Figure  8-4.  Simulated  DMSP  radiance  (mWcm  ym  sr  )  vs. 

total  optical  depth  (at  0.55  ym  wavelength)  for 

indicated  relative  humidities  (maritime  aerosol) 


Radiance 


Radiance 


• 

50 

I  =  0.440  +  1 . 375t 

0.999 

70 

I  =  0.451  +  1.2281 

0.998 

• 

80 

I  =  0.440  +  1 . 220t 

0.995 

90 

T  =  0.424  +  1 . 259x 

0.994 

<• 

95 

I  =  0.402  +  1 . 309t 

0.992 

all 

I  =  0.441  +  1 . 24  8r 

0.992 

• 

Table  8-4. 

Straight  line  fits  to  radiance  simulations 
(maritime  aerosol) . 

• 

<• 

•# 

1# 
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 J 

h 


B 


RH 

Linear  Regression 

2 

r 

50 

I  =  0.444  +  1.300T 

0.998 

70 

I  =  0.444  +  1.274T 

0.998 

80 

I  =  0.441  +  1.237T 

0.997 

90 

I  =  0.429  +  1 . 244t 

0.996 

95 

I  =  0.417  +  1.265T 

0.995 

all 

I  =  0.445  +  1.232T 

0.995 

Table  8-5.  Straight  line  fits  to  radiance  simulations 
(rural  aerosol) . 
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Urban  Aerosol 


RH 

Linear  Regression 

r2 

50 

I  =  0.540  +  0 . 336t 

0.997 

70 

I  =  0.529  +  0 . 419t 

0.998 

80 

I  =  0.511  +  0 . 564t 

0.999 

90 

I  =  0.500  +  0 . 695x 

0.999 

95 

I  =  0.478  +  0 . 810t 

0.999 

all 

I  =  0.432  +  0 . 789t 

0.969 

Table  8-6.  Straight  line  fits  to  radiance  simulations 
(urban  aerosol) . 


<# 


t# 


(normalized  to  an  irradiance  of  unity)  on  optical  depth  for 
surface  reflectance  values  of  0.02  and  0.43  (Kaufman  and 
Joseph,  1982) .  This  calculation  was  done  using  the  radiative 
transfer  solution  given  in  Section  6  assuming  a  nadir  viewing 
sensor  with  solar  zenith  angle  of  60°  and  an  aerosol  model 
with  single  scatter  albedo  of  unity  (i.e.  non-absorbing)  and 
asymmetry  parameter  of  0.63.  The  aerosol  optical  properties 
represent  a  low  relative  humidity,  rural  aerosol  and  the 
simulation  is  appropriate  for  Landsat  2-MSS6.  Note  that 
for  low  surface  reflectance  (such  as  over  the  ocean)  the 
radiance/optical  depth  behavior  is  very  linear  over  a  wide 
domain  of  optical  depth  values. 

It  should  be  noted,  however,  that  there  are  discernible 
differences  among  these  linear  relationships  dependent  both 
on  aerosol  model  and  relative  humidity.  The  most  notable 
of  these  may  be  distinguished  by  examining  the  dependence 
of  the  slopes  of  the  I  vs.  x  relatior ships  for  each  aerosol 
model  on  relative  humidity.  For  the  maritime  model,  for 
example,  radiance  increase  with  optical  depth  is  more  rapid 
at  low  than  at  high  relative  humidities.  For  the  rural 
aerosol  model,  this  trend  is  noted,  however,  the  slopes 
(AI/At)  are  almost  independent  of  relative  humidity.  Finally, 
distinctive  behavior  is  noted  for  the  urban  model.  There  a 
separate  regression  line  is  required  to  characterize  each 
relative  humidity.  Furthermore,  the  slope  variation  with 
relative  humidity  is  opposite  to  that  noted  for  the  maritime 
and  rural  models,  i.e.  the  rate  of  charge  of  radiance  with 


i# 
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increases  at  higher  relative  humidities. 

Tnese  differences  may  be  explained  by  recalling  that  the 
radiative  transfer  theory  relating  aerosol  optical  properties 
to  simulated  radiance  requires  a  minimum  of  three  input  vari¬ 
ables  including  optical  depth,  single  scattering  albedo,  and 
scattering  phase  function  (cf .  equations  6-5  through  6-12) . 
These  optical  properties  are  dependent  on  aerosol  size  distri¬ 
bution  and  index  of  refraction  which  may  ultimately  be  related 
to  chemical  composition  and  meteorological  variables  such  as 
relative  humidity  (see  Isaacs,  1980;  Sections  4  and  5).  For 
a  given  chemical  composition,  therefore,  the  optical  proper¬ 
ties  cited  above  are  related  in  a  unique  manner  allowing 
optical  depth  to  act  as  a  surrogate  for  the  others.  The 
differences  observed  in  the  radiance/optica1  depth  relation¬ 
ships  above  are  thus  manifestations  of  the  subtle  differences 
due  to  differing  composition  dependent  relative  humidity 
growth  factors  and  aerosol  absorption  properties.  This  may 
be  illustrated  by  decoupling  the  three  optical  parameters 
t,  ioQ,  g  in  a  theoretical  calculation.  Figure  8-9  illustrates 
radiance  changes  for  a  fixed  optical  depth  of  unity  due  to 
variations  in  single  scattering  albedo  and  asymmetry  factor. 

For  a  fixed  albedo  of  to  =0.95  the  crosses  (  +  )  indicate  a 

o 

decrease  in  simulated  radiance  with  increases  in  asymmetry 
factor,  e.g.,  there  is  less  backscatter  to  the  sensor.  Con¬ 
versely,  for  a  fixed  g  factor,  the  circles  (o)  indicate  an 
increase  in  radiance  with  increasing  <o  (i.e.,  decreasing 
aerosol  absorption) .  Relating  these  results  to  the  radiance/ 
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optical  depth  relationships  illustrated  in  Figures  8-4, 

8-5,  8-6,  a  slope  (AI/Ax)  decrease  (increase)  may  be  asso¬ 
ciated  with  an  increase  (decrease)  in  asymmetry  factor,  g, 
or  decrease  (increase)  in  single  scattering  albedo,  ui  . 

As  relative  humidity  increases  aerosol  growth  is  mani¬ 
fested  in  an  effective  size  distribution  change  (Nilsson, 
1979) ,  increasing  the  number  of  larger  particles.  This  is 
illustrated  in  Figure  8-10  by  the  decrease  in  calculated 
Junge  coefficient  for  the  maritime  model  with  increasing 
relative  humidity.  A  corresponding  increase  in  g  factor 
occurs.  This  effect  is  largest  for  the  maritime  aerosol 
which  responds  most  drastically  to  relative  humidity  in¬ 
creases  (cf .  Figure  4-3b) .  Since  the  maritime  aerosol  model 
is  essentially  non-absorbing  (Figure  4-1) ,  its  single  scatter 
albedo  remains  fairly  constant  with  relative  humidity  and  the 
size  distribution  (i.e.,  g  factor)  increase  effect  dominates 
decreasing  the  AI/Ax  slope  at  higher  humidities.  The  behav¬ 
ior  noted  for  the  urban  aerosol  is  due  to  its  absorption 
component.  As  relative  humidity  increases  water  is  added 
to  the  aerosol  and  its  volume  mixed  index  of  refraction  be¬ 
comes  increasingly  less  absorbing.  The  resultant  single 
scattering  albedo  thus  increases  drastically  with  relative 
humidity  (cf.  Figure  4-7b) .  Although  aerosol  growth  also 
increases  the  g  factor  (Figure  4-7c) ,  the  albedo  increase 
effect  dominates  resulting  in  the  increasing  AI/Ax  slope 
observed  in  Figure  8-6.  The  rural  model  results  (Figure  8-5) 
are  an  intermediate  case  with  slightly  less  increase  in  g 
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Figure  8-9.  Radiance  changes  with  single  scattering 
albedo  and  asymmetry  parameter  for  fixed 
optical  depth  of  unity. 
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Figure  8-10.  Effective  size  distribution  change  with 
relative  humidity. 
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factor  and  greater  increase  in  u>  with  relative  humidity 
compared  to  the  maritime  model. 

A  unique  linear  predictor  equation  relating  radiance 
and  optical  depth  would  provide  a  significant  practical 
simplification  to  the  inference  of  surface  range  from  re¬ 
motely  sensed  radiance.  The  results  presented  here  suggest 
that  a  potential  ambiguity  is  introduced  by  the  presence  of 
absorbing  aerosol  species  such  as  the  carbonaceous  component 
of  the  urban  model.  Figure  8-11  illustrates  a  composite  of 
individual  I  vs.  t  relationships  for  all  relative  humidities 
and  each  aerosol  model.  A  single  linear  relationship  may 
suffice  for  both  maritime  and  rural  models.  Since  less 
solar  radiation  is  backscattered  to  the  satellite  sensor 
by  an  absorbing  -aerosol,  however,  a  given  remotely  sensed 
radiance  for  the  urban  model  corresponds  to  a  higher  optical 
depth  and  hence  higher  boundary  layer  extinction  (i.e. 
reduced  range)  than  predicted  assuming  single  scattering 
albedos  near  unity.  Since  the  presence  of  absorbing  compon¬ 
ents  is  suggested  even  in  remote  regions  (Rosen  et  al.,  1981; 
Gerber  and  Hindman,  1981)  and  the  usual  definition  of  visual 
range  must  be  modified  when  considering  absorbing  aerosols 
(Roessler  and  Faxrog,  1981) ,  inference  of  visual  range  may 
not  be  straightforward.  Additionally,  since  carbon  is  uni¬ 
formly  absorbing  throughout  the  near  IR  and  thermal  IR  spec¬ 
tral  regions  (Twitty  and  Weinman,  1971) ,  its  presence  has  a 
significant  impact  on  multispectral  range  prediction.  For 
these  reasons  a  method  of  inferring  both  optical  depth  and 
single  scattering  albedo  may  be  advisable. 


Total  Optical  Depth 
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riaure  8-11.  Simulated  DMSP  radiances  (mWcm  um  c r  )  vs 
total  optical  depth  (at  0.55  ym  wavelength) 
for  maritime  (o) ,  rural  (+) ,  and  urban  (*) 
aerosol  models  and  all  relative  humidities. 
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9.  MULTISPECTRAL  RANGE  RETRIEVAL 

9 . 1  Rationale 

The  performance  of  electro  optical  systems  in  the  lower 
atmosphere  is  constrained  by  wavelength-dependent  extinction 
due  to  a  variety  of  mechanisms  including  aerosol  scattering 
and  absorption/  molecular  (line)  absorption,  Rayleigh  scat¬ 
tering,  and  various  continuum  absorption  features.  Assuming 
that  the  extinction  (in  km  at  wavelength  due  to  the 
summed  contributions  of  the  above  phenomena  is  known  and 
given  by  8T(A^) ,  a  reliable  assessment  of  system  effective¬ 
ness  can  be  obtained  by  evaluating  the  inverse  extinction 
length  or  range  (km) : 

R(A±)  =  Ci/8T(Ai)  (9-1) 

For  example,  at  0.55ym,  equation  (9-1)  is  a  statement  of  the 
Koschmieder  (1924)  equation  with  =  3.9  and  R(0.55ym)  de¬ 
fining  the  visual  range.  One  immediate  application  of  the 
DMSP  radiance/aerosol  optical  depth  relationships  derived 
previously  (Section  8)  is  in  the  estimation  of  the  aero¬ 
sol  extinction  contribution  to  total  extinction  in  equa¬ 
tion  (9-1)  above.  Utilizing  this  capability,  a  prototype 
range  retrieval  algorithm  is  formulated  in  this  section 
to  infer  optical  range  in  the  near  surface  environment  from 
satellite  observations.  By  utilizing  available  models  of 
molecular  and  continuum  extinction  features ,  the  domain  of 
wavelength  applicability  may  be  extended  from  the  visible 
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region  to  include  the  4  and  10ym  windows. 


9.2  Aerosol  Contribution 

The  aerosol  optical  depth  x(Aj)  at  wavelength  A  ^  inferred 
from  the  previously  derived  DMSP  radiance/optical  depth  rela¬ 
tionships  is  fundamentally: 

OO 

tUj)  =  /  [^(A.fZ)  +  B2(Xj/Z)]dz  (9-2) 

o 

where  B-^  and  B2  are  contributions  to  total  extinction  due  to 
the  accumulation  and  coarse  modes,  respectively.  It  is 
assumed  that  (9-2)  may  be  approximated  as: 

r(Xj)  ~  61(Xj)  Hx  +  B2(Xj)  H2  (9-3) 

where  B1,  B2  are  modal  extinctions  representative  of  the 
near  surface  environment  and  H^,  H2  (km)  are  mixing  heights 
corresponding  to  each  mode.  This  supposition  requires  either 
a  uniform  extinction  layer  near  the  surface  or  a  vertical 
distribution  characterized  by  an  exponential  with  scale 
height,  .  Alternatively,  (9-3)  may  be  replaced  by  an 
appropriate  mixed  layer  model  of  the  vertical  extinction 
profile  to  relate  the  inferred  x(A^)  value  with  surface 
value  of  6^  and  82-  It  is  further  assumed  that:  (1)  the 
coarse  mode  contribution  is  calculable  from  the  Munn-Katz 
oceanic  wind  speed/relative  humidity  dependent  size  distri¬ 
bution  at  each  required  wavelength,  and  (2)  the  accumulation 
mode  size  distribution  is  fixed  by  relative  humidity  and 
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requires  only  normalization  to  a  number  density  (cm  ) 
using  the  input  optical  depth.  Thus,  given  input  of  satel 
lite  radiance  inferred  optical  depth,  t(X^),  relative  humi 
dity  (RH)  ,  windspeed  (WS) ,  and  the  mixing  heights  ,  H2, 
the  aerosol  extinction  contributions  at  any  set  of  X^  wave 
lengths  may  be  evaluated  based  on  the  following  algorithm 
steps : 

9.2.1  Calculate  extinction  coefficient  of  oceanic  mode 
at  X.,  65  (X.),  from  input  RH,  WS . 

9.2.2  Calculate  optical  depth  at  X ^  of  oceanic  model 
using  input  oceanic  scale  height,  H2 : 

W  *  H2B2(V 

9.2.3  Calculate  non-oceanic  (i.e.,  accumulation  mode) 

optical  depth  at  X ^  from  input  satellite- inferred 

optical  depth  t_(X.): 

S  J 

Tl(Xj}  =  Ts(Xj)  "  T2(Xj) 

9.2.4  Calculate  fine  mode  extinction  at  X.,  0-(X.)  using 

J  £  1 

input  A._.ie  mode  scale  height  and  equation  (9-3) 
above : 

6i(X.)  =  [ts(X.)  -  H2B2(X.)  lH^1 

9.2.5  Calculate  number  density  in  fine  mcde  by  ratio 
of  inferred  fine  mode  extinction  to  a  model  fine 


mode  extinction  per  unit  number  density  B^(Xj): 


S^X  ) 

N1  =  -M 
1  S“(  Xj) 


4# 
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9.2.6  Calculate  fine  mode  extinction  at  i  required  wave¬ 
lengths  using  and  model  extinctions  at  i  wave¬ 
lengths  for  input  relative  humidities: 

6i(Xi,RH)  =  U^RH) 

9.2.7  Calculate  ranges  at  each  required  X^  using  equation 
(9-1) . 

In  practice,  required  values  of  the  relative  humidity/ 
windspeed  dependent  oceanic  mode  extinction  are  obtained 
by  interpolation  between  a  suf f iciently-dense  set  of  stored 
extinction  values  using  a  modified  Kasten-type  (1969)  rela¬ 


tion  of  the  form: 


9.3  Non-Aerosol  Contribution 

In  addition  to  the  aerosol  contributions  3^,  B2,  th© 
total  extinction  may  include  additional  non-aerosol  contri¬ 
butions  due  to  Rayleigh  scattering,  BR,  molecular  line 
absorption,  Bq,  and  continuum  absorption  due  to  water  vapor, 
8^,  and  nitrogen,  ,  i.e.: 


0T(Xi)  =  B  <X±)  +  B2(Xi)  +  BR(Xi)  +  BG(Xi) 


+  B  ( X±)  +  BjjtXj.) 


(9-5) 


The  relative  importance  of  each  contribution  is  dependent 
on  wavelength. 


9.3.1  Rayleigh  scattering 

The  Rayleigh  scattering  coefficient  is  calculated  using 
the  L0WTRAN5  (Kneizys,  et  al.,  1980)  expression  suggested 
by  Shettle  et  al .  (1980)  based  on  a  least  square  fit  to  the 
values  of  Penndorf  (1957)  : 

Br(X)  =  [926.759A4  -  10.7123X2]  1  (9-6) 

Calculated  values  appropriate  for  sea  levels  are  shown  in 
Table  9-1. 

9.3.2  Molecular  line  absorption 

The  extinction  contribution  at  each  wavelength  due  to 
molecular  line  absorption  is  evaluated  by  summing  absorption 
coefficients  due  to  relevant  gaseous  absorbers  (from  among 
H20,  C02,  0^ ,  N20,  CO,  CH^,  and  02)  using  spectral  line 


parameters  obtained  from  the  most  recent  version  of  the 
AFGL  compilation  (McClatchey  et  al.,  1973).  For  each 
spectral  line  2.  and  absorbing  species  i,  the  absorption 
coefficient  0^  is  given  by: 

$^(v,P,T)  =  ni(P,T)  S£(T)fJl(v/P,T)  (9-7) 

.  -3 

where  is  the  species  l  number  density  (molecule  cm  )  at 

pressure  P  and  temperature  T.  is  the  line  strength  (cm-1/ 

_2 

molecule  cm  ) ,  and  f ^  is  the  line  shape  (cm) .  The  total 
gaseous  absorption  coefficient  is  obtained  by  summing  over 
all  lines  and  species: 

6g(v,P,T)  =£S8i(v'p'T)  (9-8) 

i  l 

The  specific  methodology  employed  is  described  in  Susskind 
and  Searl  (1977)  who  provide  details  of  the  model  atmosphere 
and  cut  off  criteria  for  the  inclusion  of  center  and  wing 
line  contributions.  Conceptually  the  calculation  is  similar 
to  that  done  in  the  LASER  computer  code  (McClatchey  and 
D'Agati,  1978) .  For  the  calculations  described  here,  the 
model  atmosphere  was  modified  to  input  relative  humidity 
to  determine  water  vapor  number  density  and  output  values 
were  weighted  with  a  square  instrument  filter  function  of 
half  width  0.5  cm  1.  Calculated  values  at  3.8ym  and  10.6um 
of  0.0276  km  1  and  0.0328  km  1,  respectively  [see  Table  9-1] 
correspond  primarily  to  CO2  and  1^0  contributions  (some  CH^ 
at  3.8ym) . 
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9.3.3  Continuum  absorption 

Water  vapor  continuum  absorption,  B„,  at  both  3.8  and 

n 

10.6ym  and  the  pressure  induced  nitrogen  continuum  absorp¬ 
tion,  Bn/  in  the  4ym  region  are  included  in  Table  9-1.  The 
nitrogen  continuum  is  based  on  the  measurements  of  Shapiro 
and  Gush  (1966)  as  described  in  Susskind  and  Searl  (1977) . 

In  the  8-llym  region  the  water  vapor  self  broadening  is 
evaluated  based  on  the  values  from  Bignell  (1970) .  Thes 
calculations  are  consistent  with  those  using  Burch  (1970) 
and  Roberts  et  al.  (1976)  as  applied  in  LOWTRAN5.  Near 
4ym,  the  water  vapor  continuum  is  based  on  the  values  given 
by  Burch  et  al.  (1971) .  Water  vapor  continuum  effects  are 
dependent  on  its  abundance  (partial  pressure,  mass  density, 
etc.)  and  thus  an  additional  input  parameter  is  required. 

As  described  previously  the  partial  pressure  of  water  vapor 
was  input  as  relative  humidity  for  a  fixed  layer  temperature 
(288.1  K  at  approximately  standard  pressure).  Thus  the 
additional  parameter  here  is  the  layer  temperature.  Figure 
9-1  illustrates  the  variation  of  the  total  water  vapor 
continuum  absorption  coefficient  B  at  10.6ym  with  relative 

ri 

humidity.  The  absolute  water  vapor  abundance  ranges  from 
2.0  gm  ^  (RH  ,2,15.7%)  to  10.9  gm  ^  (RH^.85%).  Also  illustrated 

in  Figure  9-1  are  the  individual  contributions  due  to  self, 

s  f 

8„,  and  foreign  broadening,  B„. 

ri  ri 
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9.4  Multispectral  Range  Retrieval  Results 


By  using  equation  (9-1) [where  the  total  extinction  is 
defined  by  (9-5) ]  and  evaluating  the  aerosol  contributions 
based  on  algorithm  steps  (9.2.1  -  9.2.7)  and  non-aerosol 
contributions  as  described  in  Section  9.3,  optical  ranges 
are  retrieved  for  the  wavelengths  and  conditions  listed  in 
Table  9-1.  These  results  are  presented  in  Table  9-1  under 
heading  (c) .  Note  that  while  range  is  almost  exclusively 
determined  by  aerosol  extinction  at  visible  and  near  IR 
wavelengths,  the  significance  of  non-aerosol  effects  in¬ 
creases  with  wavelength.  At  3.8ym  aerosol  extinction  and 
non-aerosol  effects  (primarily  line  absorption)  are  about 
equal  while  at  10.6ym,  optical  range  is  almost  solely  a 
function  of  water  vapor  abundance  due  to  continuum  absorp¬ 
tion. 

While  this  example  is  quite  case  specific,  it  illustrates 
a  methodology  to  retrieve  optical  ranges  at  desired  wave¬ 
lengths  using  satellite  based  visible  radiances  as  an  aid 
in  the  inference  of  aerosol  extinction  contributions.  For 
wavelengths  in  the  visible  region,  optical  range  may  poten¬ 
tially  be  inferred  from  the  satellite  radiances  alone  with 
some  supporting  meteorology  (RH,  WS ,  H) .  At  longer  wave¬ 
lengths  where  non-aerosol  contributions  to  extinction  be¬ 
come  significant  (such  as  the  water  vapor  continuum  near 
lOym) ,  additional  input  data  (such  as  absolute  humidity) 
is  required.  Such  information  may  also  be  potentially  avail¬ 
able  from  remote  sensing  data.  In  any  case,  the  technique 
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described  here  provides  the  framework  to  evaluate  a  priori, 
the  impact  of  various  meteorologically  dependent  extinction 
mechanisms  on  optical  range.  This  capability  may  be  exploited 
to  provide  a  multispectral  range  retrieval  algorithm  when 
supported  by  the  requisite  satellite  data  source.  Further 
discussion  of  data  needs  may  be  found  in  Goroch  (1982) . 


t# 
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10.  SUMMARY ,  CONCLUSIONS,  AND  RECOMMENDATIONS 

♦ 

10.1  Report  Summary 

This  report  has  described  the  results  of  the  second 
stage  of  model  development  and  application  undertaken  in 
the  in  restigation  of  the  sensitivity  of  remotely  sensed 
radiances  in  the  DMSP  VHR/LF  spectral  region  (0.4  to 
l.lym)  to  near  surface  meteorological  parameters  such  as 
relative  humidity,  wind  speed,  and  aerosol  characteristics 
including  mixing  height  and  number  density.  A  research 
grade  computer  code  developed  in  an  earlier  study  (Isaacs, 
1980)  to  approximately  simulate  DMSP  visible  and  near  infra¬ 
red  radiances  has  been  upgraded  and  augmented  with  additional 
capabilities  desired  in  the  analysis  of  imagery  for  tactical 
purposes  including:  (a)  treatment  of  continental  (i.e., 
urban  and  rural)  aerosol  size  distribution  and  composition 
models  in  addition  to  previously  implemented  maritime  models 
for  application  to  coastal  regions  or  open  ocean  areas  under 
the  influence  of  continental  origin  air  masses,  (b)  data 
simulation  calculations  based  on  realistic  sun/sensor  geo¬ 
metries  through  the  implementation  of  simplified  satellite 
subtrack/sensor  scan  line  and  solar  ephemeris  subroutines, 
and  (c)  approximate  ocean  surface  reflection  modeling  for 
potential  sunglint  regions  based  on  both  the  classical 
Fresnel  relations  and  the  empirical  wind  dependent  statis¬ 
tics  of  sea  surface  slope  orientations.  Additionally,  the 
approximate,  analytical  radiative  transfer  submodel 
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incorporated  into  the  code  is  fully  described  and  the 
results  of  comparisons  with  accurate  numerical  treatments 
are  provided  for  a  variety  of  scattering  situations.  Com¬ 
prehensive  program  documentation  is  included  consisting  of 
listings  of  the  algorithm,  sample  output,  and  operating 
instructions  for  implementation  on  a  desk  top  minicomputer. 

A  variety  of  simulation  studies  employing  the  algorithm 
are  discussed  illustrating  the  effect  of  relative  humidity 
and  mixing  height  on  radiances  in  non  sunglint  situations 
for  both  rural  and  urban  aerosol  models  as  compared  to  the 
maritime  models  previously  examined  (i.e.  in  Isaacs,  1980) . 
Radiance  calculations  in  sunglint  regions  using  the  maritime 
aerosol  model  illustrate  the  effect  of  relative  humidity, 
wind  speed,  and  aerosol  scale  height  on  the  intensity  gra¬ 
dient  observed  along  simulated  sensor  scan  lines  in  the 
vicinity  of  the  primary  specular  point.  Additionally,  the 
model  is  exercised  to  simulate  a  swath  of  DMSP  imagery  for 
illustrative  purposes.  Results  of  these  simulations  con¬ 
firm  the  behavior  generally  observed  in  relevant  imagery 
for  corresponding  meteorological  analysis  case  studies  as 
previously  reported  (cf .  Fett  and  Mitchell,  1973) . 

Comparison  of  model  predicted  values  to  measured  data 
are  presented  in  two  contexts:  (a)  multiwavelength  aerosol 
extinction  coefficient  profiles  calculated  from  field  mea¬ 
sured  size  distributions  are  compared  to  model  predicted 
values  based  on  measured  relative  humidity  profiles  as  the 
required  meteorological  input.  Predicted  and  measured 


values  are  equated  at  one  wavelength  thus  providing  a  measure 
of  relative  wavelength  dependence.  In  general,  the  maritime 
aerosol  model  performs  well  except  at  high  relative  humidi¬ 
ties.  (b)  Model  calculated  radiances  along  a  simulated 
scan  line  in  sunglint  were  compared  to  data  points  obtained 
from  sensometry  for  the  corresponding  hardcopy  imagery. 

The  model  input  data  included  the  appropriate  satellite  and 
solar  specific  data  (AN,  TN,  SD,  etc.)  and  relevant  meteor¬ 
ological  data  supplied  with  the  imagery.  Due  to  lack  of 
overall  data  calibration,  the  modeled  results  were  norma¬ 
lized  to  one  data  point.  The  resulting  comparison  is 
promising,  considering  the  assumptions  invoked. 

In  addition  to  these  imagery  simulation  exercises 
potentially  related  to  meteorological  analyses,  application 
of  the  algorithm  to  the  inference  of  surface  propagation 
parameters  was  explored.  Utilizing  the  model's  capability 
to  provide  simulated  radiances  for  specified  meteorology 
and  choice  of  aerosol,  DMSP  radiance/aerosol  optical  depth 
relationships  were  computed.  Results  provide  a  theoretical 
fit  of  radiance  to  aerosol  optical  depth  for  maritime, 
rural,  and  urban  aerosols  and  relative  humidities  of  50 
through  95%.  These  may  be  applied  to  the  retrieval  of 
aerosol  optical  depth  (at  visible  wavelengths)  from  appro¬ 
priately  calibrated  DMSP  radiances. 

Since  surface  propagation  parameters  are  desired  at 
other  than  visible  wavelengths,  a  prototype  algorithm  was 
devised  to  infer  optical  range  multispectrally  using  the 


DMSP  radiance/optical  depth  relationships  described  above 
to  provide  the  aerosol  contribution  to  extinction. 

10.2  Conclusions 

The  following  conclusions  and  observations  are  based 
on  the  work  summarized  above: 

•  Although  limited  in  scope,  the  validation  effort 
associated  with  the  MAGAT  data  set  indicates 
that  the  AFGL  maritime  aerosol  model  performs 
relatively  well  where  expected. 

•  A  variety  of  comparisons  performed  between  the 
model  adopted  approximate  radiative  transfer 
scheme  and  accurate  numerical  treatments  for 

a  range  of  optical  depths,  angular  scattering 
functions,  and  observer/sun  orientations  leads 
to  the  conclusion  that  simulated  radiances  are 
accurate  to  within  10%  overall.  This  level  of 
accuracy  is  consistent  with  the  accuracy  of  input 
data  and  the  intended  purpose  of  the  model,  i.e. 
the  simulation  of  meteorological  influences  in 
the  imagery.  The  results  of  specific  sensitivity 
tests  conclude  that  adoption  of  simple  analytical 
phase  function  representations  and  the  hybrid 
modified  6 -Eddington  flux  parameterization  are 
appropriate  choices. 

•  The  wind-driven,  rough  ocean  surface  reflection 
model  implemented  within  the  radiance  simulation 


code  provides  a  good  representation  of  the 
sunglint  phenomenon  as  compared  to  accurate 
numerical  treatments,  simulated  radiance 
comparisons  to  hardcopy  imagery,  and  quali¬ 
tative  observations. 

Model  simulations  in  sunglint  regions  suggest 
a  decrease  in  radiance  in  the  vicinity  of  the 
primary  specular  point  due  to  increases  in 
aerosol  optical  depth  (related  to  increased 
relative  humidity  or  scale  height) .  Away  from 
the  primary  sunglint  region,  radiances  increase 
with  increasing  optical  depth  as  predicted  in 
previous  studies . 

A  comparison  of  the  optical  properties  of 
the  maritime  and  recently  implemented  rural 
and  urban  models  suggests  that  some  care  must 
be  taken  in  the  interpretation  of  meteorological 
influences  on  individual  air  masses.  For  example 
the  maritime  aerosol  coarse  mode  is  much  more 
hygroscopic  than  the  others,  while  the  carbon 
component  of  the  urban  aerosol  is  critical  in 
determining  its  radiative  properties. 

As  a  consequence  of  the  above,  each  aerosol 
composition  model  and  relative  humidity  range 
exhibit  unique  radiance/optical  depth  relation¬ 


ships.  While  those  for  maritime  and  rural  aero- 


predictor  equation  to  be  used  operationally, 
those  for  the  urban  aerosol  introduce  a  potential 
ambiguity.  It  appears  that  it  may  be  important 
to  know  whether  the  air  mass  aerosol  under  investi¬ 
gation  has  any  carbon  component  present  (anthro¬ 
pogenic  or  natural) .  This  may  be  important  in 
the  vicinity  of  inland  seas,  coastal  regions,  or 
open  ocean  in  situations  of  long  range  transport. 

A  parameter  which  may  potentially  be  retreivable 
in  this  regard  is  the  single  scattering  albedo 
(Kaufman  and  Joseph,  1982) . 

•  A  sample  multispectral  range  retrieval  exercise 
illustrates  both  the  dominance  of  aerosol  extinc¬ 
tion  in  constraining  range  at  visible  wavelengths 
and  the  complementary  significance  of  non-aerosol 
extinction  contributions  in  the  infrared.  With 
respect  to  the  latter,  calculated  aerosol  extinc¬ 
tion  accounts  for  less  than  7%  of  the  total  at 
10. 6pm,  while  that  for  water  vapor  continuum  ab¬ 
sorption  is  75%.  While  visible  radiance  data  may 
be  useful  in  predicting  visible  range,  it  must 
be  augmented  by  absolute  humidity  data  from  other 
sources  to  evaluate  IR  range. 

10.3  Recommendations 

Although  limited,  the  validation  exercises  described 

in  this  report  associated  with  both  the  MAGAT  data  and 
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various  radiance  simulations  give  some  confidence  that  even 
though  simple  techniques  are  employed,  the  overall  simula¬ 
tion  is  capable  of  providing  realistic  depictions  of  the 
sensitivity  of  remotely  sensed  radiances  to  the  various 
meteorological  data  inputs.  Based  on  the  simulation  and 
validation  studies  reported  on  here,  it  is  appropriate  to 
recommend  the  computer  model  documented  above  as  a  potential 
tool  to  aid  in  the  analysis  of  DMSP  visible  imagery.  Exer¬ 
cising  the  code  in  the  course  of  imagery  analysis  provides 
the  opportunity  for  validation  through  experience.  Addi¬ 
tionally,  it  is  probable  that  the  data  available  from  future 
experimental  programs  will  provide  additional  opportunity 
for  testing  and  evaluation. 

In  the  course  of  the  study  it  became  clear  that  an 
analogous  imagery  simulation  code  would  be  potentially 
useful  in  the  analysis  of  corresponding  infrared  window 
(i.e.,  DMSP  T  L/F)  imagery.  Since  case  study  analysis  is 
most  often  performed  with  both  data  sets  side-by-side  in  an 
"eyeball"  multispectral  comparison,  it  is  reasonable  to 
suggest  quantification  of  the  observed  phenomenon.  Such 
a  combined  simulation  tool  may  be  helpful  in  determining 
the  difference  between  low  level  moisture ,  sea  surface 
temperature,  and  low  level  aerosol  effects. 

Another  rationale  for  recommending  further  simulation 
in  the  IR  window  is  the  need  to  quantify  low  level  water 
vapor  effects  to  supplement  the  IR  range  retrieval  problem 
discussed  previously.  Further  development  of  the  prototype 
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range  retrieval  algorithm  proposed  above  is  recommended 
including  investigation  of  the  potential  use  of  the  IR 
T/LF  data  to  provide  the  required  water  vapor  amounts  in 
the  near  surface  environment. 

Finally,  it  should  be  noted  that  the  simulation  code 
as  documented  is  configured  for  a  specific  purpose,  that 
of  interactive  DMSP  radiance  simulation.  With  suitable 
modification  the  basic  algorithm  will  provide  a  simulated 
wavelength  dependent  radiance  sensed  at  a  given  location 
for  an  input  data  vector  of  site  specific  meteorological 
variables.  Thus  the  potential  exists  to  use  the  basic  algo¬ 
rithm  as  a  module  within  a  larger  model  to  simulate  fields 
of  either  radiance  or  flux  or  in  an  inverse  mode  with  sup¬ 
plementary  input  to  infer  the  data  fields  from  radiance 
fields.  Transformation  fields  may  also  be  computed  to  sim¬ 
ulate  one  visible  sensor's  data  from  another's  radiance 
field.  These  uses  remain  to  be  explored. 
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APPENDIX  A:  DERIVATION  OF  APPROXIMATE  ANALYTICAL 

SOLUTION  TO  RADIATIVE  TRANSFER  EQUATION 

Following  Kaufman  (197  9)  ,  the  analytical  solution  to  the 
radiative  transfer  equation  (6-1,  6-2)  used  in  this  study  is 
derived  by  approximating  the  integral  term  in  the  source 
function  J  (6-2)  by  assuming  that  the  upward  (y  >  o) 
and  downward  (y  <  o)  diffuse  (i.e.  scattered)  intensity 
fields  are  independent  of  direction  (y,c|>)  and  given  by  I+ 
and  I~,  respectively.  Mathematically,  this  approach  may  be 
considered  as  a  two-point  quadrature  of  the  angular  integra¬ 
tion  over  upper  and  lower  hemispheres.  Physically,  the 
method  is  interpretable  as  an  approximation  of  the  multiply 
scattered  contributions  to  the  intensity  field.  For  incident 
solar  irradiance  irF,  Equation  (6-2)  is  approximated  as: 
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Upon  substitution  of  (A.l),  the  radiative  transfer  Equation 
(6-1)  is  written  separately  for  each  of  the  two  streams,  li¬ 
as  : 
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The  angular  integrations  in  (A.2a,b)  are  performed 
the  phase  function  normalization: 
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by  noting 


(A.  3) 


and  defining  the  azimuthally  symmetric  function,  P°;  the 
backscatter  fraction  for  monodirectioi’.al  radiation,  8;  and 
the  backscatter  fraction  for  isotropic  radiation,  8'  as: 


.  2tt 

p° (y / y ' )  =  2i ~  f  P(y*v'  <*$' 


3(y) 


-  sib  /  p0<>^')  aw 


8'  =  J  3  (y  ' )  dy ' 

o 

Applying  these  definitions,  the  intensity  Equations 
become : 


y  =  l+  -  l+  0)O  [1-6 (p)  3  -  l’o)0B(y> 


-  |  P(y  ,-y0;<j),(j>o)  e  T//y° 


-  y  =  l“  -  I+  03OB(y)  -  i"  <*><,[  1-P  'y)  ] 


J  P(y,y0;4>»4>0)  e”T//yo 


<• 


i# 


Assuming  that  the  scattered  diffuse  intensities  are 
by  the  fluxes  Fi(x)  such  that: 
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(A. 4) 


(A. 5) 


(A. 6) 


(A. 2) 


(A. 7a) 


(A. 7b) 


given 


I*  (t  ,y  ,4>)  =  F^tJ/tt 


(A. 8) 


Equation  (A. 7)  may  be  written  as: 


dl  _+ 

y  ar-  -  1  = 


-  ^2 
IT 


F  [ 1-6  (y) ]  +  F  3  (y ) 


J  P(y,-y0;4>,4>0)  e"T/yo 


dl 

y  a=r  •  1 


^0  l 

IT 


F'rB(u)  +  F  [l-S(y)  ] 


P(y ,u0;<j>,<J>0)  e“T//y° 


Integration  of  equations  (A. 9)  over  the  angle 
coupled  equations  for  the  fluxes: 


*  2F+  [  l-u)0  (1—6  ’ )  ]  -  2F“  o>03  '  -  TTFw0B(y0)  e 


gf-  -  -2F- [  1— u>0  (1—0  1 )  ]  +  2F+  oj03'  +  7TFcj0[l-B(y0 


(A. 9a) 

(A. 9b) 

yields 

-t/u0 

(A. 10a) 

n  e-T/,J» 

(A. 10b) 


Applying  the  boundary  conditions  [(6-4)  in  the  text]: 


F  (T=o)  =  0 


F+(T  =  T*)  =  [IT  u0  e”T  //y°  +  F~  ( T  =  T  ' 


)  ] 


(38) 


the  solutions  for  F—  are 


F+(t)  =  AekT  +  Be"kT  +  Ce'x/y° 


F  (T)  =  — 
^2 


A(Yi“k)ekT  +  B(y1+k)e  kT+  Ye“T</y° 


(A. 11a) 


>  (A. lib) 


where  the  appropriate  constants  are  given  by: 


YX  *  2  [l-u)0(l-e')  ]  (A. 12) 

Y2  =  2  u06'  (A. 13) 

k  =  (Yx2  -  Y22)  35  (A. 14) 

A  =  [B(Yx+k)  +  Y]/(k-Y1)  (A. 15) 

B  =  (Exekx  +  E2e"x  /y0)/(E3ekT  +  E4e“kx*)  (A. 16) 


C  =  tt  oj0 


-  Y1«(y0)  -  y2[1-0(uo)  ] 

x  (  u°2  ) 

V l-k^u0 v 


(A. 17) 
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3B 

i-Y 

l 


Y  =  c  <Yi  +  ^  )  _7rF  ^0&(U0) 

Ex  =  Y  [l/(Yl-k)  -  A2/y2] 

A^Y 

E2  =  [-  C  +  TTFy0  Al  +  — ] 

e3  =  (Yi  +  k)[l/(k-y1)  +  Al/y2] 
E4  =  [1  -  Al  (Yi+k)/Y2l 


(A. 18) 

(A. 19) 

(A.  20) 

(A. 21) 

(A. 22) 


* 

■ 


The  solution  for  the  upward  (i.e.  emergent)  intensity, 
I+ (x  ,y  ,({>)  ,  is  obtained  by  integrating  equation  (A.  9)  subject 
to  the  boundary  conditions: 


I  (t=o,  y ,  <j>)  =  0 


I+  (x=x* ,  y ,  c|) )  =  F+  (X=X*)/lT 


to  obtain  [compare  to  text  (6-5) ]  : 


(43) 


(44) 


kx  „  _-kx  „  _-x/y0 


U2e 


U3e 


I  (x  ,  y  , <j> )  DX  e  /y  +  1-yk  +  1+yk  +  i+y/y0 


(A. 23) 


where : 


Di  = 


7T 


Uxe 


kx*  „  -k 


U?e 


U^e 


-T*/y< 


1-yk  1+yk  '  l+y/y< 


-r*/y 


(A. 24) 


A-6 


AUJq 

IT 


[1-6  (U)  +  (Yi“k)  ] 


(A. 25) 


€ 


% 


!  •) 


u 


1 


OWq  Qf  \ 

U2  =  —  [l-e(u)  +  SlHl  (y1+k) ]  (A. 26) 

U3  =  —  [  1-6  (u)  ]  +  —  Y  +  |  P  (y,-M0;<j>,<j>0)  (A. 27) 

The  required  value  of  F+(x*)  in  (A. 24)  is  available  from 
equation  (A. 11a).  The  values  of  S',8(y)/  and  6(y0)  in  the 
above  expressions  are  obtained  from  Figure  6-3  after  calcu¬ 
lating  the  appropriate  asymmetry  parameter,  g,  using  equa¬ 
tion  (6-13) .  Values  of  x*  and  are  obtained  using  the 
modeled  optical  data  for  gases  and  aerosols  as  described 
in  Sections  3  and  5  of  Isaacs  (1980) . 
and  5,  respectively. 


-  • 
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APPENDIX  B:  PROGRAM  DESCRIPTION 


B . 1  Technical  Approach 

A  major  goal  of  this  research  program  has  been  the 
development  of  an  analysis  tool  relating  in  situ  meteor¬ 
ological  variables  and  remotely  sensed  radiances.  This 
section  describes  and  documents  the  research  grade  computer 
algorithm  developed  during  the  course  of  this  work  and  sub¬ 
sequently  applied  to  the  simulations  of  DMSP  VHR  and  LF 
imagery  presented  in  previous  sections.  The  general  tech¬ 
nical  approach  adopted  in  formulating  the  simulation  model 
code  is  summarized  in  Figure  B-l  reproduced  from  Isaacs 
(1980)  . 

Requisite  input  parameters  include  the  meteorological 
variables  characterizing  the  lower  atmosphere:  relative 
humidity  (RH) ,  wind  speed  (WS) ,  number  density  (N) ,  or 
visual  range  (Vr) ,  and  aerosol  scale  heights  for  both  fine 
and  course  aerosol  size  ranges  (H^,  H 2,  respectively).  The 
input  required  for  specification  of  the  sun /atmosphere- 
ocean/sensor  geometric  configuration  consists  of  spacecraft 
and  solar  ephemeris  data  including:  the  longitude  of  the 
ascending  node  (AN) ,  the  time  of  the  ascending  node  (TN) , 
the  solar  declination  angle  (SD)  for  the  desired  date,  and 
the  time  of  the  scan  line  desired  (T) .  The  sensor  zenith 
angle,  solar  zenith  angle,  and  azimuth  angle  difference 
which  were  previously  input  data  are  now  computed  based  on 
the  above  input  using  the  approach  described  earlier  in 
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Section  5.  The  results  obtained  for  each  set  of  input  param¬ 
eters  are  an  intensity  (or  radiance)  spectrum  for  wavelengths 

between  0.4  and  1.1  ym,  1(A),  and  a  DMSP  VHR  response  func- 

_  2 

tion  weighted  radiance,  I(DMS),  in  units  of  mW/cm  -ym-sr. 

The  process  linking  input  parameters  and  desired  output 
as  illustrated  in  Figure  B-l  may  be  divided  into  four  basic 
modeling  efforts:  (a)  atmospheric  transmission  (i.e. ,  ex¬ 
clusive  of  aerosols) ,  (b)  physical  modeling  of  aerosol  prop¬ 

erties  including  relative  humidity  and  wind  speed  dependence 
of  aerosol  size  distribution  and  index  of  refraction,  (c) 
determination  of  wavelength-dependent  aerosol  optical  prop¬ 
erties,  and  finally,  (d)  radiative  transfer  theory.  A  de¬ 
scription  of  the  approaches  to  incorporate  these  processes 
within  the  simulation  model  algorithm  is  described  in  Isaacs 
(1980) .  The  present  version  of  the  code  contains  several 
additions  and  changes  subsequently  implemented  and  described 
in  previous  sections  of  this  document.  These  include:  (a) 
incorporation  of  rural  and  urban  aerosol  models  (Section  4) , 
(b)  calculation  of  sun/sensor  geometry  directly  from  space¬ 
craft  and  solar  ephemeris  data  (Section  5) ,  and  treatment  of 
surface  reflection  in  potential  sunglint  regions  (Section  7 ) . 


•  B.2  Program  Attributes 

These  programs  simulate  the  radiances  potentially  avail¬ 
able  from  DMSP  VHR  and  LF  sensors.  Input  and  output  are 

•  .  formulated  to  support  meteorological  analyses  of  hard  copy 

imagery.  Calculations  include  the  following  features:  (1) 
evaluation  of  solar  subpoint,  satellite  subpoint,  and  scan 
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line  surface  element  pixel  location  (i.e.,  latitude/ longitude) 
based  on  ascending  node  time  and  location,  desired  scan  line 
time,  and  solar  declination  angle,  calculation  of  solar 
zenith  angle,  sensor  zenith  angle,  azimuth  angle  difference, 
and  azimuth  relative  to  prevailing  wind  direction,  (3)  effec¬ 
tive  ocean  surface  reflectivity  defined  using  sun/sensor 
geometry,  input  wind  speed  and  direction  based  on  the  Fresnel 
relations  and  Cox-Munk  sea  slope  statistics,  and  (4)  radiative 
transfer  calculations  based  on  the  above  to  simulate  sensor 
incident  radiances. 

As  currently  configured  (Version  81/3) ,  the  program  is 
coded  in  the  extended  BASIC  language  (Hewlett-Packard,  1979) 
and  is  implemented  on  a  desk  top  (HP- 85)  computer  system. 
Conversions  to  higher  order  languages  should  be  relatively 
easy  if  desired.  The  version  documented  here  is  in  a  user- 
oriented,  interactive  format  suited  for  individual  case  an¬ 
alyses  and  thus  requires  realtime  operator  intervention  to 
supply  input  data.  The  software  is  organized  into  a  main 
program,  three  subprograms,  and  twelve  supporting  data  files. 
These  are  identified  and  described  by  function  or  contents  in 
Table  B-l.  The  relationship  between  these  .program  elements 
is  grossly  summarized  by  the  program  logic  flow  chart  pic¬ 
tured  in  Figure  B-2. 

In  order  to  optimize  computational  efficiency  the  pri¬ 
mary  iteration  loop  is  within  the  RADTRA  subprogram  and  per¬ 
forms  the  radiative  transfer  calculations  necessary  to  simu¬ 
late  the  radiances  for  each  pixel  in  the  chosen  scan  line. 
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Program  (P) 


Name 

or  Data  (D) 

Bytes 

Description 

SLGP 

P 

256 

Scan  line  generator  (sun/ 
pixel/ sensor  geometry) 

RUFF  3 

P 

256 

Specular  surface  reflectivity 

RADTRA 

P 

256 

Optical  properties  and 
radiative  transfer 

FLUX 

D 

288 

Solar  flux 

FAD ATI 

D 

6912 

Tropospheric  aerosol  model 

FAD AT 2 

D 

6912 

Maritime  aerosol  model 

URBDAT 

D 

6912 

Urban  aerosol  model 

RURDAT 

D 

6912 

Rural  aerosol  model 

CA-OWS 

D 

4320 

M-K  coarse  mode  (0  ms 

CA-5WS 

D 

4320 

M-K  coarse  mode  (5  ms 

CA-7WS 

D 

4320 

M-K  coarse  mode  (7  ms-1) 

C A- 10WS 

D 

4320 

M-K  coarse  mode  (10  ms  ) 

LTRAN4 

D 

2880 

Clear  atmosphere  optical 
depths 

PHI 

D 

288 

DMSP  VHR  response  function 

ALBEDO 

D 

288 

Lambert  sea  surface  albedo 

Table  B-l.  Simulation  Model  Software 
File  Attributes  and  Purpose 


FIGURE  P-2:  SIMULATION  MODEL  FLOW  CHART 


Before  entering  this  loop,  sun/ sensor  geometry  is  evaluated 
(in  SLGP)  for  the  entire  scan  line  supporting  both  the  deter 
mination  of  effective  reflectivity  (in  RUFF 3)  for  each  pixel 
and  the  transfer  calculation.  Optical  data  files  are  read 
once  before  the  main  iteration  loop.  As  a  consequence  of 
this  configuration,  input  variables  may  not  change  along  a 
scan  line.^ 

B. 3  Program  Listings 

Current  listings  of  the  simulation  model  computer  code 
are  given  in  Tables  Br-2  through  B-5.  Included  are  the  main 
calling  driver  (MAIN,  Table  B-2a,b)  and  the  three  subroutine 
programs  identifed  above  in  Table  B-l.  These  are:  the  scan 
line  generator  program  (SLGP,  Table  B-3a,b) ,  the  wind-driven 
ocean  surface  model  (RUFF3,  Table  B-4a,b,c,d) ,  and  the  radi¬ 
ative  transfer  program  (RADTRA,  Table  B-5a,b,c,d,e) . 

B.4  Sample  Calculation 

As  a  guide  to  program  operation,  this  section  presents 
a  sample  calculation  from  the  user's  vantage  point.  For 
illustrative  purposes,  it  is  assumed  that  the  program  code 
(as  listed  in  Section  B.4)  and  the  associated  data  files 
(see  Table  B-l)  are  stored  in  the  cartridge  format  and  the 
program  is  run  on  an  HP- 85  computer.  The  case  analyzed  in 
this  example  is  identical  to  that  discussed  in  Section  7.7 

This  restriction  is  easily  removed  by  reprogramming. 
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1  !  :r**:m<EPRF/ONR-MfiIN******* 
Main  P rog . -Sa t e 1 1 i t e  Imagery 
**R .  Isaacs/ftER* Inc . /lSSl** 

2  !  ***********  *************** 

3  !  ******Version  8 \ /3******** 

4  !  These  programs  simulate  th 
e  imagery  potentially  availa 
ble  from  the  OMSP  VHR  and  LF 
sensors . 

5  !  Calculations  include  the  f 
ollowing  features: 

6  !  <1)  evaluation  of  solar  su 

bpoint*  satellite  subpoint* 
and  scan  line  pixel  location 
< 1  at/ Ion) 

7  !  based  on  ascending  node  ti 
me  and  1  oca t i on* des i red  scan 
line  time  and  solar  declinat 
ion  angle 

8  !  <2)  calculate  solar  zenith 

ang 1 e * senso r  zenith  angle* re 
lative  azimuth* and  azimuth  d 
i f  f  erence 

9  !  with  respect  to  prevailing 
wind  d i rec t i on . 

10  !  <3)  based  on  sun/sensor  ge 
ometry  *  input  wind  speed  and 
wind  direction*  an  effective 
ocean 

11  !  surface  reflectivity  is  co 
imputed  using  Cox-Munk  sea  si 
ops  statistics  and  Fresnel  r 
elat ions . 

12  !  <4)  computed  geometry  and 

surface  reflectance  are  inpu 

.  t  to  radiative  transfer  calc 
ulat ions . 

13  !  The  radiative  transfer  sch 
erne  is  described  in  Isaacs*! 
980  <NflVENVPREDRSCHFflC  CR  80 
-06) 

14  !  Due  to  the  approx  nature  o 
f  the  calculation*  results  mu 
st  be  carefully  interpretedC 
see  text) 

15  !  Required  are  meteorolosica 
1  data  to  characterize  the  o 
pt icalpropert ies  of  the  atmo 
sphere 

16  !  including  relative  humidit 
y*aerosol  scale  heights<fine 
ycoarse)*and  choice  of  one  o 
f  four 

17  !  aerosol  mode  1 s : < t roposphe r 
ic*mari time*  rural*  urban) 

18  !  <5)  Results  include  radian 

ces<mW/cm**2-.um-sr )  in  <0.4- 
l.l^m) range  and  DMSP  weishte 
d  values 

19  !  for  each  selected  scan  lin 
e  pixel  . 

Table  3-2a.  Listing-MAIN  Program 


<0 


<» 


20  !  X************************* 
31  !  Software  is  organized  into 
three<3)  subprograms(P)and  t 
wel ve< 12)data  files <D)  liste 


c 

f  below1 

32 

SLGP 

p 

scan  line  genera 

33 

tion  program 

34 

RUFF3 

p 

specular  surface 

35 

reflectivity 

36 

RADTRA 

0 

optical  proper 

37 

ties  and  radia 

33 

/ 

tive  transfer 

39 

FLUX 

0 

solar  flux 

40 

FADAT1 

D 

t rop  aer.  mod. 

41 

FADAT2 

D 

mari t .  aer .  mod . 

42 

URBDAT 

D 

urban  aer.  mod. 

43 

RUROAT 

0 

rural  aer.  mod. 

44 

CA-0WS 

D 

coarsemode(0m/s) 

45 

CA-5WS 

D 

coarsemode(5m/s) 

46 

CA-7WS 

D 

coarsemode (7m/s) 

47 

CA10WS 

0 

coarsemode< 10m/s 

48 

LTRAN4 

0 

clear  atm.  model 

49 

PHI 

D 

DMSP  response 

50 

ALBEDO 

D 

sea  surf .  albedo 

51 

(Lambert  reflect 

52 

7 

(U  ^  ^  > If  4/  ^  ^  ^  ^  ^  \Lf  ^  ^  ^  <1/  -a  ^  *1/  ^  ^  ^  ^ 

^p* 

79  PRINT  “********************* 

11 

^  ^  ^  ^r*  g*  ^ 

80  PRINT  "tOMSP  VHR/LF  Imagery 
Simulat ion#" 

81  PRINT  "*  R.G. ISAACS/RER/ I 

nc./1981  *" 

32  PRINT  "*  Version  31/ 

3  *" 

33  PRINT  "X******************** 

txxtttxxxxt” 

90  PRINT  "  RADIANCES  FOR  SCAN  L 
INE  AT= 

110  OPTION  BASE  1 

120  COM  U (11)/ U0 (11)/ D3 <11) 

130  COM  W9 

140  COM  F2( 1 1 ) »  Y9 

150  Y9=0 

160  CHAIN  " SLGP" 

170  END 


'§ 


t 


Table  B-2b.  Listing-MAIN  Program  (continued) 
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1  !  *X**1*NEPRF>'0HR-&LGP**X4:X* 
*Scan  Line  Generator  Pros  4:* 
**R  .  Isaacs^AER,  Inc  .  /1381*** 

2  DISP  "Scan  Line  Geometry  Pro 
sram;  Assumes  sun  sync  orbit 

w/  al t . =833 km, per . =181 . 35m 

M 

3  DEG 

4  OPTION  BASE  1 

5  DISP  "LONG<RN>=Des" ; 

6  INPUT  LO 

7  COM  U  <  1 1  >  ,  IJ0 <  1 1  >  ,  D8 <  1 1  ) 

3  COM  M3 

3  COM  F2< 1 1 > , Y3 

10  DISP  “AN<GMT>=HH, MM, SS" ; 

11  INPUT  H0 ,  M0 , S0 

12  Y9=l 

15  DISP  “Desired  T  (GMT  v  =HH  ,  MM ,  S 
S 11  / 

16  INPUT  HI, Ml, SI 
20  P0=101 . 35 

25  DISP  "Solar  Decl  Ans  = 

26  INPUT  03 

27  DISP  "Scan  Line  Print?  < 1 > Ye 
sy<0>No" 

23  INPUT  P7 

30  PRINT  ,,L0NG<AN)  =  "  ,L0 

31  PRINT  "  AN(.GMT>  ="  ,  h0,  M0,  SO 

32  PRINT  " Scan < GMT >=" , HI , Ml , SI 

40  !  U  *  ARRAY  OF  SAT  ZEN  ANG  C 
OS 

41  !  U0=  ARRAY  OF  SUN  ZEN  ANG  C 
OS 

42  !  D8  =  ARRAY  OF  AZI  PUG  DIFF 
90  !  GENERATE  Sat. Sub  Point  T ra 

ck 

100  T0=H0+M0/60+S0y3600 
110  Tl=Hl+Ml/'60  +  Sl/'3600 

115  L9=<T1-T0>*15 

116  LS=<T1-12>*15 

US  PRINT  “Sol  Decl.  =  ",D9 

119  PRINT  "Sol  Lons .  =  ",L8 

120  P=P0^60 

130  C=<T1-T0>*360/P 
140  B=90 
150  A  1=3.7 

160  T1=C0S< .5*<B-C>>*SEC< .5*<B+C 
>>*C0T<  .  5*A 1 ) 

170  T2  =  SIN<  .  5*<B-C>  >*CSC<  .  5*<B+C 
) ;*C0T  < .5iftl) 

ISO  B1=ATN<T1 >+ATN<T2> 

130  C1=ATN<T1 >-ATN<T2> 

200  T3=TAN< . 5*<B-C>  >*SIN< . 5*<B1+ 
C1>)*CSC< . 5*  <  B 1 -C 1 >  > 

210  A=2*ATN<T3;> 

220  L 1 =90-A 
230  L2=L0+C 1 +L9 

235  PRINT  "Sub.  Pt.  Lat/Lons  =  ", 
LI, "y " , L2 
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500  !  SCAN  LINE  GENERATOR 

510  E=90-Ll 

520  01=81 .3 

530  FOR  N= 1  TO  11 

540  F=<N-1 >*13 . 303/10 

550  T5=CGS< . 5*(E-F>>*SEC< . 5*<£+F 

>  >*COT< . 5*D  1 > 

560  T6=SIN< . 5*<E-F) )*CSC< . 5*<E+F 

>  >*COT< . 5*01  > 

570  L 1 =ATN<  T5>  +  ATNCT6 > 

580  F1=ATN<T5>-ATN<T6> 

590  D=2*ATN<  TAN< . 5*<E-F> >*SIN<  .  5 
*< E 1+F 1 >  > *C3C < . 5*<E 1 -F 1  )  >  > 
600  L3-90-Q 
610  L4=L2-F1 
615  IF  P?=0  THEN  630 

620  PRINT  "Scan  Lin©  Pixel  =  %N 

621  PRINT  "Pixel  Lat/Lon©  =“,L3, 
"/" ,L4 

630  X=<N-1>* 1479/10 

631  2=ATN<X/833> 

750  !  SOLAR  POSITION  DATA  FOR  SC 
AN 

751  !  D9=DECLINATI0N  ANGLE 

752  !  H9=SUNS  HOUR  ANGLE 

753  »  29=S0LAR  ZENITH  ANGLE 

754  !  A9=S0LAR  AZIMUTH  ANGLE 
762  H9=L4-L3 

770  Z9=90-ASN<SIN<L3>*SIN<D9>+CO 
S<L3>*C0S(D9>*C0S(H9>> 

730  Q 1 =90-L3 

782  Q2=Z9 
784  Q3=90-D9 

736  Q4=. 5*<Q1+Q2+Q3> 

783  Q5=SIN<Q4-Q2>*SIN<Q4-G1>/«.SI 
NtQ2>*SIN<Ql >  > 

790  A9= ACS  < 1 -2*Q5 ) 

795  IF  H9 <0  THEN  A9=360-R9 
808  IF  P7=0  THEN  900 

801  PRINT  "Sat.  Zen  An©  =  ",Z 

802  PRINT  “Hour  fine.  =  "  ,H9 

803  PRINT  "Sol  Zen  An©  =“,29 

804  PRINT  "  Sol  Azm  An©  =",A9 

900  U<N>=C0S<Z> 

901  U0<N>=COS<Z9> 

902  D8<N)=261 . 3-A9 
1000  NEXT  N 

2000  PRINT  "*********0CEAN  SURrA 
CE********‘‘ 

2001  PRINT  -  REFLECTAN 

CE" 

3000  CHAIN  "RUFF3" 

9999  END 
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10  !  ****NEPRF/0NR-RUFF3**** 
**Wind-d riven  ocean  surf,  tt 
tt  R . Isaacs/RER, Inc . /1982  tt 
20  !  tttttttttttttttttttttttttt 
tt 

38  OPTION  BRSE  1 
40  COM  U<11>,U0(11>,D8<11>,W9,F 
2< 1 1 > , Y9 

50  OISP  “WIND  SPEED(M^S) = i(0.0> 
,2(5.0), 3(7. 0>, 4(10. 0>"i 
60  INPUT  W9 

70  PRINT  "Wind  Speed  =  ",W9 
80  ON  M9  GOTO  90,110/130,150 
90  W= . 005 
100  GOTO  160 
110  W=5 
128  GOTO  160 
130  W=7 
140  GOTO  160 
150  W= 1 0 

160  OISP  “WIND  RZ . < ie  N=000,E=09 
0 , S= 1 80 , W=270  >  " ; 

170  INPUT  W1 

180  PRINT  "Wind  Di r . (de*>=“ , W1 
198  FOR  Q=1  TO  11 
280  DEG 

210  T0=RCS(U0(Q>> 

220  T2=ACS(U(Q>> 

230  D=D8(Q) 

240  RRD 

250  B=Wl+98 . 7 

260  ?  ttttt  DEFINE  ttttttttttt 
270  !  Z1  =  UPWIND  SLOPE  COMP(ZU) 
280  !  Z2  =  CROSSWIND  SLOPE  COMP< 
ZC) 

290  !  U8  =  SOLRR  ZENITH  RNG  COS 
300  »  T0  =  SOLRR  ZENITH  RNG 
310  !  D  =  SOLRR  AZIMUTH  DIFFER 

ENCE 

320  !  B  =  AZIMUTH  DIFF.  BETW . 

SUN  AND  WIND 

330  *  T2  =  OBSERVER  ZENITH  ANGLE 
340  !  VEC(n)=SURF.  NORM. 

350  !  VEC(s)=SOLAR  DIR. 

360  !  VEC(o)=OBSERV.  DIR 
370  ‘  W  =  WIND  SPEED  (M/'S) 

380  U0=COS(DTR(T0) ) 

390  U=C0S(DTR(T2) ) 

400  PRINT  U0,U 

410  \  CARTESIAN  COORDS.  OF  S0 
420  S1=SQR<  1-U0“*2>*COS(DTR<  188+B 
>) 

430  S2=SQR<1-U0A2)*SIN(DTR(180+B 

)) 

440  S3=U0 

450  !  SPHERICAL  COORDS.  OF  N0 
460  D0=180+B 

470  !  CARTESIAN  COORDINATES  OF  0 
BSERVER  DIRECTION 
480  O1=SIN(DTR(T2))*COS(DTR(180+ 
B+D)  ) 
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430  O2=SIN<OTR<T2>>*SIN<DTR<180+ 
B+D>> 

500  03=U 
510  D2=180+B+D 

520  !  *****  ROTATE  AXES  ***** 

530  !  EVALUATE  z^‘=<s  x  o  ) 

540  C1=S2*03-02*S3 
550  C2=01*S3-S1*03 
560  C3=S1*02-S2*01 
570  !  COMPONENTS  OF  x~'=<Si,S2,S 
3) 

580  R 1 =S 1 
590  R2=S2 
600  R3=S3 

610  !  COMPONENTS  OF  x^‘=(zA,)x<x 
A '  ) 

620  B1=A3*C2-A2*C3 

630  B2=A1*C3-A3*C1 

640  B3=A2*C1-A1*C2 

650  *****  DOT  PROD.  CHK .  ***** 

660  A=A1*B1+A2*B2+A3*B3 

670  B=B1*C1+B2*C2+B3*C3 

680  C=A1*C1+A2*C2+A3*C3 

630  !  *****  MAG.  CHK . ***** 

700  E=SQR< A1 a2+A2a2+A3a2) 

710  F=SQR<B1 a2+B2a2+B3a2> 

720  G=SQR<C1 A2+C2A2+C3A2> 

730  !  *****  NORMALIZE  ******* 

740  Al=AlxE 
750  A2=A2/E 
768  A3=A3/E 
770  B 1 =B 1 xF 
788  B2=B2/F 
730  B3=B3xF 
800  C1=C1/G 
810  C2=C2xG 
820  C3=C3/G 

330  H=SQR< A1/s2+A2a2+A3^2> 

840  I=SQR<Bl~2  +  B2-'2  +  B3A2> 

850  J=SQR<C1^2+C2A2+C3A2> 

860  !  ***COMPUTE  <  BETW .  s-  and 
=  A0  ***** 

870  A0=RTD<ACS<O1*S1+O2*S2+O3*S3 

>  ) 

880  !  *****COMPS.  OF  n~ ' ******* 
830  Nl=COS<DTR<A0x2> ) 

300  N2=SIN<DTR<A0/2>> 

310  N3=0 

320  !  *****ROTATE  BACK  ******** 

930  M1=A1*N1+B1*N2 
940  M2=A2*N1+B2*N2 
950  M3=A3*N1+B3*N2 
960  !  *****  SPHERICAL  REP . ****** 
970  T3=RT0 <ACS<M3XSQR< Ml ''2+M2-2  + 
M3~2>  >> 

980  D3=RTD<ATN<M2xM1 > > 

990  !  *****FINO  SLOPES,  Z1 >  Z2**** 
* 

1000  Z1 =-M 1 xM3 
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1010  22=-M2/M3 

1020  IF  ABS<21 >  < . 001  THEN  1030  E 
LSE  1040 
1030  21=0 

1040  IF  flBS<22>  < . 001  THEN  1050  E 
LSE  1060 
1050  22=0 
1060  »  CONTINUE 
1070  I  ****  FRESNEL  REFL . **** 

1 030  Q 1 = 1 
1090  Q2= 1 . 333 
1100  fl4=OTR<FI0/2) 

1110  A5=fiSN<Ql*SIN<fl4>^Q2> 

1 120  Ul  =  <TflN<fl4-A5>/'TAN<fl4+fl5>  >•' 
2 

1  130  U2=<-SIN<fl4-A5>/SIN<R4+A5>  > 
-2 

1140  R9= . 5*<U1+U2> 

1150  R6=fl0^2 

1160  !  ***COX-MUNK  STRTS . ***** 
1170  G1=SQR< .00316*14) 

1130  G2=SQR< .003+.  00192*14) 

1 190  P9=EXP<-  .  5*<<21/Gl)-2+<22,'G 
2)~2> )/<2*PI*Gl*G2> 

1200  P5=<U2-U1 )/(U2+Ul ) 

1210  !  ****OEFINE  EFFECTIVE***** 
1220  !  **SPECULRR  REFLECTION  ** 
1230  !  F2  =  Effective  Fresnel 

Ref lectivitx 

1240  !  F2  =<Fresnel  Ref)*CProb 

of  2u*2c  from 

C-M) 

1 258  F2  < Q >  =R9*P9<'(4*M3M*U < Q )  > 
1260  PRINT  -PIXEL  #=",Q, "Rf <ef f e 
ct ive)=" , F2<Q> , "POL  .  =" ,  P5 
1270  NEXT  Q 
1236  CHAIN  "RROTRfl" 

1290  END 


L 
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1  !  ******NEPRF/ONR-RflDTRfi**** 
**Radiative  Transfer  Pros*** 
**R  .  Isaacs/AER,  Inc . /'1982*** 

2  PRINT  " *********RflDI ANCES*** 
*******“ 

91  OPTION  BASE  1 

92  COM  U<11),U0<11),D8<11),W9,F 
2<11), Y9 

97  DIM  U9<2),A9<10,5),U8<:9),B7< 

2) 

98  DIM  T1 <36) , T2<36) , T3<36) , T4< 

36) , T5<36) , T6<36) , G3<36),G4< 

36) 

99  DIM  I1<36),F9<36),A2<36),P9< 

36) 

100  DIM  L9 < 36 ) 

103  DISP  "REL .  HUMIDITY, RH<7>1< 

07.)  ,  2 <507)  ,  3 <707)  ,4<807)  ,  5 <9 
07) ,  6057)  ,  7  087) ,  3  097)"  ; 

104  INPUT  R9 

107  DISP  "AEROSOL  SCALE  HTS.CKM) 
=FINE, COARSE" 

103  INPUT  HI ,  H2 

110  DISP  "AEROSOL  MODEL =< 1 ) TROPG 
SPHERIC,  <2)MARITIME,  •::3)RURAL 
, <  4  >  URBAN  H i 

111  INPUT  Q9 

123  PRINT  "REL.  HUM . CRH<7)D=" , R9 

125  PRINT  “AEROSOL  SCALE  HTS(KM) 

FINE  - ",  HI , "COAR 

SE : - , H2 

126  ON  Q9  GOTO  127,129,131,133 

127  PRINT  "TROPOSPHERIC  AEROSOL" 

128  GOTO  134 

129  PRINT  "MARITIME  AEROSOL" 

130  GOTO  134 

131  PRINT  "RURAL  AEROSOL" 

132  GOTO  134 

133  PRINT  "URBAN  AEROSOL" 

134  ! 

140  !  READ  DATA  FILES  FROM  TAPE 

3 

141  !  SOLAR  IRRADI ANCE  -F9 

142  ASSIGN#  1  TO  "FLUX" 

143  FOR  1=1  TO  36 

144  READ#  1  ;  F9  < I > 

145  NEXT  I 

146  ASSIGN#  1  TO  * 

147  !  RESPONSE  FUNC .  -  P9 

148  ASSIGN#  2  TO  "PHI" 

149  FOR  1=1  TO  36 

150  READ#  2  ;  P9<I) 

151  NEXT  I 

152  ASSIGN#  2  TO  * 

153  !  FINE  AEROSOL  DATA-T2, T5, G3 

154  ON  Q9  GOTO  155,157,159,161 

155  ASSIGN#  3  TO  "FADAT 1 " 

156  GOTO  165 

157  ASSIGN#  3  TO  "FADAT 2 " 

158  GOTO  165 
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159  ASSIGN#  3  TO  "RURDAi" 

160  GOTO  165 

161  ASSIGN#  3  TO  "IJRBDAT" 

165  FOR  1=1  TO  8 

166  FOR  J=1  TO  36 

167  READ#  3  ;  A,B,C 
163  IF  IOR9  THEN  172 

169  T2<J>=A*H1 

170  T5<J>=B*H1 

171  G3< J}=C 

172  NEXT  J 

173  NEXT  I 

174  ASSIGN#  3  TO  % 

175  !  SURFACE  ALBEDO  -  A2 

176  ASSIGN#  4  TO  “ALBEDO" 

177  FOR  1=1  TO  36 
17S  READ#  4  j  A2<I> 

179  NEXT  I 

18©  ASSIGN#  4  TO  * 

181  !  LOWTRRN  DATA  -T1,T4 

182  ASSIGN#  5  TO  " LTRRN4 *' 

133  FOR  1=2  TO  6 

134  FOR  J=1  TO  36 

185  READ#  5  ;  A,B 

186  IF  IOR9  THEN  190 
137  T 1 <  J ) =  A 

188  T4  <  J  >  =B 

190  NEXT  J 

191  NEXT  I 

192  ASSIGN#  5  TO  * 

208  !  COARSE  MODE  RER0S0L-T3, T6, 

G4  -  TROPOSPHERIC  ONLY 

201  IF  Q9# 1  THEN  300 

202  ON  W9  GOTO  203/205/207,209 
283  ASSIGN#  6  TO  MCA-0WS" 

204  GOTO  210 

205  ASSIGN#  6  TO  "CR-5WS" 

206  GOTO  210 

207  ASSIGN#  6  TO  "CR-7NS" 

208  GOTO  210 

289  ASSIGN#  6  TO  "CR10WS" 

210  FOR  1=2  TO  6 

211  FOR  J=1  TO  36 

212  READ#  6  ,  A,B,C 

213  IF  IOR9  THEN  217 

214  T3< J)=A*H2 

215  T6<J>=8*H2 

216  G4< J)=C 

217  NEXT  J 
213  NEXT  I 

219  ASSIGN#  6  TO  * 

299  GOTO  315 

300  REM  DUMMY  ARRAY  VALUES 

301  FOR  J=1  TO  36 
304  T3 <  J  >  =0 

307  T6<J>=0 
389  G4< J>=0 
313  NEXT  J 

315  DISP  “ALL  FILES  READ" 

320  FOR  2=1  TO  11 
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330 

34fe 

341 

342 

345 

346 

347 

348 

349 

350 
360 
370 

390 

391 

392 

393 
399 
414 
416 

418 

420 

425 

435 

440 

442 

443 

444 

445 

446 

447 

448 

500 

501 

502 

503 

504 

505 

506 

507 

508 

509 


510 

699 


A8=F2<Z> 

'0=U0<Z) 

_  =U  <  2  ) 

D-D8<2)-180 

PRINT  "PIXEL  #  =%Z 

PRINT  "OBS.  ZENITHCv)  =  ",U 

PRINT  “SUN  ZENITH< vo>  =  “,U0 

PRINT  “AZIMUTH  ANG .<*>  =  “,  D 

PRINT  "Rf(e< f . >  =  “, A8 

FOR  L  =  1  TO  36 

IF  L>  =  16  THEN  370  ELSE  396 

A8=0 

T0=T1<L)+T2<L>+T3<L) 

W0=<T5<L>+T6<L>>^<T2<L)+T3<L 

>> 

G=<T5<L>*G3<L>+T6<L>*G4<L>>/ 
<T5<L>+T6<L) ) 

A1=A2<L) 

F=F9<L])/PI 

D=DTR<D> 

C1=-U*U0+SQR< 1-U-2)*SQR< 1-U0 
-2>*C0S<0> 

C2=U*U0  +  SQR<  l-U-B^SGRCl-UO- 
B^COSCO) 

S=ACS<-C2) 

S1=RTD<S> 

Pl=W8*<l-G-2>/< 1+G-2~2*G*C1 > 
-1 .5 

P2=W0*< l-G-2>/<l+G-2-2*G*C2> 
-1.5 

Pl  =  <T4<L>*.75*(l+C'l-2)  +  <  T5  <  L 

>+T6<L>>*Pl>/T0 

P2=<T4<L>*.?5*<l+C2-2>  +  <T5<:L 

>+T6<L>>*P2>/T0 

W0=<T4<L)+T5<L>+T6CL> )/T0 

G=<T5<L>*G3<L>+T6<L>*G4<L>>/ 

<T4<L>+T5<L>+T6<L>> 

!  H-G  PHASE  FUNCTION 
G0SU8  2000 
GOSUB  2300 
Q=4*< l-G-2*< 1-U0)> 

Gl  =  <7-3*G-2-W0*<4+3*G;>+W0*G- 

2*<4*B0+3*G>  >/Q 

G2=-< l-G-2-W0*<4-3*G>-W0*G-2 

*<4*B0+3*G-4>>/Q 

K=SQR<Gl-2-G2-2> 

C=PI*F*W0*<B0/U0-G2*< 1-B0>-G 

l<B0)<U0‘'2/<  l-K-2*U0-2> 

E1=<C*<G1+1/U0>-PI*F*W0*B0>* 

<1/<G1-K>-A1/G2> 

E2=-C+PI*F*U0*A1+A1/G2*<C*<G 

1 + 1 /U0 ) -PI *F*W0*B0 ) 

E3=<G1+K>*<-1/<G1-K>+A1/G2> 

E4=1-A1/G2*<G1+K> 

B  =  <E1»EXP<K1:T0)+E2*EXP<-T0/U 
0>>/<E3*EXP<K*T0>+E4*EXP<-K* 
T0)> 

A=<B*<G1+K>+<C*<G1+1/U0>-PI* 

F*W0*B0>>/<K-G1> 

I2=A8*PI*F*EXP<-T0/U0>*EXP<- 
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701  REM  UPWARD  INTENSITY 

710  U^=Rt<U0>/PI  >*(  1-B2+B2/'G2*<G1 

711  U2=B*<W0^PI>*<1-B2+B2/G2*<G1 
•HO  > 

712  U3=C*<W0,'PI>*<1-B2>+W0*B2,'<P 
I*G2>*<C*<G1+1/U0>-PI*F*W0*B 
0>+F/4*Pl 

713  F5=A*EXr<K*T0>+B*EXP<-K*T0>+ 
C*EXP(~T0/'U0) 

714  D2=F5/'PI-Ul*EXP<K*T0>/<l-U*K 

715  O3=-U2*EXP<-K*T0>^<l+u*tO-U3 

*EXP  < -T0'U0 > ✓ < 1 +U/U0 > 

716  O1=EXP<-T0*'U>*<D2+D3> 

719  T =8 

720  I1<L>=01*EXP<T/U>+U1*EXP<K*T 
>/<l-U*K>+U2*EXP<-K*T>/<l+U* 

K> +U3*EXP < -T/U0 > ✓ <1 +U/U0 > + I 2 

721  L9<L>=  .4+CL-1X.02 

722  PRINT  L9<L), I1<L>* 12 

749  RESTORE 

750  NEXT  L 

800  !  COMPUTE  DMSP  WEIGHTED  INTE 
NSITIES 

801  S 1 =0 

802  S2=0 

803  FOR  1=1  TO  36 

804  S1=S1+P9< I)tll(I) 

805  S2=S2+P9< I > 

806  NEXT  I 

807  I9=S1/S2 

810  PeRnI?Tty"^HDPRSS-WEIGHted  INT 

901  NEXT  2 

1999  GOTO  9996 

2000  REM  OATR  FOR  H-G  B8,B1,B2 
2930  FOR  Q=1  TO  10 

2040  FOR  S= 1  TO  5 
2050  READ  A9<Q,S> 

2060  NEXT  S 
2070  NEXT  Q 

2080  DATA  .  5, 0, 0, 0, 0,  .  5,  .  13975,39 
2889* -1 . 5989873995#  3.518411 
6963,-2.5592172257* .5*-. 120 
19000482 

2890  DATA  -.2724219928, .63859601 
23# -. 7459840146, .5, -.460171 
23188,1.18747390274,-2  0812 
230849 

2100  ORTA  .85392041407, .5,-  4868 
28796532* .49409050834,-1  01 

44699491. . 4272082373. .5,-  3 
001541656 

2110  DATA  -.35928947292, . 1589475 

781. . 00049606046. .5,-55347 

4414782. . 37397957743,-. 7476 
1782865 

2120  DATA  .42711266606* .5*-. 6267 

9466391 1 ..  18057740986, -  374 

16958202. . 32038683614,  5,- 
84678101 

Table  B-5d.  Listing-RADTRA  Subprogram  (continued) 
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<# 


« 


# 


2130 

DflTfi  .  50718036245..  -  .  3740740 

109, .2136746594, 

.5,-406823 

676662, .01406832224, . 105560 

7702 

2140 

DATA  -.2128054157 

2150 

FOR  1=1  TO  9 

2160 

READ  U8<  I  ) 

2170 

NEXT  I 

2180 

DATA  0,  . 1,  .2,  .3, 

.4, .5, .6, .8 

2200 

J 

RETURN 

2300 

REM  BACKSCATTER 

INTERP 

2310 

U9  < 1 >=U 

2311 

U9<2)=U0 

2320 

FOR  R= 1  TO  2 

2330 

IF  U9 < R >  < . 1  THEN 

2500 

2340 

IF  U9<R>  < . 2  THEN 

2502 

2350 

IF  U9CRX.3  THEN 

2504 

2360 

IF  U9<R>  < . 4  THEN 

2506 

2370 

IF  U9  <  R  X . 5  THEN 

2508 

2380 

IF  U9<RX  .  6  THEN 

2510 

2390 

IF  U9<RX .  8  THEN 

2512 

2400 

IF  U9<R)<=1  THEN 

2514 

2500 

Q=  1 

2561 

GOTO  2515 

2502 

<3=2 

2503 

GOTO  2515 

2504 

0=3 

2505 

GOTO  2515 

2506 

0=4 

2507 

GOTO  2515 

2508 

Q=5 

2509 

GOTO  2515 

2510 

Q=6 

2511 

GOTO  2515 

2512 

0=7 

2513 

GOTO  2515 

2514 

0=8 

2515 

Y5=0 

2516 

Y6=0 

2517 

Y7=0 

2518 

Q1=Q+1 

2525 

FOR  S= 1  TO  5 

2526 

Y9=A9<Q,S)*GMS- 

1  > 

2527 

Y8=A9<Q1,S)*GMS 

-1  ) 

2528 

Y3=A9<10,S)*GXS 

-1  > 

2529 

Y7=Y7+Y9 

2530 

Y6=Y6+Y8 

2531 

Y5=Y5+Y3 

2532 

NEXT  S 

2540 

B7<R)=Y7+<Y6-Y7)*<U9<R)-U8t 
QXXUSCOl  >-U8<Q>> 

2541 

NEXT  R 

2550 

B 1  =Y5 

2551 

B2=B7<1) 

2552 

B0=B7<2> 

2560 

RETURN 

9996 

PRINT  USING  9997 

9997 

IMAGE  A/ 

9998 

CHAIN  -MAIN- 

9999 

END 

Table  B-5e.  Listing-RADTRA  Subprogram  (continued) 


# 
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with  the  exception  of  the  input  wind  speed  of  7.0  ms’^  (radi¬ 
ances  in  §7.7  are  simulated  using  5.0  ms"^) . 

Running  the  program  from  tape  commences  by  loading  the 
main  program  (LOAD  "MAIN")  followed  by  a  run  command  (RUN 
key) .  As  the  printer  prints  the  heading  illustrated  in 
Figure  B-4,  the  SLGP  subprogram  is  loaded.  Figure  B-3 
illustrates  the  CRT  display  during  the  running  of  the  pro¬ 
gram.  Items  preceeded  by  an  asterisk  and  in  parentheses 
indicate  operations  during  which  the  CRT  is  blanked  out. 
Displays  followed  by  question  marks  (i.e.,  "?")  require 
operator  input.  Data  entered  in  response  by  the  user  (fol¬ 
lowed  by  an  END  LINE  key)  are  indicated  on  subsequent  lines. 

After  SLGP  is  loaded/  the  display  commencing  in  Figure 
B-3a  asks  for  satellite  and  solar  ephemeris  data  to  accom¬ 
plish  the  scan  line  geometry  calculation.  Data  required 
are:  (1)  the  longitude  (in  degrees  W)  of  the  ascending  node 

(AN)  ,  (2)  the  GMT  time  of  the  ascending  node  (TN)  in  hours, 
minutes  and  seconds  (not  input  format) ,  (3)  the  GMT  time  of 

the  desired  scan  line  (in  same  input  format) ,  and  (4)  the 
solar  declination  angle  (SD)  for  the  desired  date.  The  user 
has  the  option  of  selecting  a  scan  line  print  giving  the  lat¬ 
itude  and  longitude  (in  decimal  degrees)  and  other  information 
pertaining  to  each  of  the  eleven  available  pixels.  Figure  B-4 
illustrates  the  subsequent  printer  output  from  SLGP  for  the 
input  data  shown  in  Figure  B-3a  and  selecting  a  scan  line 
print.  (If  a  scan  line  print  is  not  selected,  output  is  not 
given  for  each  pixel,  but  the  calculations  are  done  and  stored 
for  later  use.) 


(loading  SLGP) 


* 


(a)  Scan  Line  Geometry  Prosram:  fissu 
mes  sun  sync  orbit  w/  ait.=833k 
per  .  =101 , 35m 
LONG<AN)*Des? 

94.0 

AN<GMT)=HH,MM, SS? 

15,00,08 

Desired  T<GMT)=HH, MM, SS? 

15,06, 10 

Solar  Decl  fins  =  ? 

19.2 

Scan  Line  Print?  < 1 > Yes/C0>No 
o 

i 


*  (loading  RUFF 3) 


(b)  WIND  SPEED < My S) : 1<0.0),2<5.0),3< 
7.0), 4(10. 0>? 

3 

WIND  RZ.Cie  N=060 ,  E=090,  S=180,  W= 
270)? 

000 


*  (loading  RADTRA) 


(c)  REL.  HUMIDITY, RH<5J)  :  1<0X)  ,  2<50JJ) 
,  3<70*>,4<80JO,5<90JO,6<95JO,7<9 
8*) , 8<99*>? 

6 

AEROSOL  SCALE  HTS . <KM> « FINE, COBR 
SE 

9 

!0S,0 

AEROSOL  MODEL ! < 1 )TROPOSPHERIC, <2 
)MARITIME, <  3 ) RURAL , <4>UR6AN? 

2 

*  (loading  required  data  files) 


<d>  ALL  FILES  READ 


Figure  B-3.  Display  during  simulation  run 


■ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft 
*OMSP  VHR/LF  Imagery  Simulation# 

*  R.G. ISAACS/AER, Inc./lSSl  * 

*  Version  Sl/S  * 

******************************** 

RADIANCES  FOR  SCAN  LINE  AT  * 


LONG <  AN  >  - 

94 

AN  <  GMT  >  = 

15 

0 

3 

Scan  ft!  GMT  >  — 

15 

6 

10 

Sol  Decl .  = 

19 . 2 

Sol  Lons.  = 

46 . 541666667 

Sub.  Pt .  Lat/Lons  = 
21 . 172148517 

/ 

98 . 3061061321 

Scan  Line  Pixel  = 

1 

Pixel  Lat/Long  = 

21 . 172148517 

/ 

38 . 3061061321 

Sat.  Sen  Ans  « 

0 

Hour  Ans .  = 
52.3644335251 

Sol  Zen  Ans  = 

48 . 3648006C.36 

Sol  Azm  Ans  = 
82.5050331655 

Scan  Line  Pixel  = 

2 

Pixel  Lat/Long  = 

21 . 3675008272 

/ 

97.434035738 7 

Sat.  Zen  Ans  = 

10 . 0680032343 

Hour  Ans.  = 

,  50.3523631317 

Sol  Zen  Ans  = 
47.6348009984 

Sol  Azm  Ans  = 
83.0424938691 

Scan  Line  Pixel  = 

3 

Pixel  Lat/Lons  = 

21 .5510145928 

/ 

96.0783079082 

Sa  t .  Zen  Ans  = 
19.550054325 

Hour  Ans  = 
43.5366412412 

Sol  Zen  Ans  = 
46.3048143283 

Sol  Azm  Ans  = 
83.5874013942 

Scan  Line  Pixel  = 

4 

Pixel  Lat/Lons  = 

21 . 7225478104 

/ 

34.6591113539 

Sa  t .  Zen  Ans  = 
28.0421331669 

Figure  B-4.  SLGP 

output. 

After  printing  the  heading  at  the  top  of  Figure  B-5, 
subprogram  RUFF 3  is  loaded.  When  this  operation  is  over, 
the  display  asks  for  wind  speed  and  wind  direction  informa¬ 
tion  (Figure  B-3b) .  After  the  wind  azimuth  is  entered,  the 
ocean  surface  model  calculations  are  performed  for  each  pixel. 
(The  display  screen  will  flash  with  various  intermediate  re¬ 
sults.)  Final  pixel-by-pixel  effective  reflectivities  are 
printed  as  demonstrated  in  Figure  B-5.  Completion  of  this 
operation  results  in  the  heading  in  Figure  b-6  being  printed 
followed  by  loading  of  the  radiative  transfer  program  (RADTRA) . 
This  portion  of  the  code  requests  (see  Figure  B-3c)  relative 
humidity  (only  options  2  thorugh  6  will  work) ,  fine  and  coarse 
aerosol  scale  heights,  and  one  of  the  four  available  aerosol 
optical  properties  models.  Subsequent  to  specifying  these 
input  parameters,  required  relavent  data  files  are  read. 

When  completed,  the  display  notes  "ALL  FILES  READ"  and  the 
radiative  transfer  calculation  for  each  pixel  commences. 

As  illustrated  in  Figure  B-6,  the  output  consists  of 
zenith  angle  cosines  for  sensor  (observer)  and  sun,  azimuth 
angle  difference,  and  effective  reflectivity  for  each  pixel, 
followed  by  a  radiance  spectrum  within  the  DMSP  bandpass  from 
0.4  to  1.1  ym  (each  0.02  ym)  consisting  of  first,  the  total 
radiance,  and,  second,  the  surface  contribution  only  for  com¬ 
parison  purposes.  Following  the  wavelength  dependent  radiances 
the  DMSP  response  weighted  radiance  (total  only)  is  given. 

Radiances  are  calculated  for  each  numbered  pixel  cor¬ 
responding  to  the  location  given  in  the  scan  line  print 


*********OCEAN  SURFACE******** 
REFLECTANCE 


Wind  Speed  =  3 

Wind  Di r . (dee)*  0 

.656522556696  1 

PIXEL  #*  1 

Rf (effective)* 

1 .22914923601E-4  POL.  = 

.290272498678 
.673853731358 
.984600959627 
PIXEL  #*  2 


Rf  (ef  f ect ive)  = 

1 . 74398210756E-3 
.407491354815 
.698821660688 
.942349513481 
PIXEL  #* 

Rf (ef  f ect ive)  = 

1 . 009 1 0703736E-2 
.531588891445 
.707417198008 
.882602122384 
PIXEL  #* 

Rf  (ef  f ect i ve>  = 

2 . 97354893974E-2 
.647079252209 
.72363139738 
.8153044345 
PIXEL  #= 

R f (ef  f ect i ve)  = 
5.43412865455E-2 
.743166854252 
.739455518437 
.747838184897 
PIXEL  #= 

Rf  <ef  f ect ive)  = 

7. 1 5977009 487E -2 
.816479854642 
.754831031085 
. 684405604278 
PIXEL  #* 
Rf(effective)* 

7 .55542157374E-2 
. 869162695817 
. 76989962009 
.626876843527 
PIXEL  4= 

Rf  <ef  f ect ive)= 

6 . 84923 1 38666E-2 
.905609943523 
.784503189581 
575668001212 


POL.  = 

3 

POL.  = 

4 

POL.  = 

5 

POL.  = 

6 

POL.  = 

7 

POL.  = 


POL.  = 


Figure  B-5.  RUFF3  output 


*******:URAOIANCES****:U**** 
REL.  HUM  CRH003=  6 

AEROSOL  SCALE  HTS<KM>- 


FINE  : 
rnopcc . 

MARITIME  AEROSOL 
PIXEL  *  = 

OBS  2ENI TH<v)  = 
SUN  ZENITH<po)  = 
.656522556697 
AZIMUTH  ANG  .  <$>  = 
-1.205939165 
Rf <e# f . >= 


1  22914923601E-4 
.4 

5.62023339111 
6 . 5921 4596227E-3 
.42 

6. 11269169159 
9 . 4228298362E-3 
.44 

5 . 16401864033 
1  . 10159262563E-2 
.46 

4 . 92361403792 
1 . 39160356865E-2 
.48 

3.99839627317 
1 . 46528283137E-2 
.5 

3.15216336054 
1  49552162204E-2 
.52 


2.4483416146 

1 . 47966128413E-2 

.54 

1 .96670604235 
1 . 47568724604E-2 
.56 

1 . 6391982723 

.014454008219 

.53 

1 . 38375212138 
1 . 47455989775E-2 

.6 

1 . 12580610198 
1 . 46729097 147E -2 
.62 

.968356891441 
1 . 45658123949E-2 
.64 

.824724520471 
1 . 45226991 842E-2 
.66 


.697253079743 
1 . 440031 7S787E-? 


.  68 

.584552465139 
1 .41825712095E-2 
.7 

.444821611869  8 

.72 

.359214328192  0 

.  74 

.334542917271  0 

.76 

.297924222246  0 

.78 

.259153684984  0 

.8 

.224477329539  0 

.82 

.14158028171  0 

.  84 

.158181618455  O 

.86 

.15594041316  0 

83 

.137305041092  0 

.9 

.100996069207  0 

.92 

4 . 3391S657656E-2  0 

.94 

2 . 66866606121E-2  0 

.  96 

3 . 1 4809637237E-2  0 

.98 

.057094337495  0 

1 

7 . 33573788027E-2  0 

1 . 02 

. 066103012938  0 

1  .04 

5 . 94410921 66 IE-2  0 

1 .06 

5 . 36809643037E-2  0 

1 . 08 

4 . 87745057828E-2  0 

1  .  1 

. 027649917289  0 


BANDPASS-WEIGHTED  INTENSITY 
618357806254 


Figure  B-6.  RADTRA  output  (pixel  #1  only) . 
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